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SUMMARY 


This report, an addendum to the six month report submitted 
to NASA Langley Resesearch Center in Dececember 1987, will cover 
the research performed by the Fiber and Electro-Optics Research 
Center (FEORC) at Virginia Tech for the NASA Langley Research 
Center Grant NAG-1-780, for the period from December 1987 - June 
1988. This final report will discuss the research performed in 
the following four areas that are mentioned in the proposal: 
Fabrication of Sensor Fibers Optimized for Embedding in Advanced 
Composites, Fabrication of Sensor Fiber with In-Line Splices and 
Evaluation via OTDR Methods, Modal Domain Optical Fiber Sensor 
Analysis, and Acoustic Fiber Waveguide Implementation. 


1 



DISCUSSION 


This section will discuss the four research areas mentioned 
in the NASA Langley Grant, NAG-l-780. The four research areas 
which will be discusssed are the following: Fabrication of 

Sensor Fibers Optimized for Embedding in Advanced Composite 
Materials, Fabrication of Sensor Fiber with In-Line Splices and 

i 

Evaluation via OTDR Methods, Modal Domain Optical Fiber Sensor t 

Analysis, and Acoustic Fiber Implementation. 

1. Fabrication of Sensor Fibers Optimized for Embedding in 

Advanced Composite Materials 

Optical fibers embedded in materials, i.e. composites, can 
sense position, strain, temperature, structural vibrations, and 
impact [2]. Several embedded optical fiber sensors have been 
investigated in the past six months. A Mach-Zehnder 
interferometer was used to sense tensile and composite strains in 
a composite material laminate [1], Acoustic emission in 
symmetric cross-ply composite laminates due to axial loading 
was interrogated with embedded optical fibers. This data was 
successfully compared to a piezoelectric transducer adhered to 
the composite where the transducer monitored acoustic emissions 
due to axial loading [3]. Fibers attached to an aluminum plate 
were used to assess impact damage. The various signals could be 
correlated to an impact height. Also, methods were used to 
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determine impact location [4]. The integrity of the composite 
due to embedded optical fibers has been examined by Moire 
interferometry techniques [5,11]. There has been a movement 
toward creating an 'intelligent structure.' An intelligent material 
consists of embedded actuators and sensors. Fiber optic sensors 
used to create intelligent aerospace structures have been examined 
[16]. Several papers review what smart structures are [6,7,13, 

17 ] . Problems associated with smart structures have been 
identified [8]. 

In the past six months there has been no progress in 
developing specialty fibers for improving specimen integretiy 
because the draw tower is not fully operational. 

2. Fabrication of Sensor Fiber with In-Line Splices and 

Evaluation via OTDR Methods 

Detecting localized strain with in-line air-gap splices in 
an optical fiber in conjunction with an Optical Time Domain 
Reflectometer (OTDR) was reviewed at the Smart Structures 
Workshop [9]. This experiment using the high resolution OTDR, 
has been delayed due to a defective OTDR. Dae to this delay, 
another graduate student working on this project, San jay Sudeora, 
will send a report to NASA as soon this difficulty is overcome. 

3 . Modal Domain Sensor Analysis 

* 

In the past six months the research in modal domain sensing 
has been directed towards analytically describing the movement 
of few-moded speckle patterns as perturbations are applied to an 
optical fiber [10,15,18], pursuing micro-optical fiber devices 
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which can be used with modal domain sensing systems [19] and 
real time signal processing of modal domain signals [14]. 

FEORC has analytically described the movement of two few- 
moded speckle patterns as perturbations are applied to the 
optical fiber: a three-lobed rotational speckle pattern and a 

two-lobed speckle pattern where the two lobes exchange power as 
the optical fiber is perturbed. The system sensitivity for both 
of these modal patterns has been described. FEORC has concluded 
that the two-lobed pattern is the optimum modal pattern 
[10,15,18] . 

In conjunction with modal domain sensing two in-line micro- 
optical devices can make a modal domain sensing system compact 
and packageable. An in-line laterally offset few-moded optical 
fiber splice can monitor one lobe of the two-lobed speckle 
pattern, thus performing in-line signal processing. The splice 
loss of two laterally offset few-moded optical fibers has been 
derived, because existing single mode and multimode optical fiber 
splice loss formulas can not be applied to few-moded optical 
fiber splices. The other device is a modal conversion coupler, 
which converts the input modes of the modal conversion coupler, 
a fused biconical taper coupler, to another set of modes. The 
amplitude coefficients associated with the modes propagating in a 
modal domain sensor are unknown. The modal conversion coupler 
may be a method of creating a mode with a known modal amplitude 
and superimposing several modes to create a totally 
characterizable modal domain sensor [19]. 
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A detector array can be used to process modal domain output 
patterns. This array Is linked to an IBM-PC to analyze the data 
in real time [14], 

Finally, the modal domain demonstration unit was delivered 
to NASA Langley in April 1988. This unit uses modal domain 
sensing techniques to detect quiscent strain, AC strain, 
structural vibrations, and acoustic emission in a graphite-epoxy 
composite . 

4. Acoustic Fiber Waveguide Implementation 

Acoustic waveguide sensors were reveiwed at the Smart 
Structures Workshop [11]. In the second six months of this 
contract, there has been no significant progress in the area of 
acoustic waveguides. FEORC is waiting for the arrival of some 
ultrasonics equipment and the draw tower to become fully 
operational . 
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APPENDIX A - Thesis Abstracts 



Modal Interference Techniques for Strain Detection in Few-Mode Optical 

Fibers 

by 

Bradley Dean Duncan 
Richard O. Claus. Chairman 
Electrical Engineering 
(ABSTRACT) 


Interference between the modes of an optical fiber results in specific intensity 
patterns which can be modulated as a function of disturbances in the optical 
fiber system. These modulation effects are a direct result of the difference in 
propagation constants of the constituent modes. In this presentation it is 
shown how the modulated intensity patterns created by the interference of 
specific mode groups in few-mode optical fibers (V < 5.0) can be used to 
detect strain. A detailed discussion of the modal phenomena responsible for 
the observed strain induced pattern modulation is given and it is shown that 
strain detection sensitivities on the order of 10"* can be expected Data taken 
during the evaluation of an actual experimental strain detection system based 
on the developed theory is also presented. 



Micro-Optical Fiber Devices l sed With Modal Domain .Sensing 

in 

Arris Flax 

R. O. Claus, Chairman 
Electrical Engineering 
(ABSTRACT) 

In order to develop more compact optical fiber sensing systems, modal filtering can be 
performed in-iine bv using micro-optical desiccs Tss'o such devices are a laterally offset 
fc’.v-moded optical fiber mechanical splice and a modal conversion optical fiber coupler. 
A third desia:, the air-gap splice used with multur.cde fibers, can examme the modal 
content of an optical fiber A basic thaorct.cil understanding on how these devices 
operate is rcvicsscd. A splice less calculation fer fesv-moded optical fibers is presented. 
Applications of the asymmetrical few-moded mechanical optical fiber splice, the modal 
conversion coupler, and the air-gap splice are discussed. 
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Micro-Optical Fiber Devices Used With Modal Domain Sensing 

by 

Amy Flax 

R O. Claus, Chairman 
Electrical Engineering 
(ABSTRACT) 

In order to develop more compact optical fiber sensing systems, modal filtering can be 
performed in-line by using micro-optical devices. Two such devices are a laterally offset 
few'-moded optical fiber mechanical splice and a modal conversion optical fiber coupler. 
A third device, the air-gap splice used with multimode fibers, can examine the modal 
content of an optical fiber A basic theoretical understanding on how these devices 
operate is reviewed. A splice loss calculation for few-moded optical fibers is presented. 
Applications of the asymmetrical few-moded mechanical optical fiber splice, the modal 
conversion coupler, and the air-gap splice are discussed. 
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1.0 Introduction 


Fiber optic sensing sj stems may be used to detect temperature, pressure, acoustic waves, 
strain, and a wide range of physical obsenabies The classical method for detecting 
these parameters with maximum sensiti\n> uses single mode fiber optic interferometers, 
which require two optical fibers. The reference fiber in such an interferometer must b» 
maintained in a controlled emironment for optimum detection. To eliminate this . 
of two optical fibers in a sensing system, the modal changes in one multimode optical 
fiber due to perturbations of the optical fiber can be exploited. This type of fiber optic 
sensing is called modal domain sensing. 

Modal domain sensing is most effeeme when used with a few-moded optical fibers. So 
far the few-moded optical fibers used with modal domain sensing support 2 to 10 modes 
[1,2] The optimum far field interference pattern of the fiber modes, the speckle pattern, 
is a modal field interference pattern that can be characterized by simple rotations or 
simple intensity changes m the far field intensity pattern due to modal changes from 
perturbations to the optical fiber [2] But the question now is w’hat to do with the optical 
information at the output of the sensor fiber (how to glean the relevant information 
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from the modal domain sensor). In-line micro-optical fiber devices would eliminate 
cumbersome non-optical spatial filters, realignment problems, and additional electrical 
circuitry which make far field pattern processing difficult. Micro-optical fiber devices 
create and enable compact, packageble modal domain sensing systems. 

Two such micro-optical fiber devices used with modal domain sensing systems are an 
optical fiber asymmetrical splice and an optical fiber modal conversion coupler. The 
optical fiber asymmetrical splice consists of a few-moded optical fiber spliced to a 
laterally offset few-moded optical fiber or to a single mode optical fiber for in-line spatial 
filtering. The optical fiber modal converter is an optical fiber coupler which generates a 
set of output modes from a given set of input modes. 

Another important micro-optical fiber device is the air-gap splice used with multimode 
optical fibers to detect local strain changes via Fresnel backscattered-based Optical Time 
Domain Reflectometry (OTDR) [3]. The air-gap splice causes pulse delay effects 
depending on the modes which are propagating in the multimode optical fiber. 

The principles of operation of the two devices used with modal domain sensing, the 
asymmetrical optical fiber splice and the optical fiber modal converter coupler and the 
principle of operation of the air-gap splice will be explained in Chapter 4 of this thesis. 
Finally, future applications of these micro-optical devices will be discussed m Chapte 
5. 


0 Introduction 


2 





2.0 Optical Fiber Theory 


The basic parameters that determine the nature of the electromagnetic fields which may 
propagate in optical fibers, and the basic mode theory which describes those fields need 
to be reviewed in order to understand the operation of the micro-optical fiber devices 
intended for sensor signal processing 


2.7 Basic Parameters Governing Optical Fibers 


An optical fiber typically is a cylindrical waveguide fabricated from a dielectric, silica, 
as shown in Figure 1 [4]. The center of the optical fiber, the core, has an index of 
refraction of The core has a radius, a Surrounding the core is the dadding whu 
has an index of refraction, rt^ which is slightly smaller than n, . 

The optical fiber waveguide will onl> support a discrete number of modes, having a 
unique propagation constant, /?. The propagation constants of the modes supported by 
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Frjnre !. An Opnctl Fiber |4) 


.0 Optical Fiber Theory 



the optical fiber must fall between the wave number of the core and the wave number 
of the cladding, Equation 2.1-1. 


n 2 k o <. $ <, n x k^ 


( 2 . 1 - 1 ) 


where 

= a • f ree space wavenumber, 
mkq = /i.t) v V,£, , wavenumber m the cladding, and 
n x k 2 = i x uiyj i.io E o ’ wavenumber in the core 

The amount of light that can be accepted by the optical fiber is determined by the 
numerical ipcrturc (AVI). The numerical aperture defines the acceptance cone for the 
light that will propagate m the core. The numencal aperture is defined in terms of the 
indices of refraction of the core and cladding as 

NA = (m, 2 - n 2 2 ) ] ‘ 2 = sin 0, (2. 1 - 2) 

where 6 is the cone of acceptance angle as shown in Figure 1. 

The V — number of an optical fiber determines how many propagating modes the optical 
fiber will support The V — number is defined m terms of the numencal aperture (NA), 
the optical fiber's core radius (a) and input source's wavelength of operation (2)’as 
indicatedin (2.1-3). 


V = 




(2.1-3) 
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By knowing the V — number for a few-moded optical fiber, the number of propagating 
modes can be determined analytically The results of the analysis are shown graphically 
in Figure 2 (51. Now, the propagation constant for each mode can be determined from 
the normal^ed propagation constant, b, where 


b 


m 0 ) 2 ~ "i 2 

2 2 
«, - n 2 


, and 


(2.1-4) 


p = - « 2 2 ) + « 2 2 ] i/2 . 


(2.1 -5) 


2.2 Optical Fiber Mode Theory 


After reviewing the parameters that govern the operation of an optical fiber, a rev 
of basic mode theory is essential in understanding the principles of operation oi 
micro-optical fiber devices. The electromagnetic field solutions in optical fibers are 
solutions to the wave equation for cylindrical waveguides, i.e. 



(2.2-1) 


In Equation 2 2-1, k, is the free space wavenumber and n is the index of refraction 1 - 
= core, n, = cladding) Equation 2.2-1 can be separated into six scalar equations k 
components e , , e t , e , , h„ f^, and h,. The Laplacian operator, V 2 , can be separated into 
a Laplacian operator transverse to the direction of propagation, V, 2 , and a Laplacian 
operator in the direction of propagation, V, 3 , i.e. 
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Figure 2. b vs. V -number |5] 
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(2 2 - 2 ) 


V 2 = V, : + v 2 \ 

If the electromagnetic field solutions are solved for the z -component of the 
electromagnetic field, the r and o components of the electromagnetic field can be found 
from the known r-components These relationships will be defined later, after the field 
solutions for 2 -components are established. 

The following anahsis is for the scalar field component, e, or /^, which will be redefined 
as Y, 1 e 

V, : Y + {k 0 V - fi 2 )* - 0 , (2.2-3) 

where /> is the propagation constant 


Y can be described as a function of R and <D, where R and <D can be separated, as shown 
in the following equations. 

Y = K(r)0>(d>), (2 2-4) 


+^0 = 0, and (2.2 - 5) 



where, 

q 2 = k 0 2 n 2 -fi\ (2.2-7) 

The solution to Equation 2 2-5 can be one of the following 
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T &m * 


.4, sm v<p 

cl) = j/?, cos vd> J, (2 2 - 8) 

.4, sin vd> + B, cos v<t> 

and the solution to Equation 2 2-6 inside the core is a Bessel function of the first kind, 
l e. 


R = q 2 >0{ie.n = n,). (2.2-9) 

The solution to Equation 2.2-6 outside the core is a modified Bessel function of the 
second kind, i e. 


R = BK,{q,r), q 1 < 0 (i « n = rt 2 ). (2.2 - 10) 

So, for the scalar electric field in the z-dire.tion inside the core (r < a), the solution f 
a mode is Equation 2.2-11. Similiar equations hold true for h,, i.e. 


A x sin v<£ 

^ = Bi cos vi}> J 

/J,sinv<£ + B,cosv$ 


( 2 . 2 - 11 ) 


where 


q 2 = U 2 - [k 0 2 n, 2 - /fja 2 . , (2.2-12) 

For the scalar electric field in the z-direction in the cladding region (r < a), the solution 
is given in Equation 2 2-13 Similar equations hold true for h,, Le. 
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. ![ sin \<p 

e, = *. cos vo ) 

zl 1 sm i<5 + Z?i cos v<£ 


(2.2-13) 


where 


q 1 = if 2 = (r - * 0 2 «2 2 K 


(2.2 - 14) 


It is important to note 
(Equation 2 2- 151. 


that the V — number can be also defined in terms of U and 


r 2 = u 2 + »' 2 


(2.2-15) 


Since .he Held solutions for Che .-components of .he electromagnetic held, can be 
determined, .he r and « components of .he electromagnetic fields can be found m 
Equations 2.2-16 thru 2.2-19. 


■M 

B5e. 

+ 

cr 

fid) ch } 

r d<*) ) 


' P d*, 
K T c<t> 

5h z 

j ( R 

8h z 

oc ^ 

7v 

dr 

r 6<j> )' 

J 

( P c/i z 

8e z 

q 2 

\ ' e<t> 

+ a>c £f 


(2.2-18) 


(2.2-19) 
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Lsually the modes of the optical fiber arc defined in terms of //£„ and £//.., where v 
is the order of the Bessel function and m is the root of the v order Bessel function. A 
specialized classification of these modes arc the linearly polarized ( LP) modes. These 
are the modes that wall be specified m conjunction with the micro-optical fiber dei.ces. 

Linearly polarized modes are electromagnetic fields that arc linearly polarized. LP 
modes need to propagate in fibers that are wealth guiding. Weak guidance is defined 
by Equation 2.2-20. 


A = 


2 2 

n \ ~ n 2 n \ ~ n 2 


2n, - 


(2.2 - 20) 


The linear polarization causes degeneracies among the EH„ and //£_ modes (6J. The 
LP,„ can be defined in terms of //£._ and £//.. by the following relationships' 



i = v 4- 1, 

(2.2-21) 

fori 

> 2, LP tm = HE^ Xjn , or EH^ Xjn , 

(2.2 - 22) 

fori — 

1, LP Xm = TE 3m , T.\f 0m ,or HE lm , and 

(2.2 - 23) 


fori = 0, LP Cm = HE lm . 

(2.2-2' 


The solution for the LP modes in the core yields the following: 




MUrla) icos i<f>\ 
W) ls,n ij 


and 


(2.2 - 25) 


H y = T, E x 


(2.2 - 26) 
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e 



where 


Y. 



(12 - 27) 


The solution for the LP modes m the cladding similarly yields the following: 


r K^Vrla) (coi*4>\ 


H y = n Ey 


(2.2 - 28) 
(2.2 - 29) 


where 


y - i^- 
r 2 - V n, ■ 


(12 - 30) 


The LP modes can be decomposed further into even modes denoted by LP ^ , and odd 
modes denoted by LP tm , . The odd modes are defined by the electric field component 
with the cos (/$) term, and the e\en modes are defined by the electric field comp one., 
with the sin(/<£) term. 


The basic parameters of optical fibers, the mode theory, and the understanding of LP 
modes will aid in understanding the prinicple of operation of micro-optical fiber de% ices 
for fiber sensor applications. 
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3.0 Modal Domain Sensing 


The nos el idea of usrrg a single optical fiber's modes mstead of an optical fiber 
interferometer for detecting perturbations (i.e. strain, \ibrations, etc.) was discovered 
bj Layton [7] and Kingsle\ [S]. Using the modes of an optical fiber to glean information 
about the optical fibers surroundings will be referred to as modal domain sensing. 
During the last few years at Virginia Tech a better understanding of the mechanisms 
behind modal domain sensing and practical implementations of modal domain sensing 
has been accomplished by Bennett [9], Ehrcnfeuchter [11, Srinivas [10], 
Shankaranarasanan [2], and Duncan [11]. The primary mechanism behind modal 
domain sensing is the change in the modes' propagation constants due to perturbations 
to the optical fiber. These changes in the propagation constants cause changes in the 
modal interference pattern which can be related to the applied perturbations to the 
optical fiber. Some applications of modal do main sensing thus fir are the use of 
embedded optical fibers m composites to detect acoustic emission [9], optical fibers to 
detect structural vibrations [1,11], and optical fibers to detect axial strain [2,10,11]. 
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A basic experimantal set-up of a modal domain sensing system is seen in Figure 3. The 
laser, the lens, and the optical fiber aligner can be considered to 'the mput optics" for 
the optical fiber. The optical fiber is the sensing region. At the end of the perturbed 
optical fiber is the signal processing (the filter) which may be implemented using an 
out-of-line spatial filter (no optical fibers are used ) or an in-line optical fiber as a spatial 
filter. Chapter 4 will look at some spatial filtering techniques which may be implemented 
with an optical fiber. After the spatial filtering section is the detection scheme, a 
photodiode detector feeds an oscilloscope in order to monitor the signal from the modal 
domain sensor. Then, a computer can be used m conjunction with this system for 
frequency analysis of the signal 

This chapter will examine the phase modulation m an optical fiber introduced by 
perturbations to the optical fiber, then, the optimum modal interference pattern for 
modal domain sensing will be reviewed [11]. 


3.1 Phase Modulation 


To understand the phase modulation of the an optical fiber due to perturbations, the 
phase modulation in one mode will be presented first (i.e. the change in the mode's 
propagation constant ) [12,13]. 

The transverse electric field solution is defined in Equation 3.1-1. 

E t = (3.1-1) 
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figure 3. V Modal Domain Sensing Experimental Set-Up 
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The phase of a mode is defined by Equation 3.1*2. 

4> = pz + (3.1-2) 

Here, z is defined as the distance along the as the length of the optical fiber, L . The 
phase modulation due to a perturbation is represented in Equation 3.1-3. 

A<0 = 0Az + zA0 (3.1-3) 

Equation 3.1-3 is valid for axial strain (e, ) or small acoustic waves (small vibrations ). 
A z can be expressed m terms of axial strain and the length of the optical fiber. Equation 
3 1-4. 


A 2 = t^L, for axial strain 


(3.1-4) 


0A z can be expressed in terms of an isotropic stress due to hydrostatic pressure, P, 
Poisson's Tatio, n, and Young's modulus, E, as indicated in Equation 3.1-5. 


0A 2 = fic 2 L 


0(1 - 7'j)LP 
E 


(31-5) 


The A0 term depends on a photoelastic effect and a dimensional effect. Equation 3.1-6 
U2.13]. 


where, 
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n = refractive index 

ji = longitudinal propagation constant 

D — diameter of the core. 

Equation 3 1-6 can be decomposed further by expanding the A n term. Equation 3.1-7 
[12,13]. 

A« = --y-Ml ~ (3-1-7) 

where, 

c z — axial strain 

p = Poissons ratio 

Pi\, P \2 — photoelastic constants. 


It is important to note that for small changes in the index of refraction, i.e. the difference 
between the core and the cladding indices of refraction is slight, an effective index of 
refraction, n ^ , can be considered to be equal to n, and to n 2 (Equation 3.1-8). 




(3.1-8) 


(3.1-9) 
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For the second term in Equation 3.1-6 changes in phase due to changes m the diameter, 
AD, can be expanded into the following 


AD = nc z D 


cfi 

Also, the term — — can be rewTitten as Equation 3. 1-11 
du 


IL = cb 

dD 2/5a 3 aD 


(3 1 - 10) 


(3.1-11) 


Howes er, the change in phase due to changes in the diameter is negligible compared to 

£0 

— — An, so the second term in Equation 3.1-6 can be neglected. 
an 


Finally, A d> due to changes in the propagation constant from external perturbations 
applied to the optical fiber can be rewantten as Equation 3.1-12 [12,13]. 


A<? 


= **Ly - 


p» 


CO - d)Pn - mAi3 


(3.1 - 12) 


3.2 Differential Phase Modulation 


For the few-moded optical Fibers and the micro-optical fiber devices that are used with 
modal domain sensing an understanding of th** differenual phase modulation between 
the propagating modes when external perturbations are applied to the optical fiber is 
essential. Since the optical fibers used in the experiments for modal domain sensing 
support tw o modes, differential phase modulation for two modes will be considered. The 
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reason why a mode combination of two modes is used with modal domain sensing will 
be explained in Section 3 3. 

The difference in phase under static conditions can be defined as </> u (Equation 3.2-1), 
where 1 is for the first mode and where 2 is for the second mode. 

<t> 12 = 4*1 ~ $2 ( 3 . 2 — 1 ) 

For applied perturbations to the few-modcd optical fibers Equation 3 2-1 becomes 
Equation 3 2-2. 

<f>l2 + = [A, + AAj] [$2 + A</>2] (3-2 2) 

The differential phase modulation is described by Equation 3 2-3. 

A<>i2 = [A, + AAi] — [A 2 + AA2] + A 12 (3 2 — 3) 

Equation 3 1-12 can be restated for differential phase modulation between two modes, 
and Equation 3.2-2 can be written in the form of Equation 3.1-12. 

AAn = AzAf? 12 (3.2 — 4) 


where, 


A/? - p t 2 . 

* 

P\2 = Pi ~ Pl> 

A^j 2 = £r^|^l2 — 2 ~ A*)/* 12 — (3.2 — 5) 
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This is the differential phase modulation between two modes. This differential phase 
modulation is important for the amplitude modulation of the electric Held intensity. The 
intensity wall be explained for the optimum modal pattern in the following section 


3.3 The Optimum Modal Pattern 


Recently, there has been a great deal of research to determine an easily characterizeable 
modal mte'fercnce pattern This modal interference pattern needs to be easy to 
generate. Changes in the modal interference pattern need to be correlated with external 
perturbations applied to the optical fiber. 

It has been proven b> Duncan [11] that the modes that produce the optimum 
mterference pattern is the x- polarized even LP tl and LP n modes. Figure 4 [11]. This 
modal interference pattern meets the two requirements mentioned earlier. fhis pattern 
is quite easy to generate, and as this pattern is perturbed the two lobes exchange power. 
This exchange in power between the two lobes can be correlated to external 
perturbations being applied to the optical fiber. A sensing system using this modal 
interference pattern can achieve a system sensitivity for strain of 10-’ [11] Changes in 
intensity for a modal domain sensing system using the LP„ and LP n modal interference 
pattern need to be reviewed. 
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Figure 4. Interference Pattern: LP. t and LP U Modes [II] 
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3.3.1 Intensity of the LPoi and LPn Modes 

The fields in the core of the optical fiber, the LP 0i and LP n modes, are defined 
respectively in Equations 3 3 1-1 and 3.3. 1-2. 

E x = AMU 0 rla) (3.3.1 _ j) 

Hy — YyJ Q {U Q rja) 

E x = cos 4> (3.3.1 - 2) 

Hy = A 2 Y ] J l (l\rla) cos <p 

The fields for the LP tl and LP n modes m the cladding are defined respectively by 
Equations 3. 3.1-3 and 3.3 1-4. 

E x = 't,A' 0 (l T 0 r/a) (3.3.1 — 3) 

Hy = Aj^lV.rja) 

E x ~ cos <i> (3.3.1 - 4) 

Hy — /l 2 y,A^(lK 0 r/c) cos 

Now, the total fields m the core can be described by Equation 3.3.1-5. 

E x = [A^J 0 { Uqt I a)e~ J ^ m ~ * c > + A 2 J x {U\rla) cos (3.3.1 - 5) 

Hy = (Y t EjS, 
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^ 0I and are random phase terms 


It is important to note at this point that the amplitudes of the fields A, and A„ are 
undetermined. At this point the amplitudes to the LP 0l and LP n can not be determined 
from the launching conditions. The work based on these amplitudes will be relative to 
the modes, because the modal amplitudes cannot be determined explicitly. 

Intensity for the electromagnetic fields is found in Equation 3. 3.1-6. 

/= = -i-K.lfJ 2 (3.3.1 — 6) 

B> substituting Equation 3.3. 1-5 into Equation 3.3. 1-6, Equation 3.3. 1 -7 is obtained. 

/ = I„{A t V 0 Hu,/la) + A 2 V, \l\rla\ cos <*>] 2 + 3 j _ 

A \A 2 J q { U Q rla)J l ( L\ rja) cos(A/Jz — A^)cos0} 


where. 



— $01 ~ $ll> 

B> examining Equation 3 3.1-7 the first term of this intensity equation is the intensity 
of the LP 0 . mode. Figure 5 [11], The second term in the intensity equation is the 
intensity of the LP n mode, Figure 6 |11]. The third term in the intensity equation is the 
interference of the LP tl and the LP U modes The third term is defined as the 'cross 
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term." The cross term is a cosine function as seen in Figure 7 when this term is a 
maximum [11]. 

Equation 3.3 1-7 can be rewritten as Equation 3.3.1-8 as a scries of functions dependent 
on r, <f>, and l e. 


/ = U/i(7) + / 2 ('.d>) + / 3 (r,0,z)}. 


(3.3.1 - 8) 


The intensity of the LP a: and LP U modes when I 2 — 0 is displayed in Figure 4. 


The assumption for the rest of this paper is that A 2 = 2 A,. This assumption will 

enhance the changes in intensity due to applied perturbations. If — is less than 1 then 

A i 

the LP ot term will dominate, and a single bnght spot will be present at the output. If 

~ is much greater than one, little power will exist in the LP at mode, and minimal 
A, 

interference between the two modes will occur as perturbations are applied to the optical 
fiber [11], 


3.3.2 Intensity Modulation 

In modal domain sensing the main focus is the behavior of Equation 3.3 1-7 as 
perturbations are applied to the optical fiber. The third term of Equation 3.3. 1-7 is 
affected by the perturbations to the optical fiber. When 7j = 0, the LP n and the LF n 
intensity are as displayed, respectively in Figures 5 and 6. As the optical fiber is 
perturbed to the point where /, is a maximum, there is an exchange of power between 
the two lobes as shown in Figure 8(11] In Figure S one lobe is a maximum while the 
other lobe is a minimum. 
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Figure 6. Intensity of the LP, , Mode (III 
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Figure 7. The Cross Term at a Maximum [II] 
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The cross term intensity, Itf, <*>,:), causes a sinsusoidal oscillation between the two 
lobes as the optical fiber is perturbed. This change m intensity can be related to an 
applied strain for an applied axial strain or for small amplitude sinusoidal strain. 


In the I 3 {r,4>,z) term the difTerention phase modulation between two modes from 
Section 3.2 can be applied to the A P term. The intensity equation, Equation 3.3. 1-8, can 
be rewritten as Equation 3.3.2-2 by apphing the following trigonometric identity: 

1 + cos It = 2[ cos x] 2 , (3 3.2-1) 


/=/*": l + cos(A - A*»), 


/ = 



A/Jz — A \!i 
2 



(3.3.2 - 2) 


/," consists of I, , 7,(r) and I 2 {r, A) from Equation 3.3. 1-8, and = 21,". 


The phase modulation between two modes from Section 3.2, Differential Phase 
Modulation, can be applied to the A fiz term. Az can be defined as a function of strain 
where 


Az = i^z, where (3 3.2 — 3) 

z = Lj, the length of the optical fiber, c M d 
t z — axial strain in the z — direction 

Equation 3.3.2-2 can be explained as a function of strain, i e. 
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/ = 


J 


(3.3.2 -4) 


) 1 . 





/ A/?(.r VjLf - A \l> 
C0 \ 2 


where x is optical fiber parameters related to the phase modulation, j.e. indices of 
refraction, photoclastic constants, Poisson’s ratio. 


A change in intensity. A/, can be directly related to a change in an applied axial strain 
(Equation 3.3 2-5). Assume that A^ = [1 1J. 


For axial strain- 


A/(A£;) 



/ 

-T“xA HL^ r 


(3.3.2 - 5) 


For small amplitude smsusoidal axial strain: 

peak.,) = ~-{x^pLj)^t 2 ptak COS(Ot. (3.3.2 - 6) 


A is the peak of the sinusoidal strain and co is the frequency of the sinusoidal strain. 
Equation 3.3.2-6 means that there are sinusoidal variations in intensity when the bending 
of a vibrat ng structure is quite small. 

Tkeorcticall} , Duncan has shown that the minimum detectable axail strain, A c„ is 
0 3.x 10-’ by a PIN diode and is 2.x 10“ by an Avalanche PIN diode (APD), which is the 
shot noise limit of the detector. 
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Two of the micro-optical fiber devices presented in this paper wall be discussed in 
conjunction with the ideas presented m this chapter, modal domain sensing 


3.0 Modal Domain Stiuui; 


31 



4.0 Principles of Operation of Micro-Optical Fiber 
Derices 


Three micro-optical fiber devices will be considered- the modal conversion coupler, the 
asymmetrical splice and an air-gap splice. Two of these devices, the modal conversion 
coupler and the asymmetrical splice can be used for signal processing in conjunction 
with modal do mam sensing The third micro-optical fiber device, the air-gap splice, 
detects modal effects in multimode step-index optical fibers. The air-gap splice is used 
with Optical Time Domain Reflectomctry (OTDR) techniques to detect localized 
changes in strain. 


4.1 Modal Conversion Coupler 

» 

A modal conversion coup'er is a fused biccnical taper coupler that is not made to split 
power m a 50 — 50 ratio, but uses the taper and the coupler region to generate new 
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modes at the one of the coupler's outpir arms. The modal conversion coupler that was 
fabricated in the Fiber and Elcctro-Opucs Lab by Kent Murphy is shown m Figure 9. 

A similiar deuce, a modal filter, has been fabricated by Sorin, Kim, and Shaw [14]. This 
device is shown in Figure 10. This device is an optical fiber coupler fabricated from a 
dual core optical fiber and a single mode optical fiber This device uses an applied 
acoustic wave to couple between the LP 0i and the LP n modes. This system can then 
generate the LP n mode from the LP 0l mode. However, this system is not entirely optical 
like the modal conversion cojpler and requires external perturbations to induce modal 
conversion. 

The all optical modal conversion oupler was fabneated from a 9/im core optical fiber, 
which is single mode at 1 m The modal conversion optical fiber coupler was excited 
at 633 nm using a Helium - Neon (HcNe) laser At 633/im the V— number for this optical 
fiber is 4.467. At this V — number this optical fiber supports the following modes: 
//£,„ HE*. HE a , £//,„ and IIE Jl [5]. This optical fiber also supports the degeneracies 
of these modes, the x- and v- polarizations of the these modes and the TE and TM 
modes. The modes that are launched into the optical fiber depend on the launching 
conditions. 

As seen in Figure 9, the EH tl and //£,, modes plus minor amounts of some other modes 
are launched into Arm 1 [15]. Arm 1-3 is one optical fiber, and Arm 2-4 is the other 
optical fiber, the coupled optical fiber. In Arm 3 another mode combination is 
generated. This optical fiber has tbc II mode and probably some of the HE n mode 
[15]. These two modes were deduced to be present experimentally by observing the 
concentricity of the far field pattern. Arm 4 has another mode combination, but the 
signal in this arm is about 10 dB down in power from that in Arm 3. Therefore, this 
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Fijwe 9. A Modal Conversion Coupler 
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coupler converts modes ir. the input arm, Arm I , to other modes in the output arm, Arm 

a 

How this coupler cor.serts modes is qualitatively described. In the coupler region, the 
l' -number becomes quite large, and new modes arc generated m this region. At both 
ends of the coupler is a taper region In this taper region higher modes are lost due to 
the decrease in the V — number. Due to the combination of new modes being generated 
in the coupler region and modes being lost in the t 2 pcr regions, some new modal 
combinations, which are dilTcrcnt than the input modal combinations, can propagate in 
the 9 iun optical fiber This modal conversion coupler a a method of creating other 
modes from a given set of input conditions 


4.2 An Asymmetrical Optical Fiber Splice 


This micro-optical device is used for signal processing in modal domain sensmg systems 
This splice is designed to work in conjunction with the LP m and LP n x-polarized even 
modal combinations. As strain is applied to the optical fiber, the two lobes exchange 
power. The asymmetrical splice is designed to collect the change in intensity of one lobe 
of the speckle pattern as shown m Figure 11. The intensity changes m the filtering 
optical fiber can be correlated to perturbations applied to the sensing optical fiber. 

One way this asymmetrical splice can be fabricated is shown m Figure 12. The filtering 
fiber is identical to the sensing fiber, but the filtering fiber is laterally offset with respect 
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to the sensing fiber. The offset filter fiber can detect the intensity changes of one lobe 
of the LP q{ and LP U medal pattern due to perturbations applied to the sensing fiber 

Another method of achieving this same effect is to have the filtering fiber's core diameter 
be one half of the sensing fiber's core diameter as in Figure 13. The filtering fiber's core 
then will sit inside of the sensing fiber's core. 


4.2.1 Analysis of Laterally Offset Optical Fiber Splices 

To anahzc this asjmmctncal splice problem for fevv-moded optical fibers, one may want 
to apply the standardized splice loss formulas for multimode optical fibers which are 
reviewed m Section -12 2 Alternative!}. one may want to apply standardized spbee loss 
formulas for single mode optical fibers which are reviewed m Section 4.2.3. This splice 
problem is neither in the single mode nor the multimode regime. A method for 
evaluating lateral splice losses for fevv-modes optical fibers is developed in Section 4.2.4. 


4.2.2 Laterally Offset Splices fer Multimode Optical Fibers 

For multimode optical fibers (F> 10) the anaylsis of the laterally offset splice is based 
upon the fact that many modes propagate in the sensing fiber and the filtering fiber. 
The number of modes that propagate in a multimode optical fiber can be approximated 
by Equation 4-2.3- 1 [5J 

V 2 

d/=-L-. (4.2 2-1) 
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Figure 12- Two Identical Fibers Laterally Offset 
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Sensing fiber's core 
radius is a s 


Filtering fiber's core 
radius is a , 


a* is at least 2 times 
greater than a f 


Figure 13. The Filtering Fiber's Core Inside The Sensing Fiber's Core 
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Lateral displacement between two optical fibers can be approximated as Equation 
4 2 2-2 [5] 



(4 2.2-2) 


where, 

a — core radius, 
d = lateral offset, 

A. 3mn = common core area bet.-.een the sensing and the filtering fiber 



This formula can also be applied to fro fibers of slightly different sizes except 
= -a, 2 , where a, is the cere radius of the filtering fiber. For the asymmetrical 
splice fer the LP zl and LP n modes, these formulas cannot be applied. The few-moded 
optical fiber is not operating m the true 'multimode regime". The V — number that 
supports the propagation of only the LP Zi and LP n modes falls between 2.405 and 3 9. 


4.23 Laterally Offset Splices for Singlemode Optical fibers 

For singlemode optical fibers ( V < 2.405) the propagating electromagnetic wa\e can be 
modeled as a Gaussian beam. The splice losses can be defined in terms of the beam 
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waist, w. The loss in a single mode splice due to a transiersc oflTsct can be defined by 
Equation 4.2 3-1 [16]. 


n = 4.34(4) m. 


(4.2.3 -1) 


where, 

d = trans\erse offset distance between the filtering fiber and the sensing fiber, and 
•.v = waist of the Gaussian beam. 


If two different optical fibers (,i.c. different core radii ) are used for the asymmetrical 
splice, mode field radius mismatch needs to be considered. The mode field radius is 
defined as ;,.e beam width of the Gaus^.an beam [16]. This splice loss is defined m 
Equation 4.2.3-2. 


n = 



(4.2.3 - 2) 


where, 

w r = waist of the Gaussian beam of the filtering fiber, and 
w R = waist of the Gaussian beam of the receding fiber. 

For the asymmetrical splice that supports the LP :i and LP n modes cannot be modeled 
as a Gaussian beam These modes do not fall in the single mode regime of operation 
(F< 2.405). Another method needs to be developed for few-moded optical fibers. This 
method will be presented in the following section. 
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4.2.4 Laterally Offset Splices for Few-Moded Optical fibers 


rhis problem needs to be solved b> working with the individual modes. The methods 
for splice loss evaluation developed for single mode or multimode fibers cannot be 
applied to few-moded fibers Another method for evaluation is presented in this section. 

For the few-moded sensing fiber, it is assumed that the LP a , and the LP {1 modes are 
propagating m the core 

= AJ 0 {U 0 n u)e~^ : - A.J.iU r/a) cos ve^' 2 ]^ (4.2 4-1) 

The fields that arc propagating m he sensing fiber's cladding are defined by Equation 
4.14-2. 

- [BM W t r!a)e~ J ^ - 5, A' (IF, r/a) cos <t>e~ jP 'a x . (4 2.4 - 2) 


The mterface of the splice is shown m Figure 14. Both fibers are identical in this 
configuration. At this point each fiber will be referred to its own coordinate system. 
The sensing fiber will be defined in terms of unprimed coordinates: r, 4 >, The filtering 
fiber will be defined in terms of primed coordinates: r\ </>', z' . 


A method of determining the amount of power which is received by the filtering fiber is 
to determine the modes propagating in the filtering fiber from the sensing fiber, and then 
an mtensitv in the filtering fiber can be determined. A ratio of intensity between the 
filtering fiber and the sensing fiber can be expressed as Equation 4.2.4-3. 


(v twwiPtr 
fZbtr 


(4 2.4 - 3) 
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This ratio is the efficient} of the splice. 


Another loss across the interface is reflection losses. 1 e. loss due the to index of 
refraction mismatch Since , the transmission across the interface will be 

approximated as 1 . 

Since weak guidance is assumed because the LP modes are dependent on the weak 
guidance condition, the modal amplitudes in the second fiber can be determined from the 
propagating modes in the sensing fiber [17] 

Y) cos [fVYdr'dV 

(4.2.4 - 4) 

\fj\r') co s(/£) 1 7 r'dYd<f>‘ 


i. — 


j -4 


where. 


m - 


A a J s (U/la),r<a ' 1 
B t K s {W s r}a),r> a\ 


(J f (U f r’ia), r < a "J 
lK f {P,}Y!a), r>a J 


S = the surface area to be integrated. 

In Equation 4.2.4~4 f(r ) cos (sd>) are the modes propagating in the sensing fiber. A, and 
B, are the amplitude coefficients associated with the modes in the sensing fiber. 
fJ^Y) are the modes propagating in the filtering fiber. A f is the resultant modal 

amplitude associated with iheftr") cos^') mode. The denominator of Equation 4.2.4-4 
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is a normalization factor of the J\r‘) cos(/</>') mode The area which is to be integrated 
o\er is the filtering fiber or the primed coordinates The unpnmcd coordinates need to 
be converted ;o primed coordinates in order to perform this integration Hus operation 
wall be performed later 

As seen in Figure 14, there arc three distinct areas of integration’ 

1 the sensing core - the filtering core interface, 

2 the sensing cladding - the filtering core interface, and 

3 the sensing cladding - the filtering cladding interface. 

It is important to note that each of these three integrations needs to be performed twice. 
Since the LP zl and the LP n modes are propagating in the sensing fiber, some of these 
modes will be coupled into the LP a and the LP n modes in the filtering fiber. So a modal 
amplitude integration needs to be performed for the LP ai and LP U modes m each of these 
three areas. These integrations are defined as fodows: 

1 The sensing core - filtering core interface 



C^jci^o(^c// a )^o(^r// a ) 


4 - 


'1/01 = A sU MU ]s rja) cos ^{Uoflaf/dr'dy} 

[ \J,(U 0 /Ia)\ 2 r'dr’d<t>- 


(4.2.4 - 5) 


* 
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^fgiin 


jj [.^ 1 V 0 (L/' 0 //a)J,( i\/la) cos <y + 

4n = 4. x J { (U^la)UL\,r'la) cos <ft cos <ft’jrVrV<y} 

J | J x (L\jf'la) cos <y 1 2 r'dr'd<t>’ 


2. The sensing cladding - filtering core interface 


{l 


[# A K c( t1, Cs r i a )M u o/la) + 


4’). = 4-11 n( , 4W fl M,( t V'7«) cos 


] |y 0 (l'- c //a)I 


{ 1 , 


C4oi *o( IVosrlaV^/la) cos 0' + 


4n ~ 


4rll^t(^ \s r l a )AiUijr'la) COS $ COS 

J |y,(t/,//a) COS 0' 1 2 r'dr'd<j>' 


3. The sensing cladding - filtering cladding interface 


l J s 


Bjn = 4n i 4( l V\s r iriK 0 { IVof'ld) cos ^r'dr'd^'} 

f IWVtf'latfr-dr'dV 

J s 


(4 2.4 - 6) 


(4.2.4 - 7) 


(4.2.4 - 8) 


(4 2.4 - 9) 
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C^,A- 0 (I F 0 /W(**V'/«) cos * 

s 

Byn - B JU cos 4> COS (4 2.4-10) 

J \KAn\f, a) cos <t>'\ 2 r' dr' d<*>' 


From these three regions, only the sensing core - filtering core interface and the sensing 
cladding - filtering cere interface contnbute to the filtering cere's fields. 

In order to sohe this problem the unprimed functions need to be converted to primed 
coordinates A method of conserting the unpnmed coordinates to the pruned 
coordinates is to apph the Summation Theorem for Bessel functions [18]. The 
summation theorem is the following 


S* Zj,mr ) - ) J k {md)Z, + k {mr')e^ ~ ” (4.2.4 - 1 1) 

k=—oo 


where, 

Z, = J,_ Y„ H, w , orH™ 


w here. 


J. is a Bessel Function of the first kind. 


Y. is a Bessel Function of the second kind. 


4 0 Principles of Operation of Micro-Optical Fiber De»kes 


48 



H , 01 is o Hankel Function of the first kind. 


II , 05 is a Hankel Function of the second kind. 

/. is a modified Bessel Function of the first kind. 

K,isa modified Bessel Function of the second kind, 
m is arbitrary. 

The relationship between the primed and the unprimed cooordinatcs is shown in Figure 
15. d is the distance that the sensing fiber's core is offset from the filtering fiber. The 
angles relate to one another from simple geometry. 

.> + £ = £' (4.2.4-12) 

^ = <y - (4 2.4- 13) 

Equation 4.2 4-1 1 can be rewritten in terms of primed and unprimed coordinates. 


= \ J t (^Z, +t (nr’)e /i( " fl (4 2.4-14) 

k =* — oo 

To convert JfU^la) into the pnmed coordinate system via Equation 4.2 4-14, the 
following relationships are true. 

» 

v = 0 
2 . = Jo 
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nt = L' 0 la 


In order to write Equation 4 2 4-14 in a more compact form, the following recurrence 
relationship will be cm eked. 


J { . k) = (-l)V* 


(4 2.4- 15) 


Now, J : ( L\rja) can be written m terms of the primed coordinates. 


J 0 (i 0 dja)J 0 (U 0 r’la) + 

oo 

J * {V * r!a) = 2 £ ( - 1 ) k UU,Jia)J k (iyia) cos(k<t>') 


(4 2 4- 16) 


= i 


To convert J,(l/>/a) cos p into pruned coordinates can be rewritten as Equation 
4 2.4-17. 

J.{L\rla) cos (4.2 4 - 17) 


J\lj\rla) cos m terms of the primed coordinates can be expressed as Equation 
4 2.4-15 


J^U^r/a) cos <}> 



J 0 {Uidla)J x {US la) cos <}>' + 

J k (L\dla)( -l) k LJ k+l (iyia) cos[(* + 1)*'] 
- 7* _,({/,»» cos[(A - l)d>']]} 


(4.2.4-18) 


K^lV^/a) needs to be expressed in terms of the primed coordinates. It is important to 
note the following relationships. 
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(4 24- 19) 


//,'V) = 


= /4M 


So Equation 4 2.4-14 can be rewritten in terms of the Hankcl of the first kind. 


WVodlaV^'Xnyia) + 


H 0 { '\lV 0 rla) = 


*7.<- 


1)V X ( 1 i\Jla)H K (,) (lV 0 r’la) cos(k<t>’) 


(4.2.4 - 20) 


k = 1 


where 

W 0 = jir 0 

By applying the relationships of Equation 4.2.4-19 to Equation 4.2.4-20, /4(lf' 0 r/'a) can 
be written m terms of the primed coordinates. 

I 0 (l V^la^V.yia) + 

W - 2 £ ( - D k I k {iy c dla)K k {lV 0 r'la ) cos (*<*>') (4 ‘ 2 ' 4 _ * 21) 

*= l 


Finally, A',(lF,/-/a) cos ^ needs to be expressed in terms of the primed coordinates. 


A, ( IF, r/a) cos £ 


/ 0 (lF,rf/ c )A,(lt>7a)cos <V -I- 



/.(IF^/fljCA^^lFjr'/aJcosK/: + 1)<J'] + 
A* _,(»»>•/<!) cos[(* - 1 )<£']]} 


(4.2.4 - 22) 


Now all of the unprimed coordinates arc expressed in terms of the primed coordinates. 
The integrations discussed earlier can be performed. However, it is important to point 
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out that this integration is still difficult to perform. The geometry of this problem as 
shown in Figure 14 docs not lend itself to a readily doable integration m polar 
coordinates. The symmetry is lost. This problem could be modified, so the splice loss 
due to a lateral offset could be evaluated. Instead of offsetting the two identical fibers, 
the filtering fiber could be a fiber that has a core diamenter half of the sensing fiber's 
core, as shown in I igure 16. This filtering fiber's core would sit m the core of the sensing 
fiber. This configuration would still pick-up one lobe of the LP n and LP U modal 
pattern. However, the filtering fiber will onl> support the LP n mode because due to the 
decrease m the filtering fiber's core size, the l' -number nil be halved The integrations 
presented earlier will only deal with the LP r amplitude coefficient. Finally, it will be 
assumed that the power contribution to the filtering fibers core will be due to the 
sensing fiber's core 

Note. The filtering fiber s core inset in the sensing fiber s core will be considered an 
approximation for the case where two identical fibers are laterally offset. 

The modified integration can be written as Equation 4.2.4-23. 


{L 




DO ~ 

2 Y, ( -\) k J k {U 0s dl a )J k {U c /:a) cc ,(*<*') + 


* = l 


'{/oi ~ 


.4 in y 0 (LV'2/a)[y c (L 0 ///aU I (t/, 1 i'';a) cos 4>’ -1 

oo 

YWJIa'A -l)'[f H1 (V/<>) cos[(* + l)d>'] 

l 

A-i(£i//£) cos[(A. - l)o’]]} 

J \J 0 (U 0 /2!a)\ 2 r’dr'<W 


k = 1 


(4 2.4- >3) 
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r ,<> - sensing fiber 

r\ <? - filtering fiber 


Figure 16. Cross-Section of the Asymmetrical Splice: The Filtering Fiber Is Half The Core Diameter 
of the Sensing Fiber 
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This equation can be simplified because the limits of integration for <t>' are 0 to 2?r. Any 
cos (n<t>) integrated from 0 to 2 .t, where n is an integer, is zero 


I'T Uoi MU 0 /2la)J 0 (U 0 ^la)J 0 (U 0s r'laydr'd4>' 

J o J o 

fa/2 , 

\J 0 (U c /2la)\ 2 r'dr'd4>' 

n Jn 


AjOi — 


(4 2.4 - 24) 


'o ■'o 


This can be further simplified to Equation 4.2.4-25. 


fa/2 

'^soi^oi^os^la) J J 0 {b S jr’2ja)J 0 (U 0s r' la)^ dr* 

Ajoi = 7^ (4.2 4-25) 

\J 0 (U 0/ r'2la)\ 2 r’dr' 


The integration m Equation 4.2.4-2_ can be earned out after knowing the dimensions 
of the filtering and sensing fibers. This will be done in Section 4.2.5, Implementation 
of the Asymmetncal Splice. After the integration is performed, Equation 4.2.4-25 will 
take the form of Equation 4 2.4-26, 


A/oi ~ A m X 


(4 2.4 - 26) 


where A' is a number. 


The field propagating in the filtering core is defined by Equation 4.2.4-27. 

E = A M XJ 0 (U 0 r‘2la)e- J ^ 2 a x (4.2.4-27) 
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The efficiency of this splice can be expressed in terms of the intensities of the modes 
propagating in the sensing liber, and the intensity of the LP ai mode propagating in the 
filtering fiber. The intensity, of the field propagating in the filtering core is 


hi - -T-lf-f’l - -j-Aw 2 X 2 J 0 \u 0 /2la). 
The intensity of the LP C , mode in the sensing fiber is 

lyOl = ~=T ^JOl 

The intensity of the LP U mode in the sensing fiber is 

h = ~^~A sU 2 7 , ~{U ls rla) cos 2 <t>. 


(4.2.4 - 28) 


(4 2 4 - 29) 


(4.2.4 - 30) 


In determining the efficiency of the splice, Equation 4 2 4-34, the intensity functions, 
Ija X , I M , and /,„ need to be mtegrated over their respectise fiber core areas. The filtering 
fiber's intensity function, /„„ is mtegrated o\er the primed coordinates as indicated in 
Equation 4.2.4-31. The sensing fiber's intensity functions, I M and /,„, are integrated over 
the unprimed coordinates as indicated in Equauons 4.2.4-32 and 4.2 4-33. 


(mi. t ~ 



Ml. i 



2 Jq 2 (U 0 /2la)rdr'd<y 


2 J 0 2 {V 0 /la)rdrd<i> 


(4.2.4 - 31) 


(4 2.4 - 32) 
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(4.2.4 - 33) 


An i 



V, J ( 6 |//fl)c os 


The efficiency of this splice. > 7 , is described in Equation 4 2.4-31. 


V 


^JOI + ^j|J 


(4.2.4-31) 


l n 

where tj ai — — , coupling efficiency of the LP t , mode m the sensing fiber to the 

*rtH 1 t 

LP, j in the filtering fiber, rj A = — , coupling efficiency of the LP n mode in the 

*slU» 


sensing fiber the LP ai m the filtering fiber, and A m and A, n are weighting factors. As 
seen in Equation 4 2.4-25, the only mode which contributes to the filtering fiber is the 
LP ZI mode Equation 4 2.4-26 can be rewritten as 


■^/di — ^ioi^oi + ^xii^n 


(4.2.4 - 35) 


where 


*n = 0. 


The coupling efficiency of the LP tl mode in the sensing fiber to LP n mode m the filtering 
fiber, >/,„ is zero. Equation 4.2.4-34 simplifies to 


n = 


•^jOl’l'O! 


*J0] 


+ A 


sit 


(4.2.4-36) 


This splice loss calculation can be used to approximate the case where two identical 
fibers are laterally offset, as " -11 as the case where the filtering fiber's core sits inside of 
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the sensing fiber's core 10 log of the efficiency, Equation 4 2.4-36, will yield the splice 
loss. However, the splice loss can not be explicitly evaluated because A and A, n are 
unknown. 


4.2.5 Implementation of the Laterally Offset Splice 

This section will examine a system which uses a few-moded laterally offset splice for 
in-line signal processing with a modal domain sensing system. This splice could be 
fabricated with two identical fibers laterally offset, but the simplified version of this 
splice wall be considered in this section. The splice will consist of a filtering fiber which 
is half the core diameter of the sensing fiber. Thus, the mam interest will be the 
core-core transfer of power. 

Parameters of the sensing fiber arc- 


a = 4 Sum 


XA = 0.1 


n, = 1.458 
n, = 1 4543 

*-• 7S0nm 
V = 3.6249 
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*T 

Mt 


=■ 1.9571182 
U u = 3.0220145 

=* 1 300 nm 

Parameters of the filtenng fiber are: 
a — 2 nm 
NA = 0.1 
n, = 1 458 
n x = 1 4543 

= 780nm 
1.611073 
l/y = 1.57061 

1 . 1 iff i- = 633 nm 

The filtering fiber could have been singlemode at 780nm, but the core diameter is Sum. 
That is more than half of the sensing fiber's radius. This optical fiber will filter more 
than one lobe of the speckle pattern. So it seems more feasible to use an optical fiber 
which is slightly less than half of the sensing fiber's diameter. This filtering fiber has a 
diameter of Afim. 
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Equation 4.2.4-25 can be rewritten with the parameters of the optical fiber. 


rlitm 

A^Ul’oJIa,) J J 0 {U Q /la f )J 0 (U Q /layd/ 

j^r ; H-2-5 - 1) 

\J Q {V ¥ 'la f )\\'dr' 


where 

d = 2 0 (un to2 5 nm, 
a, = 2um, 
a, = 4 5/xm, 

L'„, = 1 9571182. and 
U v = 1 57061. 

After integrating Equation 4.2 5-1, can be expressed as Equation 4.2.5-2, 

A m = A m X. (4.2 5-2) 

where X is calculated from the program. Amp Fortran, which is found in the Appendix 

A. Then »j ei . Equation 4 2.5-3, can be calculated after integrating /«, and / rtl These 

* 

calculations are performed by using the fortran program, Int Fortran which is found in 
Appendix B. 
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*1oi — 


(4.2.5 - 3) 


Y 4o> 2 x 2 Jo \u Q /ia f ydr’d4>' 

-J - ^01 2y o 2 {U 0s rja t )rdrd<t> 

J 

Finally, the efFtciency of the splice can be determined from Equation 4.2.4-36. The 
assumption which was made m Chapter 3, that A, l{ = 2<4 <fl „ will be used to crudely 
examine the splice efficiency Equation 4.2 4-36 can be rewritten as Equation 4.2.5-4. 



n = 


JTO! 

3 


(4.2.5 - 4) 


The following \alues for X, and L , the splice loss, were calculated for various lateral 
offsets, d: 

d = 2.50/un: X = 0.8847932, tjm = 0.2118394, L = 11.51 

d = 2.25/im: X = 0.944511, if. = 0.2413998, L = 10.94 dB 

d = 2-00/un: JIT = 0 9994462, = 0.2676080, L = 10.496 dB. 

The filtering fiber will detect more of the sensing fiber's LP n mode, as the filtering fiber 
is brought closer to the center of the sensing fiber ( d goes to 2 fim). This is due to the 
fact that LP tl mode is a Gaussian beam. The optimum placement of the filtering fiber 
for power considerations is a lateral offset of d = 2.0/urt. Now, the propagating field in 
the filtering fiber s core for d = 2/im can be written as Equation 4.2.5-3. 
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(4 2.5 - 5) 


E = ^ lOI (O.9994462)J 0 (t; 0 //a / )a Jt 

In order to perform the splice calculation exactly, the amplitude coefficients in Equation 
4.2.5 need to be determined. This splice will enable in-line signal processing of the 
perturbations applied to the sensing fiber. 


4.3 The Air-Gap Splice 


A third micro-optical fiber deuce that is presented is the air-gap splice as shown in 
Figure 17 [3J. Unlike the last two deuces, the modal conversion coupler and the 
as\ mmctrical splice, the air-gap splice is not used in conjunction with modal domain 
sensing. The air-gap splice is used with an OTDR system to monitor localized strain 
effects. Figure 18 (3] The fiber which is used with the air-gap splice has a 109pm core 
diameter and 140pm cladding diameter, this fiber has a step-index profile. 

There is a Fresnel reflection at each of this air-gap splices, a pulse (4%reflected power) 
can be seen on the OTDR due to the air-gap splice. Due to the change in the medium 
at the air-gap splice (from glass to air ), some of the higher order modes will be lost as 
fight propagates in the air regime. The light will diverge, and the modes which travel 
closest to the core-cladding interface, the high'er order modes will be lost as cladding 
modes. 
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Figure 17. TV Air-Gap Splice |3| 
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The air-gap splice will magnify modal effects such as the modal content of the optical 
fiber and modal dispersion Modal effects due to modal dispersion are shown m Figure 
19 (3). Modal dispersion can be expressed as Equation 4.3-1 [5] 


1 mod 


= r _ r 

* rrva* * m 


n,AL 


(4.3-1) 


where 


T moe = modal dispersion in time , 


= path of the highest order mode. 


T'mm = pat/i of the lowest order node. 


n j — index of refraction, core, 



L = length of the fiber, and 
c = 3 x 10* m/s. 


is dependent upon the launching conditions, as shown in Figure 19 [3j. If mainly 
lower order modes are launched as in a of Figure 19, the pulse is not as very board. If 
higher order modes are launched into the optical fiber as in b of Figure 19, the pulse 
tends to be more boarden than tl* case where mainly lower order modes are launched 
into the optical fiber. This is true due to the later arrival times of the higher order 
modes. 
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Figure 19. Modal Dispersion Due to Launching renditions |3|: a) Lower Order Modes b) Higher 
Order Modes 
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As the fiber is strained, the pulses become boardencd and elongated as in Figure 20 [3J. 
This is due to modal dispersion As the modes travel down the fiber as the fiber is 
elongated (strained), the pulses will boarden. This effect is demonstrated tn Equation 
4 3-2 


T'mod ~ 


/i|A[L + AL] 


(4.3-2) 


where, 

A L = change in length of the optical fiber due to an applied axial strain 


The same effect can be seen with a fiber w hich has several air-gap splices. The further 
away from the beginning of the fiber each air-gap splice is located, the more boardened 
the pulse becomes. This boardenmg is due to modal dispersion. This is related to L in 
the numerator of Equation 4 3-1. The longer the light propagates the more dispersed it 
becomes. 
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5.0 Future Work and Applications 


These three micro-optical fiber devices, the modal cons ersion coupler, the asymmetrical 
splice and the air-gap splice, can help further research in fiber optic sensing and 
communications. 

In conjunction with modal domain sensing, the modal conversion coupler is a method 
of generating new modes from a given set of input conditions. The modal conversion 
coupler needs to be mathematically described. The mechanisms behind the modal 
conversions in the taper regions and the coupler region needs to be modeled. Once a 
model is established for the modal conversion coupler, a repeatable method for 
fabricating these couplers needs to be established. This method must be consistent to 
generate the same output modes from a given set of input modes. 

An application of this device is that a modal conversion coupler could generate the 
LP U mode. The LP lt mode amplitude could be measured from 'he output of the modal 
conversion coupler. This field could be superimposed with a known LP n field. With the 
correct modal amplitudes these two superimposed modes could generate the two-lobed 
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speckle pattern used for modal domain sensing. This would solve the present problem 
of unknown modal amplitudes due to variable launching conditions. 

Another application of this device is a multiplexer in a communications system. This 
device would use the method of modal multiplexing. Certain modal combinations could 
carry a 'channel of data " The coupler could be used lor generating the modal 
combination to be multiplexed onto a fiber. All of these modal channels would need to 
be superimposed onto the data carrying optical fiber. An analogous device would need 
to be fabricated to demultiplex the modal channels. This modal conversion coupler 
could increase the information earn ing capacity of a multimode optical fiber. 

Also, in conjunction with modal domain sensing, the asymmetrical splice is a means of 
performing in-line signal processing In order for this splice to be completely effective, 
a function correlating the intensity changes in the filtering fiber due to applied 
perturbations to the sensing fiber need to be derived. This function will be nonlinear, 
since the intensity is a point function, which is a linear function. This asymmetrical 
splice will allow for a compact and packageable modal domain sensing system. This 
splice will be in the near field as opposed to the far field, so more power can be coupled 
into this splice as opposed to fiber placed in the far field. In terms of the loss 
calculations associated with the asymmetrical splice, the modal amplitudes for the 
sensing fields in the core need to be determined, so this calculation can be complete As 
mentioned earlier, to solve the problem of unknown modal amplitudes, the modal 
conversion coupler could be applied for a completely characterizable modal domain 
sensing system. Once the modal amplitudes are known, the theoretical results can be 
compared to an experimental asymmetrical splice with the filtering fiber's core inside the 
sensing fiber's core as in Figure 16. Then, the theoretical approximation of the filtering 
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fiber's core inside the sensing fiber's core can then be compared to the experimental 
asjmmetrical splice with two identical fibers laterally offset as in Figure 15. 


Finally, the air-gap splice exhibits modal effects evident m multimode optical fibers. For 
applications where modal dispersion effects are undcsireable, the air-gap splice could be 
fabricated with single mode optical fibers. However, the air-gap splice could be used to 
monitor modal dispersal effects as a fiber is strained, or monitor disperson along the 
length of the optical fiber. The air-gap splice could be used to analyze the modal content 
which is launched into the optical fiber. This could be a method of determining whether 
lower, higher, or evenh distributed modes are propagating in the optical fiber. 
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6.0 Conclusions 


Three micro-optical fiber deuces, a modal conversion coupler, an asymmetrical splice, 
and an air-gap splice, were presented. Two of these devices can be used in conjunction 
with modal domain sensing, the modal conversion coupler and the asymmetrical splice. 

The modal conversion conversion coupler can be a method of describing launching 
conditions. Modal amplitudes in conjunction with modal domain sensing are unknown 
at this time. The modal cons ersion coupler can generate a mode from a given set of 
input conditions. The modal amplitude from the output of the coupler can be 
determined, and then superimposed with other known modal amplitudes. This is a 
possible solution to unknown launching conditions. 

The asymmetrical splice is a method of performing in-line signal processing The 
asymmetrical splice allows a modal domain sensing system to be a compact and 
packageable system. In ten sits changes in the filtering fiber can be correlated to 
perturbations in the sensing fiber. Finally, splice loss calculation! for laterally ofTset 
multimode and single mode fibers is not applicable to few-moded fibers. A new method 
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for evaluating splice loss for few-modcd fibers has been developed In order to determine 
the splice loss explicitly, amplitude coefficients for the modes propagating in the filtering 
fiber need to be determined from the modes propagating in the sensing fiber. 

The air-gap splices exhibit modal effects when used with multimode fibers. The air-gap 
splice is a method of viewing modal dispersion, whether modal dispersion is dependent 
upon the modal content launched into the optical fiber or is due to ler gth elongation 
due to axial strain applied to the optical fiber 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C AMP FORTRAN C 

c c 

c c 

c c 

C WRITTEN BY AMY FLAX C 

C C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C THIS PROGRAM IS DESIGNED TO AID IN CALCULATING AFOl IN C 

C C 

C EQUATION A 2.5-1, WHERE AFOl = ASOl X. THIS PROGRAM CALCULATES C 

c c 

C X. THE MAIN PROGRAM USES AN IMSL ROUTINE, DCADRE, TO DO THE C 

c c 

C INTEGRATIONS. DCADRE USES THE FOLLOWING PARAMETERS: C 

c c 

c F= FUNCTION TO BE INTEGRATED, WHICH IS CALLED EXTERNALLY C 

C A=LOUER LIMIT OF INTEGRATION C 

C B=UPPER LIMIT OF INTEGRATION C 

C AERR=INPUT ERROR C 

C RERR=RETURN ERROR C 

C ERROR=ERROR MESSAGES C 

C C=INTEGRATED ANSWER C 

C C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c c 

C DCADRE IS USED TWICE. C 

c c 

C 1) TO INTEGRATE THE NUMERATOR OF EQUATION 4. 2. 5-1. C 

C THE EXTERNAL FUNCTION TO BE INTEGRATED IS 'F*. C 

C THE RESULT IS 'C' . C 

C C 

C 2) TO INTEGRATE THE DENOMINATOR OF EQUATION 4. 2. 5-1. C 

C THE EXTERNAL FUNCTION TO BE INTEGRATED IS 'FX'. C 

C THE RESULT IS 'D*. C 

C C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


INTEGER IER 

REAL DCADRE, F, A, B, AERR, RERR, ERROR. C,D,FX, J.DI.UF.AF 
EXTERNAL F.FX 
A=0 . 0 
B=2.0 
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RERR=0. 

AERR= 1.0E-5 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C DCAORE IS CALLED TO INTEGRATE F AND FX. WHICH ARE FOUND IN C 

c c 

C SUBROUTINES AT THE BOTTOM OF THE PROGRAM. C 

C C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


C = DCADRE (F, A, B> AERR. RERR. ERROR. IER ) 
D = DCADRE ( FX. A, B, AERR, RERR. ERROR. IER) 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c c 

C THIS SECTION CALCULATES JOC UOS DI/AS) OF EQUATION A. 2. 5-1. C 
C C 
C THIS TERM ACCOUNTS FOR THE LATERAL OFFSET Or THE FILTERING C 
C C 
C FIBER. C 
C C 
C US = UOS/AS, FIBER PARAMETER CALCULATED C 
C c 
C DI = DISTANCE CENTER OF FILTERING FIBER OFFSET FROM C 
C CENTER OF SENSING FIBER C 
C C 
C J = JOS (US DI >, BESSEL FUNCTION OF THE FIRST KIND C 
C APPROXIMATED BY A POLYNOMIAL. THE APPROXIMATION C 
C OF THIS FUNCTION IS FOUND IN REFERENCE C 
C C 


CCCCCCCCCCC^CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


DI=2.00 

US=0 . A3A915159 

J=C1- 2.2A99997* ((US*DI/3.0)X*2 ) + 1 .2656208*( (US*DI/3 . 0>X*A) 
C - 0. 3 16 3866 X( (US*DI/3 0)**6) + 0 . 0A4AA79*C ( USKDI/3 . 0) *«S) - 
C0.00394AAK(CUS*DI/3.0)**10) - 0.0002100X(CUS*DI/3.0)*K12>) 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c c 
C X OF EQUATION A. 2. 5-1 IS CALCULATED. WHERE AF = X. C 
C C 
C X IS PRINTED OUT ON THE TERMINAL SCREEN. C 
C C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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AF=J*C/D 

WRITE (*.*) AF 

STOP 

END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C THIS IS A SUB ROUTINE TO CALCULATE THE FUNCTION IN THE 
C 

C NUMERATOR OF EQUATION 9.I.-5, (JO (UOS X / AS ) ) (JO (UOF X /AF)). 

C 

C 1) US = UOS/AS 

C 2) UF = UOF/AF 

C 3) THE BESSELS FUNCTION ARE APPROXIMATED BY A POLYNOMIAL 

C 9) S = JO (US X) 

C 5) FF= JO (UF X) 

C 6) F*S*FF*X, WHICH IS THE FUNCTION TO BE INTEGRATED 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


REAL FUNCTION F(X) 

REAL X-UF.US.S.FF 
UF=0. 7853052 
US=0 . 939915159 

FF=( 1- 2 2999997* ( (UF*X/3 -0)**2 ) + 1 . 2656208*( (UFXX/3. 0 )**9) 
C - 0.3163866*((UF*X/3.0)**6) ♦ 0 . 0999979*( (UF*X/3 . 0 )**8 ) - 
C0.0039999*((UF*X/3.0)**10) - 0 . 0002100*( ( UF*X/3 . 0)**12) ) 

S=(l- 2.2999997* ((US*X/3 0)**2 ) + 1 . 2656208*( ( US*X/3 . 0 )**9) 
C-0.3163866*((US*X/3.0)**6) + 0 . 0999979*( ( US*X/3 . 0 )**8 ) - 
CO . 0 0 39449* ( (US*X/ 3.0)**10) - 0.0002100*((US*X/3.0)**12)) 
F=5*FF*X 

RETURN 

END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C C 

C THIS SUBROUTINE CALCULATES THE FUNCTION IN THE DENOMINATOR C 

C C 

C IN EQUATION 9. 2. 5-1. THE METHOD IS THE SAME AS IN THE OTHER C 
C C 

C SUBROUTINE, EXCEPT THESE MODIFICATIONS: C 

C C 

C 1) FF. UF ARE THE SAME AS IN THE OTHER SUBROUTINE C 

C 2) THE FUNCTION TO BE INTEGRATED IS FX C 

C 3) FX = (ABS (FF) )**2 X C 

C C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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REAL FUNCTION FXCX) 

REAL X.UF.FF 
UF=0. 7853052 

FF=C1- 2 . 2499997* ((UF*X/3.0)**2 ) + 1 .2656208*((UF*X/3.0>**4> 
C - 0.3163866*CCUF*X/3.0)**6) ♦ 0.0444479»((UF*X/3.0>**8) - 
CO. 003994 A*( (UFXX/3. 0)**10) - 0.0002100*(<UF*X/3.0)X*12)) 

FX= (C ABS (FF) )**2 )*X 

RETURN 

END 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 

c 

c 

c 

c 

c 


I NT FORTRAN 


WRITTEN BY ANY FLAX 


C 

C 

C 

C 

C 

C 

C 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 

c 

c 

c 

r 

w 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS PROGRAM IS DESIGNED TO AID IN CALCULATING NOl IN 
EQUATION A. 2. 4-31. THIS PROGRAM CALCULATES 

NOl. THE MAIN PROGRAM USES AN IMSL ROUTINE. DCADRE. TO DO THE 
INTEGRATIONS DCADRE USES THE FOLLOWING PARAMETERS: 


C 
C 
C 
C 

c 

C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


F=FUNCTION TO BE 
A=LOWER LIMIT OF 
B=UPPER LIMIT OF 
AERR=IHPUT ERROR 
RERR=RETURN ERROR 
ERROR=ERROR MESSAGES 
C=INTEGRATED ANSWER 


INTEGRATED, WHICH 

INTEGRATION 

INTEGRATION 


IS CALLED EXTERNALLY 


C 

c 

C 

C 

c 

c 

c 

c 

c 

c 

c 


DCADRE IS USED TWICE. 

1) TO INTEGRATE THE NUMERATOR OF EQUATION 4.2.4-31. 

THE EXTERNAL FUNCTION TO BE INTEGRATED IS 'F*. 

THE RESULT IS 'C' WHICH IS INTEGRATED FROM 0 TO 2. 

2) TO INTEGRATE THE DENOMINATOR OF EQUATION 4. 2. 4-3. 

THE EXTERNAL FUNCTION TO BE INTEGRATED IS »FX*. 

THE RESULT IS 'D' WHICH IS INTEGRATED FROM 0 TO 4.5. 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


INTEGER IER 

REAL DCADRE, F, A, B, AERR, RERR, ERROR, C, D, FX. E, XI . A2, B2 
EXTERNAL F,FX 
A = 0 . 0 

B-2 . 0 
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RERR=0 . 

AERR= 1.0E-5 
A2 =0.0 
B2= A. 5 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c c 

C DCADRE IS CALLED TO INTEGRATE F AND FX, WHICH ARE FOUND IN C 
C C 
C SUBROUTINES AT THE BOTTOM OF THE PROGRAM. C 
C C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C = DCADRE (F,A.B,AERR»RERR,ERROR»IER > 

D = DCADRE CFX,A2,B2.AERR»RERR, ERROR, IER) 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C c 
C X=X1, CALCULATED FROM THE AMP FORTRAN PROGRAM IS ENTERED HERE. C 
C C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

XI = 0.994462 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C C 
C NOl OF EQUATION 4.2.4-31 IS CALCULATED, WHERE NOl * IFOl/ISOl C 
C C 
C E = NOl C 
c c 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


E=CX1 ** 2 ) * C /D 


WRITE (*,*) E 

STOP 

END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C C 
C THIS IS A SUBROUTINE TO CALCULATE THE FUNCTION IN THE C 
C C 
C NUMERATOR OF EQUATION A. 2.4-31, JO C UOF X /AF). C 
C C 
C 1) UF = UOF/AF C 
C 2) THE BESSELS FUNCTION ARE APPROXIMATED BY A POLYNOMIAL 1193 C 
C 3) FF= JO (UF X) C 
C A) F=((FF)**2)*X, WHICH IS THE FUNCTION TO BE INTEGRATED C 
C C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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REAL FUNCTION F(X) 

REAL X.UF.FF 
UF=0. 7853052 

FF=(1- 2.2499997* < (UF*X/3 . 0 )**2 ) + i .2656208*UUF*X/3.0)**4) 
C - 0.3163866*((UF*X/'3.0)**6) ♦ 0 . 0444479*( ( UF*X/3 . 0 )**8 ) - 
CO 0039444*1 (UF*X./ 3 . 0 )**10) ♦ 0 . 0002100XC (UF*X/3 . 07**12) ) 

F= ( ( FF)**2 )*X 

RETURN 

END 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS SUBROUTINE CALCULATES THE FUNCTION IN THE DENOMINATOR 
IN EQUATION 4.2.4-31. THE HETHOD IS THE SAME AS IN THE OTHER 
SUBROUTINE, EXCEPT THESE MODIFICATIONS: 

1) S AND US, WHERE US = UOS/AS. REPLACE FF AND UF 

2) THE FUNCTION TO BE INTEGRATED IS FX 

3) FX = (S)**2 * X 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


REAL FUNCTION FX(X) 

REAL X,US,S 
US=0. 434915159 

S=tl- 2.2499997* (IUS*Xy3.0)**2 > + 1 . 2656208*1 (USXX/3. 0>**4) 
C - 0 . 3163S66*( (US*X/3 . 0 )**6 ) + 0.0444479*( (US*X/3.0)**8) - 
C0.0039444*((US*X/3.0)**10) + 0.0002100*(CUS*X/3.0)**12)) 

FX= ( <S)**2 )*X 

RETURN 

END 
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(ABSTRACT) 


Interference between the modes of an optical fiber results in specific intensity 
patterns which can be modulated as a function of disturbances in the optical 
fiber system These modulation effects are a direct result of the difference i>. 
propagation constants of the constituent modes. In this presentation it is 
shown how the modulated intensity patterns created by the interference of 
specific mode groups in few-mode opt 1 si fibers (V < 5.0) can be used to 
detect strain. A detailed discussion of the modal phenomena responsible for 
the observed strain induced pattern modulation is given and it is shown that 
strain detection sensitivities on the order of 10 -9 can be expected. Data taken 
during the evaluation of an actual experimental strain detection system based 
on the developed theory is also presented. 
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1.0 INTRODUCTION 


The past several years have given rise to many optical fiber sensing 
techniques which may be employed in systems to detect such physical 
observables as temperature, pressure, displacement and strain [1], As 
sensing devices, optical fibers are attractive in part due to their small size, 
flexibility, ruggedness and intnns'c dielectric nature. In addition, their 
proclivity to be highly sensitive to external perturbations is quickly making 
optical fibers the sensing devices of choice for many passive non-intrusive 
sensor applications. 

Simply put. optical fiber sensors exploit the effects of perturbations and 
external factors on the light in the fiber. Sj»*jcifically, the intensity, phase and 
polarization are among those properties of light which are available for 
alteration and the detection of these changes is the basis of most fiber optic 
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sensors. Unparalleled in sensitivity are the phase modulated sensors, a 
classic example of which is the dual-fiber Mach-Zehnder interferometer [2], 
Because it is based on an absolute phase modulation scheme, however, such 
a sensor tends to be quite complex, with the reference phase and quadrature 
point often needing to be actively controlled. 

By comparison, sensor systems utilizing single-fiber modal interference 
schemes, such as those to be discussed in this thesis, essentially have the two 
arms of the interferometer within the same environment. This provides for a 
high common mode rejection capacity while affording the system the further 
advantages of simplicity and increased ruggedness, with no need for couplers. 
Though remaining highly sensitive, as will be shown, some sensitivity will be 
sacrificed in a single-fiber interferometer, or modal domain sensor, due to the 
fact that the sensor mechanism will now be based on differential phase 
modulation. 

Previous work with modal domain sensors has involved the detection of such 
observables as quasi-static strain [3,4], vibration [5], and acoustic emission [6]. 
Often, though, this work was done without the prior benefit of a rigorous 
theoretical basis for interpreting the sensor output. Also, no.effort as yet has 
been made to theoretically predict and experimentally verify the ultimate 
sensitivities achieveable with given modal domain sensor systems. Finally, 
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no attempt has yet been made to use anything but large discrete optical 
components in a modal domain sensor system design. This last point will, of 
course, limit the practicality of single-fiber modal interference based systems. 

The objectives of this thesis are therefore several. First the theoretical 
considerations of the behavior of the intensity patterns caused by the 
interference of specific mode groups in few-mode optical fibers (V < 5 0) will 
be thoroughly discussed and it will be shown how these intensity patterns, 
when modulated as a function of strain, can be used as a sensor mechanism. 
Secondly, the ultimate sensitivities achieveable with modal domain sensor 
systems will be investigated both theoretically and experimentally. And 
finally, a discussion of how modal interference techniques can be used in 
simple, practical sensor systems which avoid the use of large discrete optics 
will be presented. To begin, however, a brief review of several topics relevant 
to the understanding of fiber optics will be presented. 
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2.0 REVIEW OF BASIC FIBER OPTICS 


This chapter primarily consists of a senes of discussions of topics relevant to 
the understanding of basic fiber optics and is meant as a review More 
detailed considerations of all topics may be obtained from their respective 
references. Subjects tc be addressed include linearly polarized (LP) modes 
in weakly guiding fibers, photodetectors and their associated noise factors and 
the photoelastic effect in optical fibers under strain. We will begin, however, 
with a brief review of optical fiber classification. 
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2.1 Optical Fiber Classification 


Fibers used for optical communications and sensing are waveguides made of 
transparent dielectrics whose function is to guide visible and infrared light. A 
typical optical fiber consists of an inner cylinder of glass, called the core, 
surrounded by a cylindrical shell of slightly lower index of refraction called the 
cladding. Optical fibers may be classified in terms of the refractive index 
profile of the core and whether one mode (single-mode) or many modes 
(multimode) are allowed to propagate in the waveguide. If the core has a 
uniform index of refraction ru it is called a step-index fiber. If the core has a 
nonuniform index of refraction that gradually decreases from the center 
toward the core-cladding interface the fiber is called a graded-index fiber. The 
cladding surrounding the core typically is of a uniform index of refraction n». 
where 


n 2 = "i(1 ~ A) (2-1-2) 

in a step-index fiber. This relation is also approximately true for graded-index 
fibers where n » would then represent the index of refraction at the center of the 
core. The parameter A 'is called the core-cladding index difference and 
nominally A £ 001. Figure 1 shows typical dimensions of commonly used 
optical fibers. And finally, as the chapters to follow will be dedicated to 
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describing modal interference as a sensor mechanism in step-index optical 
fibers we note here that a typical value of the core index in a step-index fiber 
is n i = 1.458 [7], We now turn our attention to the guided modes allowed to 
propagate in step-index optical fibers. 


2.2 Linearly Polarized Modes 


Put simply, a mode in an optical fiber waveguide (or any other waveguide) is 
an allowable field configuration, for a given waveguide geometry, that satisfies 
Maxwell's equations tor the derived wave equations) and all the boundary 
conditions of the problem [7]. In an optical fiber a guided mode is one in 
which electromagnetic energy is earned along the fiber axis only and whose 
transverse field components decay to zero in an approximately exponential 
manner at distances far from the core. For these modes the integral 


P = \E\ 2 dxdy, (2.2.1) 

J z=ccnst 

which is a measure of the power flow, is finite. Also, at any given operating 
frequency the number of guided moies is finite with each mode being 
associated with a cutoff frequency That is, for a guided mode to be allowed 

i 
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!o exist, its cutoff must be less than the frequency of operation. It can further 
be shown that a mode will remain guided as long as its propagation constant 
P satisfies the condition 


n 2 < Plk 0 < n, . (2 22) 

where in the case of the step-index fiber, as presently under consideration, iu 
and n? are the tndicies of refraction of the core and cladding, respectively, and 
where ko = 2rr//o, with Xo being the free space source wavelength. Note, 
though, that in order for a guided mode to carry power along the fiber it must 
be excited by an external source. Indeed, since modes are simply analytical 
solutions to the wave equations in a waveguide and represent allowed field 
configurations it is quite possib.e to have mode solutions of a fiber (or other 
waveguide) which do not carry power due to their not being externally excited 
[8], A detailed discussion of selective mode excitation will not be given here, 
but it will be seen later that the excitation of specific mode groups in an optical 
fiber will give rise to some very interesting sensor mechanisms. Recall, 
however, that for light entering a fiber to excite a guided mode it must in 
general satisfy the condition that it enters the core through an acceptance 
cone defined by the numerical •aperture, or NA, of the fiber. A physical 
interpretation of a fiber NA is given in Figure 2 with the analytical value of NA 
given as 
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n 1v 2A 


(2.2.3) 


NA = \ n i ~ n\ = 

in a step-index optical fiber. This result is derived from considering that light 
guided in a fiber must satisfy the total internal reflection requirements at the 
core-cladding interface according to Snell's law of refraction [9], 

In practical optical fibers the refractive indices of the core and cladding are 
nearly equal with n» > n? and A « 1 0. The condition A « 1 0 is called 
the weak guidance condition, and the corresponding optical fibers are called 
weakly guiding fibers. The concept of weak guidance is better appreciated if 
one notes that for A = 0 the infinite clad optical fiber reduces to a 
homogeneous medium with no guidance properties. For A « 1.0 but 
A # 0 the guidance of energy is considered weak. However, this does not 
imply loose confinement of power to the core: at very high frequencies the 
energy can be tightly bound to the core, even in a weakly guiding fiber. 
Furthermore, the condition A « 1.0 is mandated in both single-mode and 
multimode fibers. In multimode fibers A « 1.0 is required for low 
inter-modal dispersion since the delay time between the lowest order and 
highest order modes is proportional to A. In a single mode fiber A « 1.0 is 
necessary in order to have manageable core sizes. In addition, the weak 

A 

guidance condition leads to greatly simplified modal expressions which are 
quite valid for short lengths of fiber, as is the case with typical fiber sensors. 
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In a step-index, lossless and non-magnetic optical fiber that is translationally 
invariant along the axial direction, the guided mode solutions can be 
summarized in phasor form as follows - 


a* 

a y j r<a (2 2.4) 

h(r,o) -■ - y,(e(r, 4>) x a 2 ) 

and 

rcos(r'd)'} ra,') 
elrol = A Kj'.Vr'a) \ 

(.sini'OlJ U,J r> a . (2 2 5) 

h(r.o) = -y 2 (e(r,c>) x a z ) 


rco 

e(r<p) = AJAUr/a K 

( sn 


s(^<p)i r 
> \ 


where there is an implied multiplicative factor of which represents tb 

phase of the mode and where <f = 0,1 2. . , with r, d>. and z being the 
cylindrical coordinates representative of the fiber geometry Also. JAx) is the 
Bessel function of the first kind, order and K c -{x) is the modified Bessel 
function of the second kind, order i J f (x) is plotted for orders 0.1 and 2 in 
Figure 3 The amplitude coefficient A is determined by the source input 
conditions and is assumed to absorb any signs (plus or,minus) at present, 

» while Yi and are the characteristic admittances of the core and cladding, 
respectively. In addition, the parameters U and W are given by 
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V 


jlk 


U = K 0 a x n, 2 -p 2 (2 2 6) 

W = K 0 a^jP - n\ , (2.2.7) 

where “a ’ is the fiber core radius and /) = fl/k a , with /) being the modal 
propagation constant [8] 

Notice in equations (2 2.4) and (2 2.5) that the electric fields are either x or y 
directed Since polarization is defined by the direction of the electric field, the 
modes described by equations (2 2 4) and (2 2.5) are linearly polarized and are 
thus called linearly polarized or LP modes The bracketed quantities in 
equations (2 2.4) and (2 2.5) may be taken as multiplicative factors in any 
combination with the only stipulation being that the magnetic field vectors be 
orthogonally onented with respect to the electric field vectors. Modes 
containing the cos(^) multiplicative factor are called even modes while 
modes containing the sint<?<*>) multiplicative factor are called odd modes. For 
all cases except when £ — 0 the modes have a degeneracy of four (i.e.. even, 
x-polarized; even, y-polarized; odd, x-polarized and odd, y-polarized). When 
£ - 0 it is obvious that LP modes are independent of <t> and therefore have 

9 

a degeneracy of two (i e.. x-polanzed and y-polarized) [8], 
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The propagation constants p for LP modes are found by solving the following 
characteristic equation, the solutions of which are plotted in Figure 4 [9.10] 


MU) K,(W) 

+ WK,_ X (W) 


( 2 . 2 . 8 ) 


Note that the solutions to this characteristic equation actually appear as 
functions of P vs a>. where co is related to the source wavelength /.o by 
oj/.o = 2rtc. with c being the speed of light in free space. The solutions are. 
however, plotted in Figure 4 in the normalized form of b versus V where b is 
the normalized propagation constant given by 

l 'I ik /)) 2 — Ho 

b= f — . (2.29) 

n, — r> 2 


The normalized frequency V is given in terms of the source wavelength and 
fiber parameters by 

V = ft 0 a N ; /?i 2 - nf = k 0 a(NA ) ^ f< 0 an lN 2A . (2.2.10) 

Note also that there are multiple solutions to the characteristic equation for 
each value of (. Therefore, specific LP modes are typically labeled lp { ^, where 
the propagation constant for the mode is the m th solution for p, at a given 
frequency (or source wavelength), when the principal mode number ( takes a 
particular value. For example, the LPn mode is the mode whose propagation 
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constant P is found from the third solution to the characteristic equation when 

e - 2 [ 8 ], 

Guided mode cutoff occurs when b=0 and as can be seen in Figure 4, a fiber 
is operated in a single-mode regime when V <2 405. As V increases (by 
decreasing the source wavelength /.o or by increasing the core radius a) more 
and more modes are allowed to propagate in a fiber until the fiber enters the 
multimode operating regime Strictly speaking a fiber is a multimode fiber if 
more than one mode is allowed to propagate. Practically, though, multimode 
fibers are taken as those fibers with extremely high V values, typically in 
excess of V = 50 To achieve these high V values the core radii are usually 
made very large in multimode fibers, thus making the fiber multimode for a 
wide range of source wavelengths. By comparison, single-mode fibers 
typically have very small core radii. For example, a fiber intended for single 
mode operation at a nominal source wavelength of 1300 nm will have a core 
radius of approximately 4.5 jim (see again Figure 1). Using these values in 
equation (2.2.10), with n i = 1 458, we can solve for A to find A ^ 0 003. It is 
therefore obvious that, as mentioned earlier, A must be much less than unity 

tor singie-moae fibers to nave manageaoie core sizes. 

* 

Of interest in the study of modal domain sensing techniques are those fibers 
which allow the propagation of only a few modes. Such fibers will often be 
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called few-mode fibers, typically with V < 5, and are of interest since it is much 
simpler to study and characterize the interference of a few modes than it is the 
many modes in a true multimode optical fiber. Recall that modal domain 
sensors are essentially self-referenced interferometers which exploit modal 
interference as the basis for a sensor mechanism. Few-mode Fibers are also 
attractive for use in modal domain sensors since, as in all Fiber optic sensors, 
the sensor data will be collected by monitonng the modulated intensity 
distribution at the fiber/sensor output. These intensity distributions (or 
patterns) at the output of modal domain sensors are created due to the 
interference of the guided modes in the fiber and are modulated as a function 
of the environmental disturbances being detected. For modal domain sensors 
using few-mode fibers these output intensity patterns are often qualitatively 
simple and well behaved as a function of disturbance and therefore give rise 
to simple sensor systems. 

One can cause a fiber to be few-moded by simply injecting into a single-mode 
fiber light of a shorter wavelength than that typically injected for single mode 
operation. For example, it will later be shown that a 4 5 /im core radius fiber, 
usually intended for single-mode operation at /o = 1300 nm, will have 
V = 4.616 when the source wavelength is reduced to /o = 633 nm, as with 

a 

He-Ne gas laser injection. We can thus see from Figure 4 that cj'ly the LA», 
LPm, LPit and L Pot modes and their degeneracies will be allowed to 
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propagate. Such a fiber is easily seen to be few-moded. Plots of the intensity 
distributions of the LPoi. LPu . LPn and LP 02 modes can be found in Figures 
5-8. respectively. These and other similar plots to be addressed later were 
created using the numerical expansions lorJo(x) and J i(x) given in Appendix C 
and the FORTRAN programs given in Appendix D. It should be mentioned that 
the x and y values of Figures 5-8 are given in microns and that these figures 
are valid only for r= N x--^y 2 <45 am, as these plots were created 
considering only the modal solutions in the core for the case of 633 nm light 
injection into a 4.5 urn radius fiber, as discussed above. The intensity beyond 
r = 4 5 decays approximately exponentially Finally, it should also be 
mentioned that in Figures 5-8 no interference has been considered and that 
the amplitudes of the modal intensity distributions are quite arbitrary with only 
the general pattern shape being of any importance at present. 


2.3 Photodetectors 


Data is collected from a fiber optic sensor by monitoring the modulated 
intensity distribution at the sensor output. Typically this monitoring is 
achieved through the use of a photodetection circuit. Such a photodetection 
circuit is shown in Figure 9 with its equivalent circuit. It is seen that the 
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photodiode is modeled by a senes bulk resistance R, and a total capacitance 
C tf consisting of junction and packaging capacitances. A bias (or load) 
resistance is given by R and the amplifier following the photodiode has an 
input capacitance C, and a shunt resistance R 4 . In a practical circuit, however, 
R, is much smaller than R and can be neglected. Also, C* is usually small and 
only affects the high frequency response (rise time) of the photodetection 
circuit Thus, for low frequency input signals, such as those created at the 
output of the modal domain sensors to be discussed later, the capacitance C, 
is seen essentially as an open circuit. Furthermore, for the present discussion 
it will be assumed that the amplifier is used as a high impedence buffer with 
unity gain, as would be the case at the input to an oscilloscope, for instance. 
For low frequency or DC signals the photodetection circuit is therefore seen to 
reduce q dite simply to a photodiode in parallel with a load resistance [9]. 

Optical power incident on the photodiode is converted into a photocurrent 
when the excess electron-hole pairs created by the incident light are swept out 
of the photodiode depletion region by the fields created by the high reverse 
bias voltage ( ~ 10 V for pin photodiodes). The photocurrent is related to the 
incident optical power Po by the relation 

a 

V 

I p = &Pq . (2.3.1) 
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where 91 is called the responsivity of the photodiode and has the units of amps 
per watt. is given in turn by the expression 

tiq 

a • {2.32) 

where hv is the energy of the incident photon, q is the electronic charge and 
t] is the photodiode quantum efficiency which relates the number of 
electron-hole pairs created to the number of incident photons. For silicon 
photodiodes a nominal value of responsivity is .#=0.5 A/W for common 
detection wavelengths [9] The photocurrent thus created causes a voltage 
across the load resistor given by 


V 0 = ^PoR l . (2.3.3) 

This is the quantity typically recorded when using simple photodetection 
circuits and the measured quantity Vo is easily converted to a measure of 
incident optical power if the load resistance R L and the responsivity 31 at the 
detection wavelength are known If we now consider that the incident optical 
power Po consists of an average or quadrature point component P Q and a 

signal component AP such that P 0 = P Q + AP we see that the electrical signal 

* 

power c.-eated by the signal component of the incident optical power is given 
by 


2 0 REVIEW OF BASIC FIBER OPTICS 


16 



(2.3.4) 


P £S = (5 A P) 2 R l . 

The pnncipal noises associated with photodetectors are quantum or shot 
noise, dark-current noise generated in the bulk material of the photodiode, and 
the surface leakage current noise. For small and well packaged photodiodes, 
however, the surface leakage current noise is small and can be neglected. The 
shot noise arises from the statistical nature of the production and collection 
of photoelectrons when an optical signal is incident on a photodetector and 
sets a fundamental lower limit on receiver sensitivity when all other 
conditions are optimum The RMS shot noise power is given as 


P SN = 2ql Q BR L . (2.3.5) 

where S is the detection bandwidth and I„ is the average photocurrent given 
as I 3 = 3P 0 where P Q can be measured with no signal applied. 

The photodiode dark current is the current that continues to flow when no light 
is incident on the photodiode. This is in general a combination of bulk and 
surface currents. However, the bulk dark current is usually dominant and 
arises from the electron-hole pairs which are thermally generated in the pn 
junction of the photodiode. The RMS bulk dark current noise is given as 

P DB = 2ql D BR L (2.3.6) 
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where I,, is the primary detector bulk dark current and can be found for 
specific photodiodes in various vendor catalogs [9], 

In addition to the detector noise factors just discussed the photodetector load 
resistor contributes an rtMS thermal (Johnson) noise power given as 

P T = Ak 8 TB , (2.3.7) 

where k 5 is Boltzman's constant and T is the absolute temperature. This noise 
can be reduced by cooling the detection circuit [9]. 

The power signal-to-noise ratio S/iV at the output of an optical receiver is now 
easily seen to be 


SIN = 


^£S 

P SN + P DB + P T 


(2.3 5 


or 


SIN = 


P?A P ) 2 

AkoTB 

2qB(^P 0 + I 0 ) + -|— 


(2.3.9) 


Of special interest is the case when S/N goes to unity. At this point the signal 
power just equals the noise power and enables one to determine the minimum 
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detectable optical signal A P miH . When only photodetector and load resistor 
noise factors are accounted for. we find that 


1 / BkoTB \a 

AP mi „ = 1 4q8(^?P Q + I 0 ) 4 — ) (2.3.10) 

and in the shot noise limited case we find that 

-V’m.n, s,V) = . (2.3.11) 

where a factor of 2 has been included in the noise power terms of equations 
(2 3.10) and (2 3 11) due to the fact that at the operating point the incident 
optical power is uncertain to within ±(shot noise and/or thermal noise). 

Note that this analysis has been for pin photodiodes only and that the shot 

noise limit of detection is practically unrealizable. In addition, if one must 

operate a photodetection circuit at elevated temperatures or if it is impractical 

to use a large load resistor, which would reduce the thermal noise component 

of equation (2.3.1 hi, one may have to resort to the use of an avalanche 

photodiode (or APD). An APD will also give a sharper system rise time for 

detection of high frequency signals but should not be used unless thermal 

noise is excessive This is due to the generation of excess shot noise during 

* 

the avalanche process in APD's. At reasonable temperatures and for ’ 
reasonable values of R the signal-to-noise ratio will be reduced by using an 
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APD. Finally, it should be mentioned that it has been assumed in the previous 
discussions that noise factors arising from sources other than the photodiode 
or its load resistor are negligible. This is reasonable only in a laboratory or 
other situation where environmental factors can be controlled We now turn 
our attention to a discussion of the photoelastic effect in optical fibers under 
strain 


2.4 The Photoelastic Effect 


Recall that the phase of a mode is given by <£ = /Jz - ^ where ^ is a random 
phase term. Since the modal domain sensors to be discussed in the next 
chapter will derive their output from a differential modal phase modulation 
(i e., modal interference) as a function of strain, it is important to discuss how 
the phase of a mode changes as a function of fiber strain. 

As a fiber is strained three principal factors act to change the modal phase. 
First, the fiber length changes. Second, the index of refraction changes due 

• a 

to the photoelastic effect and finally, the fiber diameter reduces. The change 
in phase due to these effects is thus 
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CQ i o 

Aq = — -Az ^ - — An, - —AD 
cz on, cD 


Ad =r plz z— -An, , 


where the effect of the change in fiber diameter is considered negligible and 
where c/J/cn, = /J/n, [1] Also from [1], for low order modes, we have 


An, = | - y n, 3 [c,(t - v)p, 2 - vc^,,] | . (2.42) 


where c. is the axial strain x is Poisson s ratio. n> is the index of the core and 
p»i and p<t are the photoelastic constants. Reasonable values of these unitless 
constants are given for silica as follows [4] - 


= 1.458 


By substituting equation (2 4 2) into equation (2.4 1) and rearranging terms we 
find that 


A(j> = /?Az|l -yn, 2 [(1 - v)p 12 - xp n ]| 
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Aq = PAz{z). 


(2.4 3) 


where the relation z x e r = Az has been used. Using the constants given 
above we see that a = 0.78. Although this is a reasonable value of a, the true 
value may change from fiber to fiber Experimentally determined values have 
been found as large as j = 0.92, as will be seen later. In any case, though, 
it is clear that due to the photoelastic effect a fiber under strain will see a 
significant reduction of the core index of refraction and therefore a significant 
reduction in the effective optical pathlength of the modes in the fiber. Th.. 
results in a change in modal phase in a fiber under strain which is less than 
that expected if the fiber simply undergoes pure elongation 

When a fiber supports more than one mode it is also important to investigate 
how the differential phase of the modes is influenced by strain. Consider two 
modes with phases <i>, = fl,z — i//, and <p* = /Lr — i > 2 , respectively. The 
differential phase is given as 

$12 = ~ Pt) z ~ ($1 — $ 2 ) = ~ • (2.4.4) 

We may now consider that <*> 1 * = 4>o + A$ where <po is the initial differential 
phase of the modes in a no strain situation and where A<t> 1 * is the differential 
phase change induced as the^fiber is strained. In an exactly analogous 
manner that was used to prove the result of equation (2.4.3) we can show that 
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A<p )2 = A0 Az {x} . 


(2 4 5) 


so that 


0,2 = A pcM*) + 4>o - (2 4 6) 

where e z L 1 =Az and L, is the gage length of, say. a fiber optic sensor. Note 
also that <po = -A0. As will be evident in the next chapter the relationship 
given in equation (2 4 6) is extremely important in the study of fibei optic strain 
sensors. 

One last comment should be made before concluding this discussion That is. 
note that in the previous analyses it has been assumed that the amplitudes of 
the modes in a fiber under strain remain constant and that only the modal 
phase is influenced as strain is varied. This is not exactly the case since the 
amplitude of a mode is a function of U = A 0 a v ' n, 2 — /) 2 , in the core region, and 
both n, and p are functions of strain. However, only very large amounts of 
strain vary n » or /? appreciably so that U will essentially remain constant with 
strain. For instance. 

/7, = n,(^) as 1.458 + An(eJ = 1.458 - (0.32)^ (2.4.7) 

9 

using the previous calculation for An in equation (2.5 2) It is thus seen that for 
strains even on the order of one percent (which is near the breaking point of 
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most fibers) n» remains nearly constant. A similar result can be obtained for 
]} as a function of strain, making the assumption of constant field amplitudes 
quite valid. 
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Figure 1. Optical fiber geometry summary 
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Figure 2. Numerical aperture (NA) of an optical fiber [16] 
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Figure 6. Intensity distribution ot the LP jt modo 
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Figure 8. Intensity distribution of the LP ct mode 


2.0 REVIEW CF BASIC FIBER OPTICS 


32 


4 
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Figure 9. A simple photodetector and its equivalent circuit [9] 
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3.0 MODAL DOMAIN SENSING 


Modal domain sensing is a fiber optic sensing technique based on the 
differential phase modulation occuring as a function of disturbance between 
the interfering guided modes of an optical fiber As we know, specific mode 
groups will propagate in a fiber and interfere to cause characteristic intensity 
distributions, or patterns, at the fiber output, as is shown in a very simple 
fashion in Figure 10. When the fiber is subjected to a disturbance, these 
intensity patterns are modulated as a direct function of the external 
disturbance and by monitoring this pattern modulation, often by using simple 
spatial filters, one can obtain information about the disturbance. 


Strictly speaking, any sensor based on the interference of modes in an oDtical. 
fiber can be called a '‘modal domain sensor." However, of special interest are 
those sensors based on modal interference in few-mode optical fibers with 
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V<5 0. Such sensors are qualitatively simple in that since only a few modes 
are allowed to propagate it is possible to obtain simple, closed form solutions 
for the modulated output patterns which may manifest themselves. This is 
obviously an advantage when one wishes to design a predictable sensor. As 
a result, only modal domain sensors utilizing few-mode fibers will be 
discussed in the remainder of this chapter 

Specifically, this chapter will be dedicated to descnbing the electromagnetic 

basis of the observed intensity pattern modulation, as a function of axial strain, 

of modal domain sensors using 9um core diameter fiber at source 

wavelengths of 633 nm and 850 nm For simplicity we will use the notation 

MDS/9-633 and MDS/9-850 to describe these modal domain sensors at the 

source wavelengths of 633 nm and 850 nm, respectively. Also, note that 

although only axial strain considerations will be described here, modal 

domain sensors are not limited to measuring only axial strain. Axial strain is 

simply one of the most fundamental disturbances to which a fiber can be 

subjected and a detailed analysis of how axial strain affects an optical fiber 

sensor will provide the basis for designing sensors to measure other 

disturbances. For instance, if a sensing length of fiber is embedded in a 

material with known thermal expansion characteristics, a thermal sensor 

based on modal interference could be designed since the thermal expansion 

* 

of the material can be related to induced axial strain in the embedded fiber 
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sensor Sensors designed to detect many other environmental influences are 
also possible We now proceed to a detailed discussion of the MDS/9-633 
sensor. 


3.1 The MDS/9-633 Sensor 


It has been observed that when 633 nm light from a He-Ne gas laser is injected 
into a 9 pm core diameter fiber a symmetric three lobe pattern can be created 
at the output which exhibits pure rotation when the fiber is subjected to axial 
strain Using a Photon Kinetics FQA-2000 fiber optic analyzer, the second 
mode cutoff for the above mentioned fiber, which is designed for single-mode 
operation at 1300 nm, has been determined to be 1215 nm (i.e., V = 2.405 at 
/.o = 1215 nm). Also, this fiber has a measured numerical aperture of NA = 0.1 
(refer to Figure 1 1 for plots of the data generated by the FOA-2000 unit) At a 
source wavelength of 633 nm this fiber is easily seen to have a normalized 
frequency of 


V = 2.405 = 4.616 . 

€33 nm 

This allows the propagation of the LPoi, LPn, LPu and LPoi modes and their 

9 

degeneracies, as discussed in section 22. However, as the three lobe pattern 
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under discussion exhibits a central null, it is evident from Figures 5 and 8 that 
the LPoi and LPoj must not be excited in the fiber for the three lob «2 pattern to 
be created. If we further require that our source emits polarized light, say 
x-polarized, we are left with only four modes which may interfere to cause the 
rotating three lobe pattern of interest. Specifically, these modes are the 
x-polarized odd LP »». even LPn. odd LP2, and even LPi « modes and it can be 
shown that with a proper choice of amplitude coefficients, determined in reality 
by the launch conditions of the fiber, that these modes do interfere to cause 
the observed three lobe pattern and its well defined rotational behavior [11]. 
From equation (2 2 4 ) we see that the electric fields in the core of the above 
mentioned modes are given as 

ff, = Af J 1 (U u rla)cos<p e~ j(P " z ~ v,,) a x 

E° } = A° J,(Uy } r/a)smQ e~ liP " z ~ v " ] a x 

x ( 3 -I 1) 

£21 = A%J 2 [U 2 irla) c os26e~ Jil!zZ ~ U2 ' ) a x 
ff, = A2J 2 (.U 2 irla)sm2$e~ j{t,yZ ~ v *' ) a x 

Since the sensor behavior can be completely described by investigating only 
the interference of the mode solutions in the core, no further consideration will 
be given to the modal solutions for r > a. 

* 

The intensity pattern present at the fiber output is found from 
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I = Re[(E x a x x Y y E x a y ) • a 2 ] = ± Y,\E X \ 2 . 


( 3 . 12 ) 


where 


E x = (E^ - £?,) - [Ef, - £ 2 c ,) (3.1.3) 

and Yt is the characteristic admittance of the core region 

We now assume that the amplitude coefficients are real and that 
Af = Af = A , and At = A? = A 2 with the further assumption that the odd 
terms of equation *3 1 1i are 9C 3 out of phase with relation to their respective 
even terms, so that 


*11 = *11 + 90° = *11 ~ 90° 

and 

*2t = *21 - 90 ° = ^21 ~ 90 ° 

Adding the even and odd ternis of the LPu mode, it is easy to see that 

{Eft Eft) = A, J\ ( j ( cos <t> + jsm&e-M'' 1 -*"' 

* 

or 


f 

* 
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(£f, + £f,) = /»i*/i(-^r)e 


(3 1 4) 


-]y$uZ - V,| - ^ 


Similarly, adding the even and odd terms of the LPzi mode it is seen that 

(£21 + £|,) = * 2 J 2 (~T-) e~ fl **'~ Vi ' + 2 * } (3 1 5) 


With these last two results in mind we now find the output intensity distribution 
using equations (3.1 2) and (3 1 3). The result is 



cos (Apz 


-Ail/ — 3 <i>) 1 , 


(3.1.6) 


where Afl = /?„ - /?,, , Ail/ = i>,, - 1 J/ 2i and I 0 = Y.J2. 


It should be evident that the mtens'ty distribution of equation (3.1.6) results in 
a pattern consisting of three equal lobes displaced azimuthally by 120°. It 
should also be evident that this pattern rotates as the modes propagate 
through the fiber. The first two terms of equation (3.1 6 ) are simply 
background intensity terms and are constant. The sum of these two terms is 
plotted in Figure 12, where it is assumed that A< = Aa = 1.0 and where the x 
and y values are 7 'ven in microns The third term of equation (3.1.6) is the 
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term which gives rise to the pattern shape and rotational behavior as the 
modes advance in the z-direction and is plotted in Figure 13. The entire 
intensity distribution of equation {3 1.6) is plotted in Figure 14. As with Figure 
12. Figures 13 and 14 are plotted assuming Ai = Aj = 1.0 , with the x and y 
values being given in microns. The choice of At and Aa was made after 
evaluating several plots of equation (3 1 6). with the amplitude coefficients 
taken as parameters For At Aa it was seen that there were approximate 
radial intensity nulls in the calculated intensity distrubution, displaced by 120°, 
which agree very closely with experimental observations. In addition, the 
plots in Figures 12-14 are valid only for the core region of the fiber (i.e., the 
plots are valid only for r= N> x--^y ; s 4 5 ^m) and were plotted using the 
Bessel function mumerical expansions given in Appendix C and the FORTRAN 
programs of Appendix D. Notice also that for At = A*, the amplitudes of i. 
plots of Figures 13 and 14 are approximately equal in the radial region where 
the plotted functions are maximum, indicating that for this region in equation 
(3.1.6) 


+ . (3 ,. 7) 


Using this result we can simplify equation (3.1.6) at the radial position where 
the function goes through a maximum to find 
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I = r 0 ri + cos(A/?z - Av - 3<p)] 



cos 


21 Apz - Ai If ~ 3o 


(3.1.8) 


where 1' 0 has incorporated into it the terms of equation (3 1.7) after they are 
factored from equation (3.1 6) In addition to obtain the result of equation 
(3 1 8), the trigonometric identity 1 - 1 - cos(2x) = 2 cos : (x) has been used, with 
the peak lobe intensity l 3 taken as I P = 21' 0 Now let us discuss how this 
rotating three lobe pattern can be used to detect axial strain. 


If we observe the three lobe pattern as it emerges from the end face of the 
fiber, we will see an intensity distribution similar to that shown in Figure 15(a). 
If we then increase the optical pathlength of the fiber by inducing an axial 
strain, we will observe a pattern rotation. Due to the cosine squared nature 
of this pattern along the circular path which contains the lobe peaks, as given 
in equation (3 1.8) and as shown in Figures 15(a) and 15(b), it is seen that the 
most sensitive and linear position for monitoring the pattern rotation for the 
detection of axial strain is along this path at a position where the intensity is 
half that of the intensity of a lobe peak, as shown in Figure 15(b). In other 
words, if in a no strain situation we choose a fixed quadrature point in .space 
to monitor on the three lobe pattern that both has an intensity of exactly half 
that of the peak intensity and is located on the circular path that contains the 
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lobe peaks, then for small increases in strain the intensity at this point will 
increase or decrease, as the patterns rotates, in direct proportion to the 
increase in strain. 


If we define the azimuthal coordinate o to take on the value <p = 0 at our 
quadrature point, then equation (3 1 8) becomes 


I = 


I p cos 


i 


Aflz — AO 


) 


(3 19) 


By now noticing that tre A/?z — term in equation (3 1 9) represents an 
ensemble differential phase of the mooes contributing to the rotating three 
lobe pattern, we can use arguments analogous to those used to prove the 
result of equation (2.4 6) so that the photoelastic effects may be included in 
equation (3.1 9). as the fiber is strained, to give 


I = I p cos 


2 /' ^ L f 2 * \ 


(3 1 10) 


where c, is the axial strain, L, is the gage length of the fiber sensor, a is a 

constant less than unity that arises due to photoelastic effects and <t> 0 is the 

* 

ensemble differential phase of the modes in a no strain situation. If we now 
choose the quadrature point, or Q-point, such that the intensity at this point 
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increases with strain from an unstrained value of l = L/2. we must require 
that (f> 0 = - r.,'2 so that at the O-point 



21 &P h L r x 
C0S (— 2 — 



(3 1.11) 


In practice however, once the quadrature point has been set it will often be 
of more use to measure the change in intensity from the O-point as a function 
of strain. This is found simply to be 

Al(Ar r ) = = I c (SpL f a)Ac z . (3 1 12) 

where I 0 = I s /2. Note that this result was obtained from an assumed no strain 
initial condition: however, it is still valid if the quadrature point is set with the 
fiber initially under strain. This is useful in a practical situation so that ± Ae t 
can be detected. Also note that this result is valid only for small changes in 
axial strain and that for small 1 Ac, i . the function A^Ae,) is linear 

Some other observations should be made at this point. First of all, since the 
response AI^eJ is linear, we can use the same sensor arrangement to 
measure static strain or small amplitude sinusoidal strain, where 

AI(Ac pea „.f) = l Q {SP L, x) Ac pe3k cos(uf) . (3.1 13) 
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if A£. m , is the amplitude of the sinusoidal axial strain. This result v.'ll also hold 
for small amplitude vibrational strain such as would be caused if a fiber was 
bonded to a vibrating structure. “Small amplitude vibrational strain ' is meant 
to indicate that the peak deflection of the vibrating structure is such that the 
fiber is tensioned without being appreciably bent. Second, since for small 
amplitude sinusoidal perturbations in the fiber the resulting change in intensity 
is also sinusoidal at the same frequency, the sensor has obvious uses as a 
vibration sensor. And finally, it should be noted that in an actual sensor 
arrangement, a detection circuit will be used to convert intensity to voltage. 
If a small spatially filtering pinhole is placed between the fiber end face and 
the photcdetector to monitor the intensity pattern at the Q-point this voltage 
will be directly related to the area of the spatial filter, the responsivity of the 
photodetector and the load resistance of the detection circuit such that 

AV(Ac z ) = V Q (Ap L f 7) M z 


or 


^ Vq (A/J Lf x) 


(3 1 . 14 ) 


for the case of static strain in the fiber and 

SV(Se pe3k .t) = V Q (ApL f 2)AI peak cos (tot) 
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or 


±V peak cos(wf) 
= V 0 WL,,) 


(3 1.15) 


for the case of small amplitude ti.ne harmonic strain and small peak 
displacement vibrational strain. Note that in equations (3.1.14) and (3.1.15) the 
notation Si, has been used to indicate the average measured axial strain in the 
fiber This has been done since strain is a localized quantity and fiber opti r 
strain sensors are only able to detect strain as averaged along their gage 
length. See Appendix A for further details of strain averaging In addition, the 
results of equations (3.1.14) and |3.1 15) will likely be of most use in an actual 
system with V 0 measured at quadrature and and o measured, say, using 
a spectrum analyzer. 

Before discussing sensitivity let us discuss some of the practical aspects of the 
MDS/9-633 sensor. First, it should be noted that the interest in studying the 
rotating three lobe pattern of the MDS/9-633 sensor arose due to its ability to 
be created using simple optical components. Also, once created, this pattern 
is easy to work with due to its well defined rotational behavior and due to the 
fact that 633 nm light is visible and requires no special viewing equipment. 
Creating the rotating three lobe pattern, however, is not at all an easy task. 
Indeed, it often seems that obtaining this pattern requires a fair amount of luck 
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Once created, though the modal content of the pattern is unique and can be 
studied, as has been done in the past few pages, leading to the simple strain 
measurement results of equations (3 1.14) and (3.1.15). These equations state 
very simply how small changes in strain can be measured in the optical fiber 
of a MDS/9-633 sensor once the quadrature point has been properly set. 

Setting the O-point of the MDS/9-633 sensor, however, is another difficult task 
to accomplish Typically, in order to position a small spatial filter at the proper 
quadrature point, one must apply a small signal sinusoidal disturbance to the 
sensor via a piezoelectric device and use linear positioners to move the 
spatial filter until one locates the position of the pattern that gives the largest 
amplitude and least distorted sinusoidal signal at the output of the 
photodetection circuit, as displayed on. say, an oscilloscope. Note that the 
spatial filter must be placed between the output end of the fiber sensor and the 
photodetector and must be small so that any intensity variation, which is a 
point function, can be accurately simulated and monitored. One may also set 
the quadrature point by aligning the spatial filter at the peak of a lobe and by 
then causing a static strain in the fiber, via some static strain control device, 
so that the pattern rotates until the quadrature point is property aligned with 
the spatral filter The use of both of these techniques for setting the Q-point 
will be discussed in chapter 4. 
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It should be obvious now that the MDS/9-633 sensor has practical limitations. 
In addition, due to the expected ease of use of the MDS/9-850 sensor to be 
discussed in the next section, no analysis has been given for the injection 
conditions necessary for the rotating three lobe pattern to exist in a MDS/9-633 
sensor Therefore, from this point on it will be assumed that if one wishes to 
use a MDS/9-633 sensor that a rotating three lobe pattern can be obtained and 
that there is available equipment which allows the quadrature point to be set. 
The results of equations (3.1 14) and (3.1.15) can then be used With this in 
mind the sensitivity of the MDS/9-633 sensor can now be addressed 


Although in a typical measurement situation one would measure Ac, as a 
function of the change in voltage at the output of a photodetection circuit with 
respect to its Q-point value, we may still relate Ac, to the change in optical 
power through the spatial filter with respect to its Q-point value This is done 
simply by dividing the numerator and denominator of equation t3 1.14), say. 
by &R l , where 5? is the responsivity of the photodiode and R L is the load 
resistance of the detection circuit This yields 


Ac, = 


A P __ 
pQ&l}L f cx 


( 3 . 1 . 16 ) 


By refenng to Figure 4 we can see that b,, = 0 55 and b,, = 0 22 . Using 
equations (2 2.9) and (2.2.10) and Figure 4 we can therefore find that A/J for the 
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MDS/9-633 se.’.sor is A/j' = 1 2x10 '• m ' when V = 4 616, n, = 1 458. a = 4 5 pm 
and /„ = 633 nm In addition we will reasonably assume for our detection 
circuit that for pm diode detection the fiber gage length L, is 0 5 m , the bulk 
dark current I, of the photodiode is approximately 1 nA, the detection 
bandwidth B is 1 Hz. the responsivity Ji of the photodiode is approximately 0.5 
A/W at 633 nm, the operating temperature T is 300°K and the load resistance 
is R = 10 kQ We further assume that the optical power through the spatial 
filter at quadrature is P ; = 1 uW Using these values we may calculate the 
minimum detectable change in optical power through the spatial filter using 
equation (2 3.10) We thus find 

AP min = 3 8x10 - ' 2 VJ 

If we now substitute this value into equation (3 1.14) and assume a = 0.78 as 
discussed in section 2 4 we find that a reasonable minimum detectable strain 
expected for a typical MDS/9-633 sensor is 

AZ min s 10 -3 m/m 

The MDS/9-633 sensor is thus seen to be a highly sensitive strain detector It 
is not, however, as sensitive as say a Mach-Zehnder interferometer since 
modal domain tensors are based on differential modal phase modulation, 
whereas the interferometer is based on absolute phase modulation [2). By 
sacrificing some sensitivity, though the modal domain sensor will likely gain 
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increased stability and immunity to environmental noise over the 
interferometer, making modal domain sensors, in general, a more practical 
approach to strain detection in many real world situations. 


3.2 The MDSI 9-850 Sensor 


Similar to the MDS/9-633 sensor, it has been observed that when 850 nm light 
from a semiconductor laser diode is injected into a 9 ptm core diameter fiber 
a symmetric two lobe pattern can be created at the fiber output which 
oscillates as a function of axial strain Considering the same fiber as discussed 
in section (3.1) we see that the normalized frequency at /, = 850 nm is 


V = 2 405 


1215 nm 
850 nm 


3 44 . 


From Figure 4 we find that this allows for the propagation of only the lp„ and 
u>,, modes, including all degeneracies We will now show for the case of 
x-polarized injection that the oscillatory two lobe pattern is caused by the 
interference of the lp 0 , and even lp„ modes 


From equation (2.2.4) we see that the electric f>plds in the fiber core of the 
x-polarized lp 0 , and even lp„ are. respectively 
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(3 2 1) 


£ 01 = A,J 0 (U Qi na)e- J,p ^- v --\ 


and 


= 4 2 J,(U n r/a) cos <> e J<p " z v,, 'a x 


(3 2 2) 


Using equation (3.1.2) we can find the intensity distribution created by the 
interference of these modes to be 



In equation (3 2 3) A/? = /?„, — /?„ , Ai 1/ = \p 0 , — ii/ u and 1, = Y,/2 . In addition, 
from Figure 4 we can find that b 0l = 0.730 and b H = 0 325, so that if equations 
(2 2.9) and (2 2 10) and Figure 4 are used with V = 3.44 , n, = 1.458, 
a = 4 5 /im and /. e = 850 nm we calculate A/f to be A/3 = 1 .1x10’ m~'. 


Note that in equation (3 2.2) the first term is the intensity of the lp 0 , mode alone 
and the second term is the intensity of the ir n mode alone while the third term 
arises due to the interference of the »p 01 and u>„ modes. For discussion 

*e 

purposes let us rewrite equation (3.2.3) as 
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I = UM r ) -r M r . <p) -r i 3 lr.<p. Z)\ 


(3 2 4) 


Notice that I,(r) is azimuthally symmetric while I 2 (r, <t > ) has a maximum at 
4> - 0 and 4> = z. The third term I 3 (r, 4>, z) also has a maximum at <t> = 0 but 
has a minimum at o — rt In addition I 3 (r, <f>. z) is modified by the 
cos {Apz — Ac/) term which varies between -1 and +1 as z is varied. Thus we 
can see that for proper choices of A, and A, we should expect to see an output 
intensity pattern consisting of two lobes, one of which gets brighter as z varies, 
while the other gets dimmer. Figures 16-22 help visualize this behavior. 
Figures 16 and 17 show the intensity distributions of I,(r) and l 2 (r, <f>), 
respectively, while Figure 18 shows a sketch of the u»,, mode, which is similar 
in the coordinate & to the output of the MDS/9-850 sensor at quadrature. 
Figure 19 shows I,(r) l-(r, <p) when I 3 (r, <t>. z) = 0 , while Figure 20 shows 

I 3 (r. <?, z) at a maximum. Finally. Figure 21 shows the entire intensity 
distribution of equation (3 2.4) when l 3 (r, <t>, z) is a maximum, with Figure 22 
being a three dimensional plot of I,(r) + I,(r, <p) when l-(r. <p,z ) = 0. In these 
figures it is assumed that A, = 2A,, with all but Figure 18 showing plots 
normalized to I,(r) having a maximum intensity of unity. The choices of A ( and 
A 2 are somewhat arbitrary and give indications as to the ratio of the input 
power existing in each mode. Here the choice of = 2A, was done to 
- enhance the visualization of the behavior of the two lobes as a function of 
strain, which has been observed. This choice was made after evaluating 
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several plots of equation (3 2 3), with the amplitude coefficients taken as 
parameters With A 2 s 2 A, it was seen that the calculated output intensity 
distribution of the MDS/9-850 sensor most closely agrees with experimental 
observation One must realize, however, that the ratio AJA, will be 
determined by source input conditions and must be some moderate value 
greater than one If AJA, is less than unity the lp., mode will dominate and 
only a single lobe will present at the output. If A.,' A. is much greater than unity 
only negligible power will exist in the u> gi mode and thus no interference will 
take place For convenience we will proceed under the assumption that 
A, = 2 A, 


Some further comments concerning Figures 16-22 should be made. First, it 
should be noted that all line graphs are plotted in the plane where I ,(r) , 
I : (r, <f>) and I 3 (r. <t>, z) have maxima. Second, all the units of the spa* - 
dimensions in Figures 16-22 are in microns. Third, the 3-D plot in Figure 22 
is only valid in the region where r=^x* + y : <4.5 pm And finally, the 
sketches of Figure 18 are included to give a better visualization as to how the 
pattern caused by lp qi and lp,, modal interference will appear as seen looking 
toward the fiber end face As with similar plots in earlier sections, these plots, 
with the exception of Figure 18, were generated with aid of the material in both 
Append.* C and Appendix D. 
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At this point it should be easily seen how this modal interference phenomena 
can be used as a sensor mechanism If one adjusts the input conditions in 
such a way that I 3 (r, 4>, z) is zero at the fiber output end face, a maximum 
sensitivity quadrature point will have been c r eated If one then monitors one 
of the peaks at this quadrature point by spatially filtering out all but a small 
portion of the output pattern at this peak, a change in optical power through 
the aperture of the spatial filter will be directly related to any perturbations in 
the fiber 


Notice now that in Figure 21 if we monitor the peak of a lobe, then a maximum 
intensity at some amount of strain will occur at r s 2.38 um in the maximum 
intensity plane If we choose the lobe at say <p - 0 and if it can be shown that 



+ 4 J 





(3.2.5) 


at r s 2.38 pm, for A 2 ^ 2A,, then we can obtain the reduced form for the core 
field intensity, equation (3 2 3), as 


I = I' 0 {1 + cos (&P z — Ai/0} 


or 

I = Ip cos ^ 2 J 
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(3.2.6) 



where F 0 has incorporated into it the terms of equation (3 2 5) after being 
factored from equation (3 2.3) and where l p = 2l' 5 represents the peak intensity 
of a lobe. Indeed, by recalling the fiber parameters and using Figure 3 we can 
show for 2.38 um that the quantities of equation (3 2 5) are the same to 
within 4 8 percent, with A 2 = 2A ,, so that the above simplification of the core 
intensity function is valid in the region of a lobe peak at either <p = 0 or 
o — n The only difference, recall, in these two lobes is that one will increase 
in intensity while the other decreases for small amounts of strain. Also, since 
the peaks of I 2 (r, o). I 3 (r. <$>, z) and I.(r) + 1 2 (r, <?. z) don't precisely align with 
one another, with the peak radius of Mr. o) being the largest, it is expected 
that the above simplification of the core intensity function should be valid for 
reasonable ratios of A 2 jA , > 2 m the regiorrs~of the lobe peaks at <? = 0 and 
$ = n. Therefore from now on it will be assumed that we are discussing a 
sensor system in which the input has been adjusted so that two lobes of equal 
intensity are observed at the output under a no strain situation Further, it will 
be assumed that the ratio A : jA, is some reasonable value so that the 
interference takes place and that in the region about the peaks of the lobes the 
intensity can be approximated by equation (3.2.6). 

We recall that by using equation (2.4.6), we can include in equation (3.2.6) the 

photoelastic effects of the fiber under strain, at the radial coordinate of a lobe 

* 

peak, to give 
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(3 2.7) 


I = Ip cos 


Ap t z L f 



where <£„ = - n/2 if we consider that we are monitoring the lobe which 
increases with strain from I = L,2 under a no strain situation. We recall that 
at quadrature I 3 (r, q z) = 0 and we have two lobes of equal intensity, one of 
which increases while the other decreases as a function fo strain. Using the 
same arguments as were used to arrive at the results of equations (3 1 14) and 
(3.1.15), we can show for small changes of strain, as averaged over the gage 
length of the sensing fiber, that AZ, can be measured from quadrature with the 
MDS/9-850 sensor as 


Sc f 


AV 

V Q (Ap L,*) 


for the case of static strain and 


(3 2 8) 


AVpe 3 * c °sM) 

Vq (AP L, J) 


(3.2.9) 


for the case of small amplitude time harmonic strain and small peak 
displacement vibrational strain. Note the similarity of these results (only AP 
takes on a different value) and those of equations (3.1.14) a.^d (3.1.15) even 
though the sensing techniques are rather different The sensitivities of the 
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MOS/9-850 and MDS/9-633 sensors are also comparable and equations (3.2 8) 
and (3.2 9) may be used to detect ~ At, from quadrature if there is an initial 
strain in the fiber 

As expected, using a laser diode in a modal domain sensor provides several 
advantages. First, due to the fundamental mode (LPot) not needing to be 
extinguished for proper sensor operation, it is much easier to align the input. 
Also, since information will be taken from the sensor by monitoring the 
intensity variations of a lobe peak, alignment of the Q-point of the sensor is 
expected to be relatively simple. In addition, laser diodes are light weight, 
physically small and it is possible to find highly linearly polarized laser diodes 
operating at /„ = 780 nm that are relatively inexpensive 

Although the preceding analysis has been for /. = 850 nm. it is still valid for 
any source wavelength such that V < 3.8 so that the lp 2 , mode is still cutoff. 
For X„ = 780 nm, \ / = 3 75 In addition, for /. 0 = 780 nm, A/5 s 1 1x10'* m~' as with 
850 nm injection so that sensitivity is preserved. 

The only apparent drawback to using laser diodes as sources is that they 
require special viewing equipment to observe the infrared radiation. This 
problem may, however, be overcome through the use of small, light weight 
in-line detection techniques presently under investigation. And finally, as 
laser diodes may be operated in a pulsed fashion, it may be possible to 


3 0 MODAL DOMAIN SENSING 


55 



simulate the use of a beam chopper in a detection system so that 
signal-to-noise ratio may be preserved for low frequency noise as would be 
present if there was, say, excess thermal drift. 
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Figure 10 A simplistic look at model interference 
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Figure 16. Th* LP w mod* intensity distribution of the MDS/9-8S0 sensor 
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Figure 17. The L P J? mod* intensity distribution of the MDS/9-850 sensor 
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Figure 19. MDS/3-850 Intensity variation of l,(r) + l^r.^). at quadrature 
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Figure 20. Intensity varictlon of the MDS/9-850 cross term I,(r,^.z). at a maximum 
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Flgur* 21. Total intensity pattern ol the MDS/9-850 output when ljfr.^^z) is maximum 
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4.0 MDS/9-633 SENSOR EVALUATION 


Although an MDS/9-633 sensor system will have practical limitations, it is easy 
to construct a worhng sensor system based of the MDS/9-633 technique using 
readily available optical laboratory equipment In addition, by evaluating 
MDS/9-633 sensor one can obtain insight about the practical aspects ana 
performance characteristics of modal domain sensors in general. We begin 
with a description of the experimental system used to evaluate the MDS/9-633 
sensor technique. 

The experimental setup is shown diagrammatically in Figure 23 As seen. tr.. 
emission from a He-Ne gas laser is focused through a microscope objective 
onto the end face of a 9^m core diameter fiber held in a fiber positioner. This 
fiber is of the same characteristics as the fiber discussed in section 3.1 and is 
attached in sequence to a static strain control deVice, a piezoelectric PZT 
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cylinder, a cantilever beam and a holding post at the output. The static strain 
control device simply consists of two plastic disks on which the fiber is wound, 
one of which has an attached push rod and is mounted on a rotateable base. 
As seen in Figure 24. tension is applied to the fiber between the two disks as 
force is placed on the push rod via a micrometer screw. In addition, so as to 
avoid appreciable bend loss in the fiber, the radius of the disks was 
approximately three inches 

The PZT cylinder following the static strain control device had wound upon it 
several turns of the 9 ^m core diameter fiber and obtained its excitation from 
an external signal generator As with the static strain control device, the 
purpose of the PZT cylinder was to aid in aligning the syctem quadrature point. 
For reasons to be discussed later, however, only the PZT cylinder was able to 
help achieve this goal. 

After the PZT cylinder, the fiber was attached, along with a resistive strain 
gage, to a cantilever beam. Appendix A presents a detailed analysis of the 
strain characteristics of a loaded cantilever beam, as well as a discussion of 
the geometry of the fiber path on the beam necessary so that the strain 
induced in the optical fiber is identical to that induced in a resistive strain gage 

S 

placed at the geometric center of the beam surface. The reader is encouraged 
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to review Appendix A before proceeding. Referring to Appendix A, the 
dimensions of the cantilever beam and the optical fiber path were chosen as 


L « 0 288 m 
L, = 0 0254 m 
b - 0 071 m 
h = 0 002 2 m 
r = 0 0159 m 


where L is the beam length, L-U is the length of the straight fiber segment^ 
and r is the radius of the curved portion of the fiber, with b and h being the 
beam width and thickness, respectively With this geometry the change in 
average strain induced m both the optical fiber and the resistive strain gage 
is given as a function of load by 


SI = 


3mgL 

Ebh 2 


(4.0 1) 


as shown in Appendix A, where m is the change in initial load mass and g is 
the acceleration due to gravity. Young's modulus E for the beam depends on 
the beam material and as our beam was constructed of alumin. 

E = 6 9x10 10 Pa. In addition, for this beam the Poisson's ratio is v=0.33, while 
the gage length of the Tiber on the beam is 

L f — 2 [L — Lj) -f- nr = 0.575 m . 
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Further, the optical fiber was mounted to the beam unstripped to prevent 
breaking and both the fiber and resistive strain gage were bonded to the beam 
using thermal-set epoxy 

The output end of the fiber sensor was attached to a holding post positioned 
a few centimeters from the photodetection circuit The detection circuit 
consisted simply of a photodiode, reverse biased at nine volts, mounted on an 
x-y-z positioning stage Load resistance was provided via a decade box In 
addition, epoxied to the front of the photodiode (UDT model PIN-10D) was an 
iris which served as a spatial filtering device. Fully open the iris was 1 2 cm 
in diameter and when closed a 1.0 mm aperture remained. Typically, spatial 
filtering was done with the iris closed; however, the actual size of the spatial 
filter aperture is mostly irrelevant as long as it is small with, respect to th. 
output intensity distribution being monitored. 

Not shown is the Wheatstone bridge and amplifier circuit used to detect small 
changes in resistance of the resistive strain gage. The output of this circuit 
and the photodetector were both connected in parallel to a digital voltmeter 
an oscilloscope and a spectrum analyzer for data collection purposes, 
addition, the entire system, except for measurement equipment, was mounted, 
to a floating (pneumatically supported) optical table. 
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Anally, throughout the discussions to follow no account will be given to any 
system noise arising from sources other than the photodetection circuit and its 
associated load resistor In addition, it will be assumed that the bonding of 
both the optical fiber and the resistive strain gage to the cantilever beam is 
perfect with no shearing of the fiber tacket occunng as strain is applied These 
assumptions will not likely be valid outside of a laboratory environment 


4.1 Resistive Strain Gage Evaluation 


As the resistive strain gage is a commonly used and well accepted instrument 
for measuring strain it was decided that a resistive strain gage would be used 
as a standard with which to compare the performance of the MDS/9-633 
sensor. It was therefore important to evaluate the performance of the resistive 
strain gage used in our sensor system prior to all other measurements. As 
mentioned previously, the resistive strain gage was placed at the geomef’"" 
center of the cantilever beam surface, so that as a function of load, the inc. 
strain characteristics of the gage follow equation (4 0.1). In addition, as shown 
in Appendix B, the change in average strain in the beam as measured by the 
resistive strain gage is given by 
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t - fyo 


(4 1 1) 


-Vqu* 4 *S 
V, n S s R s 


1 — k t v 


where V :ut is the negative voltage measured (for positive strain) at the output 
of the inverting amplifier following the Wheatstone bridge circuit. The bridge, 
recall is used to allow small changes in strain to be detected. In addition the 
input V,„ to the bridge circuit, as supplied by a regulated power supply, w?- 
measured 10 be 10 0 V, while the gain of the amplifier was measured to be 
| R s //7 3 | = 179.78. Individual values of R» and Rt were measured to be 
102 834 and 572 0 kH, respectively. These measurements and other 
similar ones to be presented later were made on a Solartron model 7061 
digital voltmeter that was capable of displaying seven significant-digits. 


The specific resistive strain gage used was manufactured by 
Micro-measurements Corporation and had a nominal resistance Ro of 120.0 
± 0 18 H. In addition, the gage factor S, was given as 2 095 ± 05 at 300°K 
while the transverse sensitivity factor k, was given as k, = 0.006. Finally u 
Poisson's ratios in equation (4.1.1) are given as v o =»0.285 and v = 0.33 and 
represent the Poisson's ratios of the factory gage calibration beam and the 
aluminum experimental test beam, respectively Refer to Appendix B for 

m 

further details. 
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In order to evaluate our resistive strain gage, the beam was initally loaded 
with a 200 gram mass using the loading technique depicted in Figure 25(b). 
This was done to take any ‘‘slack’’ out of the system. The gage resistance was 
then measured and the other resistances of the bridge circuit adjusted 
(tnmpots were used) so as to null the output. Additional load mass was then 
added to the beam with readings of V gut being taken as the excess load mass 
reached the values of 5. 10. 20. 30. 40, 50 and 100 grams. This procedure was 
repeated several times and a representative plot of data depicting the 
measured versus induced change in average beam strain is given in Figure 
26. In particular, note the regression polynomial calculated for this data. 
is seen, the calculated slone of the line deviates only slightly from the 
case of unity slope In addition, the y-axis crossing value can be shown to 
very closely correspond to the measurement uncertainty of the resistive strain 
gage. 

If we now conservatively assume that our Solartron volt/ohm meter is 
uncertain to within plus or minus one digit in the fourth decimal place (for a 
given scale factor) we may calculate our minimum change in s 
measurable with the resistive strain gage using the procedure outlined in 
Appendix B. With dV :u1 = ± 10" V, V m = 10.0 ±10 ~*V,R S = 572.0 ± 10" kQ and 
R t = 102.834 ± 10" MCI we find that our minimum detectable change in * 
average beam strain is Ae min = 0.21 nm/m, as measured with the resistive 
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strain gage. As this value also represents twice the absolute measurement 
uncertainty of the resistive strain gage (even with no applied load) it is seen 
that the y-axis crossing of the regression polynomial of Figure 26 very closely 
matches with the predicted value. We therefore see that our resistive strain 
gage provides a justifiable strain measurement standard. 


4.2 Dynamic Strain Measurements 


With confidence now gained in the performance of the resistive strain gage, 
the next task was to investigate the performance of the MDS/9-633 sensor. A 
discussion will be given in this section relating to the modal domain sensor's 
abilities to serve as a dynamic strain measurement device. Sections 4.4 ar 
4 5 will discuss static strain measurements. 

Dynamic strain was induced in the cantilever beam using a setup as shown in 
Figure 25(a). As is seen, the action of the vibration exciter (Bruel & Kjaertype 
4809) was transmitted to the beam via a push rod. The ball and socket joint 
was used so that the deflecting force on the beam remained normal to the 
surface and so that the shaft of the vibration exciter was not subjected to any 
off axis Stress. Notice, however, that this setup does not provide a way by 
which the force on the beam or the beam tip deflection can be measured as it 
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vibrates. All is known is that the beam can be caused to vibrate at the 


frequency of the signal input to the vibration exciter and that the amplitude of 
this vibration may be varied (in an uncalibrated way) by adjusting the input 
signal amplitude. For this reason it was important that the resistive strain 
gage behave properly so th~t the MDS/9-633 sensor readings could be 
compared to a reliable standard. 

It was decided that as the beam vibrated the amplitude of the sinusoidallv 
varying average beam strain, or peak average strain, would be the quantity to 
be measured For the modal domain sensor, we have from section 3.1 that the 
measured peak average strain is given as 

For the resistive strain gage we have from Appendix B that this same quantity 
is measured as 


4fl 5 r 1 — Vq ~ 
cwscpw* Vn s ^ 6 [ 1 - k fV 


(4 2.2) 


Note that in equations (4>2.1) and (4.2.2) we are considering that the beam is 
initially unloaded and that and V g correspond to the amplitudes of the 
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sinusoidal signals generated at the output of the MDS/9-633 and resistive 
strain gage sensors, respectively 

To avoid both resonance conditions and damping of the vibration induced into 
the cantilever beam it was decided that the beam would be excited at as far 
below its fundamental frequency of vibration as possible This frequency was 
determined by plucking the beam and allowing it to vibrate freely (with the 
push rod removed) and by then monitoring the output of the resistive strain 
gage sensor on a spectrum analyzer As seen in Figure 27 the fundamental 
frequency of vibration was measured to be 21.6 Hz. Since the vibration exciter 
could not be reliably operated below 10 Hz this was chosen as the excitation 
frequency for dynamic strain measurements 

To conduct the experiment the fiber launch conditions were adjusted until a 
rotating three lobe pattern was observed. This was often facilitated 
adjusting the tension in the length of fiber between the spools of the static 
strain control device, although no good explanation for this effect presently 
exists. The photodetector was then typically moved toward the fiber till the 
pattern just filled the open iris (a somewhat arbitrary act) and the Q-point was 
set using the PZT cylinder. 

•> 

It was observed that when the PZT cylinder was excited at its .resonance 
frequency of 38.5 kHz an appreciable small signal sinusoidal strain could be 
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imparted to the fiber, causing a small sinusoidal response at the photodetector 
output The Q-point of the system was then set by using an x-y positioner to 
position the closed iris over the point of the pattern that gave the highest 
amplitude and least distorted sinusoidal response. This proved to be the most 
reliable way of setting the quadrature point and Figure 28 shows a highly 
aligned MDS/9-633 output signal observed during such an alignment. The top 
trace is the modal domain signal and has a peak-to-peak voltage of 80 mV. 
The bottom trace is the input to the PZT cylinder and has a peak-to-peak 
voltage of five volts. Both traces are at 38 5 kHz. 

Conceivably, it should be possible to set the Q-point by locating the iris o._ 
the peak of a lobe and causing a static strain in the fiber so as to cause the 
pattern to rotate to the point where the iris is located half way between the 
maximum and minimum of a lobe. In practice this proved to be difficult, 
though, as it was often difficult to align the input in such a way that the three 
lobe pattern rotated properly through more than about twenty to thirty 
degrees. This was fine for small signal work but made alignment of the Q-point 
difficult using the “bulk effect" technique just descnbed. Therefore only the 
“small signal" method of setting the Q-point, using the excited PZT cylinder, 
was used. 
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With the O-point properly set. data was taken by recording the peak voltage 
out of each of the sensors on the cantilever beam, for varying amplitudes of 
the signal input to the vibration exciter, at the 10 Hz peak observed on a 
Spectrum analyzer. Recall that for all dynamic measurements the beam 
excitation is at 10 Hz. Plots of some of this data, converted to strain 
measurements using equations (4 2.1) and (4.2.2) are given in Figures 29, 30 
and 35 with several Polaroid photos corresponding to the data of Figure 30 
being shown in Figures 31-34. These shall be discussed shortly. 

For the data of Figures 29. 30, and 35, the load resistance was 50 kX2 and tb*» 
input to the Wheatstone bridge V n was taken from a nine volt battery. Tnc. 
value of a was taken to be x = 0.78 in equations (4.2.1) and (4.2.2), as per the 
discussion of section 2.4, although this value will be seen to change when 
static strain measurements are undertaken. Recall also that Ap = 1.2x10 4 m - ' 
for the MDS/9-633 sensor. In addition, so as to make the data sets 
independant, the pattern and quadrature point were readjusted before eacn 
set of data was collected. For the data sets of Figure 29, 30 and 35 the Q-point 
voltages V 0 for the MDS/S-633 sensor were 1.2 V, 0.55 V and 0.35 V, 
respectively. Since the sensitivity of the MDS/9-633 sensor directly depends 
on the Q-point value of optical power through the spatial filter aperture, it is 
easy to see why the y-axis crossing is at a much higher value of strain in 

■s 

Figure 35 than it is in Figure 29, Note, however, the linearity and near unity 
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slope of all data in Figures 29 30 and 35. Especially note that the y.-axis 
crossings of the extrapolated data in each of these figures agrees well with the 
minimum expected detectable average beam strain calculated in section 3.1 

Figure 31 shows a spectrum analyzer trace for both the MDS/9-633 sensor (top 
trace) and the resistive strain gage (bottom trace) at a point where the beam 
excitation is small. Note the output of each sensor is precisely at 10Hz and 
that the MDS/9-633 sensor has a much higher signal to noise ratio than the 
resistive strain gage sensor, thus indicating that the modal domain sensor w> 
be able to provide strain measurement data long after the resistive strain gage 
has reached its minimum. Also, the spectrum analyzer gave the amplitude of 
each detected signal in dBV RMS so that peak voltage was calculated using 
the conversion 

_ . V 4SV 9KIS ■ 

V p = N 2 10 20 (42.3) 

Furthermore, the data of Figure 31 corresponds to the first blackened data 

point in Figure 30 and is seen to be at the low end of the usable range of the 

resistive strain gage. Similarly, the polaroids of Figure 32 correspond to tht. 

second blackened data point of Figure 30 and are the result of a larger beam 

excitation. All data shown in these photographs were averaged eight times for 
% 

visual purposes only. Finally, Figures 33 and 34 show the vibration exciter 
input signal and the noise floors of the strain sensors at quadrature with no 


* 


J 

4 
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excitation, respectively The bottom trace in Figure 34 corresponds, at 10 Hz, 
to the noise {-104 5 dBV RMS) of the MDS/9-633 sensor, while the top trace 
corresponds to the noise (-86.0 dBV RMS) of the resistive strain gage sensor. 
Using equations (4.2.1). (4 2.2) and (4.2 3) it is easily seen that these noise 
figures correspond to minimum detectable average beam strains of 
L = 2 8x10'* m/m and c ^ 0.1x10’* m/m for the MDS/9-633 and resistive 
strain gage sensors, respectively. This is in very good agreement with the 
predicted minimum detectable strains for each sensor and lends validity to the 
assumption of negligible excess environmental noise. 

There are two final comments to make before moving to the next section. First, 
it should be mentioned that the maximum tip deflection never exceeded more 
than a few millimeters for all dynamic measurements. Second, although no 
values for beam displacement have been given, this can be calculated if the 
average strain in the beam is known using the relation 



where d is the beam tip displacement from its undisturbed position, L is 
length of the cantilever beam and h is the beam ihickness [1]. For the modal 
domain sensor, this gives us 
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where the beam is assumed to be initially unloaded. 


4.3 Extraneous Dynamic Modulation 


At this point it is reasonable to wonder if there are any effects other than 
rotation of the three lobe intensity distribution which give rise to the signals 
discussed in the previous section. In order to investigate this, a properly 
rotating three lobe pattern was obtained and the detector positioned, with the 
iris fully open, so that the pattern Tilled the area of the detector face. The b 
was then caused to vibrate at 10 Hz using the vibration exciter, with data taken 
as peak voltage readings at the output of both the photodetector and the 
resistive strain gage sensor circuits. Figure 36 shows this data taken for a 
senes of varying signal amplitudes into the vibration exciter. The data is 
plotted as millivolts of extraneous modulation out of the photodetector versus 
resistive strain gage detected strain. 

Figure 37 shows Polaroid photos of the spectrum analyzer traces for large (top 
photo) and moderate (bottom photo) beam excitations. These correspond to 
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their respective blackened data points in Figure 36. Similarly, Figure 38 shows 
a Polaroid photo of the spectrum analyzer trace for small beam excitation and 
corresponds to the blackened data point nearest the origin in Figure 36. For 
all photos in Figures 37 and 38 the top trace shows the output of the resistive 
strain gage sensor while the bottom trace shows the extraneous modulation 
output of the photodetector. Figure 39 shows the noise floors, with no 
excitation, for the measurements of extraneous modulation. Again, the top 
trace corresponds to the noise of the resistive strain gage (-88 dBV RMS) while 
the bottom trace corresponds to the noise of the photodetector output (-99.5 
dBV RMS) 

In a typical measurement situation, recall, the iris will be closed, leaving i ~- 
a 1.0 mm aperture. This means that with the diameter of the iris open bein^, 
1 2 cm, only one 144 m of the pattern is spatially filtered in a normal 
measurement situation. This spatial filtering also takes place at approximately 
the average intensity of the pattern, at quadrature, so that the extraneous 
modulation figures shown in Figure 36 will be divided by 144 for measurement 
situations such as those described in section 4.2. For the largest extraneous 
modulation reading of Figure 36 (-62.1 cfBV RMS or 1.11 mV), the 
corresponding value after spatial Filtering is 7.7^V (-105.26 dBV RMS). This is 
easily seen to be well within the noise of the optical sensor output, even as 
shown in Figure 34. We thus conclude that the results presented in section 4 2 
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are entirely due to three lobe rotation and that extraneous modulation effects 
are entirely negligible. 

As a final comment for this section, it has been suggested that, though small, 
the extraneous modulation effects are due to bend loss. If this were true, 
however, a peak at 20 Hz would be found due to bend loss occuring when the 
fiber both goes into tension and compression as the beam vibrates, but no 
such -peak was observed. A more likely reason for the extraneous 10 Hz 
modulation is that some 10 Hz excitation is being transmitted to the floating 
table as the cantilever beam is forced to vibrate. Because the 10 H-* 
extraneous modulation figures are small, though, it is evident that the flc 
table is properly performing its function of excess vibration suppression, as 
expected, even though a vibrating structure is attached rigidly to it 


4.4 Static Strain Measurements (First Attempt) 


It was also of interest :o study the static strain measurement capabilities of tr.- 
MDS/9-633 sensor. As with the resistive strain gage evaluation measurements, 
the cantilever beam was loaded as shown in Figure 25(b), with massSs 
periodically added to the basket in excess of an initial load mass of 200 grams. 
Now, however, the change in average beam strain, as measured by the 
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MDS/9-633 sensor, was to be compared to the known induced change in 
average beam strain caused by varying the load mass. Recall that the induced 
change in average beam strain is given as 


3mgL 
Ebh 2 


(4.4.1) 


where m is the change in load mass. Further recall that this same quantity, 
as measured by the MDS/9-633 sensor is given by 


_ AV 
Z ' V Q A0 L, i 


(4.4.2) 


as was shown is section 3.1 and where A/? = 1 2x10* m-\ 


As before, the input to the optical fiber was adjusted until a properly rotating 
three lobe pattern was observed and the Q-point was adjusted, with the initial 
load mass in place, u< >g the excited PZT cylinder. Mass was then added to 
the basket and AV measured (this time through use of a strip chart recorder) 
as the excess load mass reached 50, 100, 150, 200, 250 and 300 grams. It wa 
observed, however, during early static loading tests that there was 
appreciable pattern drift over a period of a minute or so. It is believed that this 
drift is attributable to the static strain control mechanism and due in part to 
jacket shear for large amounts of induced static strain (recall it was mentioned 
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earlier that obtaining a rotating three lobe pattern was often facilitated by 
tensioning the fiber using the static strain control device). As the need fo. 
loading, unloading and settling the basket make the procedure for taking static 
strain data much more lengthy (1-2 minutes) than the procedure for taking 
dynamic strain data (<30 seconds) it was necessary to fully release the 
tension of the fiber in the static strain control device when static strain 
measurements were to be attempted 

In addition, preliminary plots of measured versus induced strain data gave 
regression lines whose slope was often noticeably different from unity. This 
was believed to be due to the physical parameters leading to the calculation 
of a, in equation (4.4 2), not being precisely known for the fiber being used. 

In an attempt to directly measure 7 le launch conditions were adjusted until 
a three lobe pattern was obtained whic' exhibited proper rotation over greater 
than 300°. With the fiber between the spools of the static strain control device 
just taut, the tension in the fiber was increased by turning the micromet?- 
screw until the pattern rotated 120° (lobe peak to tobe peak). The change in 
the micrometer setting was then recorded and the tension in the fiber quickly 
released. This procedure was repeated several times and it was determined 
that the micrometer was advanced an average of 691/im for 120° rotation. The 
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change in length AL, :o of the fiber between the spools was therefore 
approximately 

= 564.75 fim , 

where r> = 7 62 cm and r 2 = 9 32 cm reprpsent, respectively, the radius of the 
rotating spool and the length of the push rod. This indicates that a change in 
length of 141 19 /im is needed to cause a rotation from quadrature to the 
nearest lobe peak (a 30° rotation). 

From equation (3 1 11) we can solve for the change in fiber length required to 
cause pattern rotation from quadrature to the nearest lobe peak as 


691t 

= rjSi n { — 


AL 30 


2 * 

4 A/J a 


(4 4 


where, recall. AP = 1 2x10 4 m~' for the MDS/9-633 sensor. With 
AL :a = 141.19 nm we can thus solve for x to find a = 0 92. on average, for our 
fiber. 



New data was now taken for x=0.92, with R L = 50/tii. See Figures 40 and 41 
for plots of this data and note the near unity slope for each of these figures. 
In addition, the quadrature 'point voltage V a was 25.8 mV and 23.1 mV for the 
data of Figures 40 and 41. respectively. 
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It is interesting to note that the highest point on the regression line of Figure 
40 occurs at Ac f = 104 5 umjm or AL, = Ael^ = 60.1 um. This indicates that 
the response of the MDS/9-633 sensor stays linear for AV, in equation (4.4.2), 
approximately as large as 0 66t/~ The dynamic range for this linear response 
region can thus be found, in decibels, as 


OR = 10 log 


APr 


DR = 10 log 


0 66P o 


■y(4 qB{2P Q - ! 0 


8 k s TB \± 


(4 4.4) 


where AP mn was given in equation (2 3.10) and where P Q is the optical powe. 
through the spatial filter aperture at quadrature. Recall that the voltage at the 
output of the photodetection circuit is directly related to the optical power 
through the spatial filter by the photo diode responsivity & and the load 
resistance R L . Using the system parameters used in section 3.1 to determine 
the minimum detectable strain of the MDS/9-633 sensor we thus find that the 
dynamic range of this sensor, in the lirear region, is approximately 52 dB. 
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4.5 Static Strain Measurements (Second Attempt) 


Due to the drift problems encountered using the static strain control device 
shown in Figure 24 it was decided that this device would be redesigned and 
additional static strain measurements attempted. Shown in Figure 42 is the 
resulting new static strain control device. As seen, a length of optical fiber is 
attached to two axially aligned mounts, one of which exhibits axial translation 
as the motor shaft turns The motor used was a 200 step per revolution 
stepping motor and connected between the motor shaft and the movab'- 
mount was a "2 threaded rod with fifty-six threads per inch As the ga^. 
length of fiber between the two mounts was measured to be L, — 0 522 m, then 
for every step applied to the motor the change in induced average strain in the 
fiber was Ae = 4 345 fimjm 

This new static strain control device provided several advantages. First, the 
induced strain in the fiber was purely axial, as opposed to the axial and 
possibly transverse strains induced in the floer wound on the spools of the . 
static strain control device. Second, there was a very simple relation between 

the number of steps applied to the motor and the strain induced in the fiber, 

* 

as discussed above. Most importantly, though, was that there was no jacket 
shear in the fiber between the stationary and movable mounts This effect was 
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eliminated by first stripping the fiber of its |acket at the points where it was 
clamped and by then bonding directly to the bare fiber a length of stainless 
steel tubing. Due to the more rigid nature of the stainless steel tubing and the 
epoxy used to bond it to the bare fiber the, induced strain was more faithfully 
transmitted to the optical f'ber, with no relaxation of the tension occuring due 
to elastic stretching of the polymer jacket. Indeed, it was observed that 
negligible pattern drift occurred over a period of several days with the fiber 
under significant tension. The only disadvantage to using tne new static strain 
control device was that the fiber between the clamps had a tendency to break 
just at the point where the fiber exited the stainless steel tubing This is 
believed to be due to slight axial misalignments of the fiber mounts causing 
high stress concentrations in the regions of these breaks. 

As before, the input to the fiber was aligned such that a three lobed intensity 
distribution was present at the output that exhibited pure rotation. Pains were 
taken to ensure that the alignment was such that this rotation would occur over 
most of a full revolution Then, using a linear x-y positioning stage the spatial 
filter aperture was positioned over the peak of a lobe, with steps then being 
applied to the motor to cause the pattern to rotate to the point where the 
spatial filter monitored the point located half way between the maximum and 
minimum intensity points of a lobe. Due to the discrete nature of the increases 
in applied strain, though, the proper quadrature point was often overshot. 
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requiring final alignment using the x-y positioning stage With the operating 
point thus properly set. strain was measured as additional steps were given to 
the motor 

As before, the strain measured by the MDS/9-633 sensor followed equation 
(4.4 2). In addition, as before it was noticed that plots of preliminary measured 
versus induced strain data gave a linear regression line with slope different 
from unity when x was taken as x = 0.78. By fitting only the first point of . 
preliminary data to a regression line of unity slope it was determined that 
x = 0 86 was a more acceptable value. Assuming a = 0.86 additional data was 
taken and is shown in Figures 43 and 44 

Two final comments are necessary. First, the method by which a 
determined for the data plotted in Figures 43 and 44 is admittedly questionable 
due to the preliminary data used to determine a possibly not being 
independant from the data subsequently taken and presented There simply 
was no way in which a reliable and independant determination of a could be 
determined with the setup used. Indeed, an independant determination of the 
photoelastic coefficients used to calculate a is likely to be an involved task in 
itself. It is believed, however. t**at photoelastic effects must be considered as 
the fiber is subjected to strain and that the analysis leading to equation (4.4.2) 
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is correct The value of x. however, and its constituent terms is very much 






l 


t 

i 


open to question. 

The other comment necessary is that this second round of static strain 
measurement data was taken using a UDT model 81 Optometer This 
instrument directly displayed the optical power passing through the spatial 
filter aperture and incident on the face of the detector head. As per tt" 
discussion in section 3 1 leading to the result of equation (3.1.14). however. 
AP and P c may be directly substituted for AV and V Q in equation (4.4.2). The 
value of P 3 was 0.13 uW and 0.16 uW for the data plotted in Figures 43 and 44, 
respectively. 
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Figure 23. Experimental MDS/9-633 system 
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Figure 24. First static strain control davica 
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Figure 25. loading schemes of the experiments MD S/3-633 sensor system 
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Figure 28. Highly aligned MDS/9-633 sensor output (sc® tost for details) 
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Figura 30. MDS/9-C33 dynamic stialn m*a*ur«m«nt data (trial 2) 


4.0 MDS/3^33 SENSOR EVALUATION 


103 





ORIGINAL PAGE IS 

cm EQQ8, QlTAJJ-m 



4.0 MDS/9-633 SENSOR EVALUATION 


105 



Figure 33. Input signel of ih® vibration exciter 
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Figurs 33. Extraneous dynamic modulation signals for small beam excitation 
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Figure 40. MDS/9-433 static strain manure merits (first attempt; trial 1) 
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Figure 41. MOS/9-633 static strain measurements (first attempt; trial 2) 
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Flguro 42. R«d«slgned static strain control dovico 
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Figure 43. MDSS3433 stave it rain measurements (second attempt: trial 1) 
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Figure 44. MDS/S-633 static strain measurements (second attempt; trial 2) 
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5.0 CONCLUSION 


Modal domain sensors using 9 um core diameter fiber at source wavelengths 
of 633 nm and 850 nm have been thoroughly described. The mode content of 
the most commonly observed intensity distributions at the fiber/sensor 
outputs, for each wavelength, have been discussed and the resulting str 
modulated behavior of these intensity distributions has been proven from first 
principles. In addition, it has been shown how this strain modulated behavior 
may be used as a strain detection mechanism and suggestions have been 
made concerning how micro-optical components can be used to construct a 
MDS/9-850 sensor. 

The linear and highly sensitive nature of modal domain sensors has also been 
predicted and experimentally verified for small amounts of strain and the 
dynamic range of a typical modal domain sensor has been calculated. 
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Furthermore the many advantages and disadvantages of the various modal 
domain techniques have been addressed in great detail. 

It is hoped that future work with modal domain sensors will involve a complete 
experimental evaluation of a MDS/9-850 sensor and that further investigations 
will result in the realization of additional in-line micro-optical 
injection/detection techniques for modal domain sensor use In addition, an 
independant study into the values of the individual photoelastic coefficients of 
the optical fiber used in the previously described experiments should be 
undertaken with a parallel study undertaken to determine the radial changes 
of the core index of refraction as a function of strain and how this will affect the 
propagating modes. 

Finally, it was assumed in section 4 2 that for quasi-static induced beam 
vibrations at frequencies well below the fundamental natural frequency of 
vibration, with small peak beam tip displacements, that the cantilever beam 
stress-strain relations held just as if the beam was statically loaded, with on 1 '' 
a sinusoidal multiplicative factor appearing for the quasi-static strain case 
was also assumed that only tne fundamental mode of vibration was excited. 
Although the corresponding nature of the experimental results and the 
theoretical predictions as to the behavior of the MDS/9-633 sensor lend 
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validity to these two assumptions, further investigation may be necessary to 
finally prove that these assumptions are valid. 
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APPENDIX A: Strain Characteristics of Loaded 
Cantilever Beams 


One of the simplest and most fundamental ways in which a newly developed 
strain sensor may be evaluated is by attaching it to the surface of a cantilever 
beam which is placed under load. If the beam is constructed from some 
homogeneous material, such as aluminum or steel then the strain 
characteristics of the beam are quite easily described as a function of load, 
thus making the loaded cantilever beam a convenient strain standard. 
Consider for instance the cantilever beam of Figure 45. Assuming the 
is made of a homogeneous material, the stress at any location on or within 
beam is given by 


a 


M g z 

l 


[Pa] . 


(41) 
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where 2 is the distance above or below the neutral axis. M, = F x i is the 
bending moment and I is the moment of inertia given as 


h 2 b 

12 


(A 2) 


If we now consider only the surface which is in tension when the beam is 
loaded we find that the surface stress a, ur is 


Ffl/i/2) = e Ff_ 
hb 2 l 12 h 2 b 


(A.31 


Recalling that strain z is related to stress by the relation e = c t/E , where E is 
the modulus of elasticity, we see that the strain along the tensed surface of the 
cantilever beam is 


/ e\ . a a \ 

£,(-0 = 5 — • 

Eh 2 b 

Note that this result indicates that the strain along the surface of the beam 
increases linearly as a function of (, being zero at the free end and be.. 
maximum at the clamped end. The average strain in the beam thus occurs at 
£ — LI 2. Also note that by changing the sign of the strain relation in equation 
(A. 4) one can describe the strain induced in the beam surface being 
compressed. The negative sign simply indicates the compressional force. 
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Often, though, one is interested in the change in average strain in a situation 
where the loading force is changed by varying masses attached to the beam. 
We will call this change in average strain Ac which is simply given by 


Ac = 


3mgL 

Ebh 2 


(*•5) 


where m is the change in mass, g is the acceleration due to gravity and the 
fact that average strain occurs at £ = L/2 has been used. This quantity can 
easily be measured by a resistive strain gage placed at the geometric centc. 
of the beam surface [12] 


A A Average Strain induced in an Optical Fiber 


In the evaluation of fiber optic strain sensors it is often desirable to loop a 
sensing length of fiber along a cantilever beam as indicated in Figure 4P 
Since the output of an optical fiber strain sensor is actually a measure of 
average strain in the fiber sensor length (more commonly called the gage 
length) we can calculate the average change in strain. Ac, , to which the fiber 
is subjected along its path by evaluating 
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Si/ — 


A L 




1 

4 


'L, 


£, Cfif 


where AL is the total elongation of the fiber path and L, is the gage length. 
Using equation (A 4) we can evaluate AL as follows - 


A L 


2 z x df z x df 

Jab *ac 


A L 


2 1*" d£ - 2r f 2 

e/i‘h -o 


6f(L, — rcoso) 

sin <p dQ 

Eh‘b 


A L = 


6mg 

Eh 2 b 


i L 2 - L, 2 ) 


12mgr 

Eh z b 



AL = -^-(L 2 -L 2 -2L,r-r 2 ) . {A. 7) 

Eh 2 b 

Note that a simple transformation from rectangular to polar coordinates was 
necessary to integrate the strain effects in the curved portion of the fiber. Now, 
by dividing equation (A.7) by the fiber gage length and recalling that m is the 
change in load mass we find that Ac, is 
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where 


ACf 


6 mg 
Eh 2 bL f 


(L 2 -L^ + 2L x r-r 2 ) 


(A 8) 


L, = 2(L - L,) 4- rrr 


(A 9) 


Notice in equation (A 8) that if r=Li so that the curved portion of fiber just 
touches the free end of the cantilever beam then 


Sc, = 


6 mg i_ 2 
Eh 2 b L < 


(A 10) 


which, if Lf =s 2 L, is approximately the same result expressed in equation (A.5). 
More precisely, if we choose Li, L and r properly, it is possible to make the 
results of equations (A.5) and (A.8) express the same quantity. Assuming L 
and Li to be fixed, we can solve for r as follows: 


6 mg 
Eh 2 bL, 


(L* — L 2 + 2rL x 


r 7 ) 


ZmgL 

Ebh 2 


L 2 -L? + 2rL,-r 2 = 2(L - L,) + nr) 

f2 ~ r (zL x — ~~ ^ — (L X L — L 2 ) = 0 
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r= + {‘•.'■(’-t) + (^) 2 } 


(4.11) 


These last considerations are interesting in that they demonstrate that by 
using simply loaded cantilever beams an optical fiber strain sensor can be 
subjected to Known strain effects while having its performance evaluated 
simultaneously against more commonly used strain sensors such as resistive 
strain gages. Finally, it should be mentioned that due to the relatively 
insensitive nature of both resistive strain gages and optical fibers to 
transverse strain, no considerations in the previous formulations have been 
given to transverse beam contractions [13]. 


1 / 
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Figure AS. Cantilever beam geometry 
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Figure 46. Optical fiber and resistive strain gage placement parameters 
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APPENDIX B: The Resistive Strain Gage 


Lord Kelvin first observed that a change in strain imposed on a wire is 
accompanied by a change in resistance Art of the wire The relationship 
between resistance change Art and strain t can be derived by considering a 
uniform conductor of length L, cross-sectional area A . and resistivity p. The 
resistance R of such a conductor is given by 

* - • > 8 ’> 
Differentiating equation (B.1) and dividing by R gives 

dR _ d P , dL dA 

~R~ = ~ + ~L ‘ (fl - J 


However, 
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** - * «*fc*W* t«Hi 


R : 






under elastic deformations where v is the Poisson's ratio of the conductor 
material The strain sensitivity S A of the conductor is defined as 


S* = 


dR/R 


Substituting equations (B 2) and (B.3) into equation (5 4) gives us 


S* = 




The last three terras of equation (B.5) are due to dimensional changes in th° 
conductor The first term is due to changes of resistivity with strain. The last 
term is usually neglected for elastic strains since it is small I < 0 It) compared 
to the other terms ( — 1.60). The derivation of (B 5) is modified slightly for 
large strains since the conductor undergoes plastic deformation [13]. 
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B.1 Strain Gage Calibration 


In practice equation (8 5) is not very useful since it is only valid for a single 
filament conductor and because the first factor which describes the change of 
resistivity with respect to strain is often unknown. Typically the conductor of 
a resistive strain gage is formed into a grid to reduce its size as shown in 
Figure 47 This causes the gage to exhibit sensitivity to both axial and 
transverse strain The response of a surface mounted gage that is subjected 
to a axial stra.n a transverse strain e. and a shearing strain can be 
expressed by 


dR 

R 


^a £ a + S s ) at . 


(B. 6) 


where S,. S,andS, are the sensitivities of the gage to axial transverse and 
shearing strains, respectively In general, the gage sensitivity to sheanng 
strain is small and thus neglected so that 



(3 7) 


where k t = S~S t is defined as the transverse sensitivity factor for the gage. 
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Strain gage manufacturers provide the transverse sensitivity factor k, and a 
calibration constant known as the gage factor for each gage. The gage factor 
S, represents the calibration constant for a batch or lot of gages and is 
determined by testing sample gages drawn from a lot of gages in a given 
production run Resistance change dR experienced by a gage is related to the 
gage factor and the axial strain by 


dR 

R 


= V, 


(6 8 ) 


The stress field in the calibration beam used for the determination of S, is 
always unaxial. therefore the gage is subjected to a biaxial strain of 

«r = “ v 0 £ a . CS S) 

where v 0 = 0.285 is Poisson's ratio for the calibration beam material 
Substituting equation (B 9) into equation (B 7) and companng with equation 
(B.8) yields 


■Sg — S a (1 k(V o) 


or 



(S.10) 
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Therefore when the strain gage is used to measure unaxial strain in a material 
other than the calibration beam we see that 


dR 

R 


S g 

1 _ k , (h + Vi) = 


1 - Vo 


(1 


V) 


or 


fa = 


dR ( 1 ~ Vo \ 1 

r { i - v y 


(S 11) 


where v is the Poisson s ratio of the structure to which the strain gage is 
attached. It should be evident from equation (B.11) that the calibration 
procedure used to determine k,andS 7 allows transverse strain effects to be 
“calibrated out’ of the actual strain gage measurement [13] Also, it should 
be evident that the change in average strain in a loaded cantilever beam as 
described in Appendix A could be measured using a resistive strain gage 
placed in the geometric center of the beam by using the relation 


Ae = 


A R ( 1 ~ V o \ 1 
R { 1 - k t v ) S g 


( 8 . 12 ) 


where R is the initial gage resistance (often measured), v is Poisson's ratio of 
the test material and A R is the measured change in gage resistance 
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B.2 The Use of Wheatstone Bridges to Measure A R 


In order to enable small changes in strain and therefore small changes in 
resistance to be measured in a resistive strain gage, a Wheatstone bridge 
circuit with an amplifier output stage is often used. A constant-voltage 
excitation Wheatstone bridge circuit is shown in Figure 43. For this bridge 
circuit, the error voltage Eo is given by 




IR, ~ R : )(R 3 - * 4 ) 


(8.13) 


Equation (B 13) indicates that the initial error voltage will vanish (Eo = 0) if 

R^R 3 — ^ 2^4 (8.14) 


When equation (B.14) is satisfied, the bridge is said to be balanced. This 
means that the small unbalanced voltage caused by a change in resistance of 
any arm of the bridge is measured from a zaro or near-zero condition. This 
small signal can easily be amplified to significant levels for recording. 


If a single resistive strain gage is placed in the bridge circuit a indicated in 
Figure 48, a voltage AEo developes at the input to the amplifier stage when the 
nominal resistance of the gage R is changed by an amount AR Such changes 


APPENDIX B The Resistive Strain Gage 


134 



in resistance are in general due to strain or temperature variations: however 
in the following analysis the change in A R due to temperature effects will be 
assumed negligible. 

If the fixed resistors R 2 Ri and Rj are chosen to have the same resistance 
value as the nominal resistance R of the resistive strain gage, we have that 

A£ 0 = X ^" (1 ~ n)V ‘" ' (S 15) 


"/here 


•7 = 


SRIR 
AR/R -r 2 


(B 16 ) 


Substituting equation (B 16) into equation (B.15) we find that 


0 A R \ SRIR -r 2 I “ 


SE, 


= _L M. ( 2 '\ v 

4 fi \ ARIR + 2 ) ,r 


_ 1 SR ( * SR \ u 

AE„ = T — (t - ■ 
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where the binomial expansion (1 i- a) 1 =» (1 - a) has been used since 
typically A R << R Continuing. 

A£o = X ^V in - -^—v in . (B. 17) 

4 R 8 R 2 

Since the second term cf equation (B 17) is very much smaller that the first, 
we have 


- -J- • (B.18) 

The output of the inverting amplifier is *hen given by the relation 

R? 1 Rz ,\p 

v ~- -rT ae »= -T sT f v " • |S - 19) 

where Rs is chosen large enough so as not to cause significant loading of the 
bridge network and where a near-zero initial condition is assumed for AEa 
(i.e., the bridge is assumed to be initially balanced) 114]. 

By substituting equation (B.12) into equation (B 19) and then rearranging 
terms we find that the change in average strain induced in a resistive strain 
gage is given by 
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Ac = 


1 - k ,\ 0 


( 8 . 20 ) 


Vpat _f*5_ 

V,n S g R 6 


1 — k t v 


Due to the linear nature of resistive strain gages. Wheatstone bridges and 
simple inverting amplifiers, the peak change in average strain as induced in a 
resistive strain gage by a sinusoidally varying strain field is given by 


Al p = 


_^p 

V. n 


5 


1 - V 0 
i - v 


( 8 . 21 ) 


where the negative sign has been dropped in equation (B 21) since V a , which 
represents the amplitude of the sinusoidally varying output voltage, is 
independant of phase Considering once again the case of measuring the 
change in average strain in a cantilever beam we see that equation (B.20) will 
be most useful when the beam is subjected to a static load and that equation 
(B21; will be most useful when the beam is forced to vibrate under the 
influence of a sinusoidally varying load. 


Finally, in order to investigate the minimum detectable strain for a resistive 
strain gage used in a detection network such as that of Figure -48, we 
differentiate equation (B.21) and divide by Ai„ to obtain 


c/A£ p dV p dR s dV in dS g dR s 

~A^~ = ”v£" + ~R~S K, ^ «T 


( 8 . 22 ) 
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'r^ '•« >• j y, % ^ 


i 


where in equation (B 21) k \o and v are assumed to be constant and where 
dR s , dR s and dS g represent manufacturer stated tolerances on the values Rs, 
R$ and S J( respectively. Also. dV p and dV m represent the inaccuracies of 
measuring V p and V n on. say. digital voltmeters, oscilloscopes or spectrum 
analizers By then setting the ratio c/AEj/AZ, equal to unity and solving for 
V. we will have a minimum value of V B needed for reliable measurements 
which when substituted into equation (B.21) will give the minimum reliable 
value of A I p Minimum detectable strain values on the order of 10~* are not 
uncommon for resistive strain gages used in typical measurement situations, 
with higher sensitivities expected in cortrcled laboratory environments [13] 


t- 

t 

f S' 
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Figure 47. A resistive strain gage 
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( APPENDIX C: Bessel Function Numerical 
Expansions 


1 

i 

The following numerical expansions for Jo(x) and Ji(x) can be used to aid in 
plotting modal intensity patterns [15]. 


where 

A = 2.2499997 
8 = 1 .2656208 


i 
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C = 0 3163866 


D = 0.0444479 
E = 0.0039444 
F = 0.0002100 

and 


where 


G = 0.56249985 
H = 0 21093573 
/ = 0.03954289 
J = 0.00443319 
K — 0 00031761 
L = 0 00001109 

These expansions give eight significant digits for — 3 < x <, 3. For reasonable 
departures of I x | ^3 these expansions also give very good results for plot* 
purposes. 
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APPENDIX D: FORTRAN Programs 


The following FORTRAN programs are useful for generating the data needed 
to create the line and surface plots of the various intensity distributions of 
interest in modal domain sensing. 
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CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


THIS PROGRAM GENERATES THE VALUES WHICH WHEN PLOTTED RESULT 
IN THE LINE PLOTS WHICH DESCRIBE THE TWO LOBE PATTERN OF 


INTEREST IN MODAL DOMAIN SENSING 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCI 



c 

REAL JO(IOO). Jl<100),XVAL(100),I,J,K,l 
C 

A = 0.56249985 
B = 0.21093573 
C = 0.03954289 
D = 0.00443319 
E = 0 . 000317ol 
F = 0.00001109 
G = 2 . 2499997 
H = 1 . 2656208 
I = 0.3163866 
J = 0.0444479 
K = 0 . 0039 AAA 
L = 0.0002100 
C 

RBAR = -A. 59 
C 

DO 10 N = 1,100 
C 

RBAR = RBAR +0.09 

XVAL(N) = RBAR 

XI = 62SKR3AR/3.0 

XO = . 400XRBAR/3 0 
C 

JOIN) = 1.0 - G*(X0**2) + Hx(X0xx4) - lx(X0xx6) 

$ + JX(XOXXS) - Kx(XOxxiQ) + LX(X0XX12) 

C 

J1(N) = (0.5 - AX(X1XX2) + BX(XIXXA) - C*(Xlxx6 ) 

$ + DX(XlxxB) - £*(X1XX10) - FXCX1XX12)) 

$ X3.0XX1 

C 

10 CONTINUE 
C 

DO 20 N = 1.100 
C 

OUT = <J0(N)xx2) 

$ ♦ (Jl(N)xx2) 

$ ♦ (2.0XJ0(N)XJ1(N)) 

C 

WRITEt 08, *) XVALCN'.OUT 
C 

20 CONTINUE 
C 

STOP 

END 
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CCCCCCCCCf 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


:ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THIS PROGRAM CALCULATES THE VALUES WHICH WHEN PLOTED CREATE 
THE THREE LOBE PATTERN OF INTEREST IN MODAL DOMAIN SENSING. 
NOTE THAT THIS ROUTINE EXHIBITS HIGHEST ACCURACY WHEN THE 
ARGUMENTS OF THE BESSEL FUNCTIONS ARE IN THE RANGE (-3. SI. 

(THIS ROUTINE USES THE LP11 AND LP21 EVEN AND ODD MODES!!!) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
XX XX 

xxx go DOWN 38 LINES TO CHANGE PLOT PARAMETERS'" XXX 

XX XX 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 

IMPLICIT DOUBLE PRECISION CA-Z) 

INTEGER N, M 


A = 3.56249935D0 
B = 0.21093573DO 
C = 0 . 03954289DO 
D = 0. Q0443319D0 
E = 0.00031761D0 
F = 3 . 0Q001 1 09DO 
G = 2 . 2999997 0D0 
H = 1 . Z656208CD0 
1 = 3 316 3866CD0 
J = 3.34A44790DO 
K = 3 . 30394990D0 
L = 3 . 00021 000D0 


PI = 3.191592659D0 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


CHANGE THESE PARAMETERS TO MODIFY THE PLOTS COMPLETELY. 

A1 AND A2 ARE THE AMPLITUDE TERMS. DTHETA IS THE RELATIVE 
PHASE TERM AND Z IS SOME DISTANCE ALONG THE FIBER. 

C PHOTOELASTICITY IS NOT CONSIDERED IN THIS PROGRAM.) 


A1 

A2 

DTHETA 

Z 


1.000 
1 .ODO 
O.TODO 
O.OOD-4 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

B1 = 1.2004 

c 

X = -4.750100 

Y = -4.7501D0 
C 

DO 10 N = 1,37 

X = X + 0.2500 
DO 20 M = 1,37 

Y = Y + 0.25D0 
C 

RSQR = (X**2) + CYXX2) 

RBAR = DSQRT ( RSQR) 

PHI = DATAN2CX, Y) 

XI = ,688D0*RBAR/3.0D0 
X2 = . 9 06 DQxRBAR/3 . 000 


C 

C 

C 

C 


J02 


1.000 - G*(X2**2) + HX(X2XX4) - I*CX2**6) 

+ J *( X2**S ) - K*CX2**10> ♦ L«(X2**12) 


J12 = (0.500 - A*(X2**2) ♦ Bx(X2XX4) - C*(X2**6) 

+ Dx<X2xx8) - EXCX2XX10) ♦ Fx(X2xxl2)) 
X3.000XX2 

J2 = (J12X2.0D0/CX2X3.0D0)) - J02 

J1 = (0.5D0 - AX(X1X*2) ♦ Bx(Xlxx^) - Cx(Xlx*6) 

+ DX(XIXXS) - Ex(XlxxlO) + Fx(xix<12) ) 
X3.0D0XX1 

INTENS = 

( ( AlxJl )XX2) 

♦ (CA2XJ2JXX2) 

+ (2 . ODOXA1XA2XJ1XJ2XDCOS(3 . ODOXPHI 

- B1XZ ♦ DTHETA 


WRITE(8,X) X,Y, INTENS 
20 CONTINUE 

Y = -4.750100 


10 CONTINUE 
C 

STOP 

END 
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Introduction to Smart Structures 

Dr. Richard O. Claus 

Director Fiber Eectro-Optics Research Center 
Protessor. Electrical Engineering 

Smart," ''intelligent” or 'sense-able” structures are structures which contain their own 
sensors, actuators and computation and control hardware. Although smart structure 
concepts may be applied to me design and implementation of buildings, dams, bridges, 
pipelines and ships, recent research efforts have been concentrated on potential 
aerospace applications in advanced aircraft, launch vehicles, and large space-based 
platforms. 

Virginia Tech has been involved in smart structures research since 1979 with sustained 
support from the NASA Langley Research Center, Simmonds Precision/Hercules and the 
Virginia Center for Innovative Technology. In 1979 and 1980, Virginia Tech participated 
in the first documented smart structure experiments, conducted at NASA Langley, which 
demonstrated the use of embedded optical fiber sensors for the measurement of strain in 
low temperature composite materials. Since then Virginia Tech has worked on the 
development of optical fiber interferometric, blackbodv, evanescent, modal domain and 
time domain sensors for the evaluation of composite cure, in-service structural 
component monitoring, nondestructive materials evaluation, and damage detection and 
evaluation. 

At the same time, Virginia Tech faculty. students and staff have cooperatively developed 
and demonstrated effective embedded actuators for structural shape modification, 
experimentally analyzed the micromechamcal problems associated with embedded sensor 
and actuator elements in advanced composites, developed computational and cortrol 
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system architecture concepts for smart structure applications, and considered the 
integration of these interdisciplinary problems in specific aerospace stiucture designs. 

This Smart Structures Workshop reviews the major non-proprietary' accomplishments in 
this area by Virginia Tech during the past several years and in particular highlights the 
recent results of the past several months More than a dozen brief technical discussions 
are followed by a working panel discussion lunch period and afternoon laboratory 
demonstrations of implemented smart structure concepts. 
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Materials Issues for Smart Structures 

Dr. Garth L Wilkes 
Professor 
Chemistry 


The term "smart structures" carr as the implication that a given material or "structure" 
when provided a particular stimulus, will generate a specific response. What the specific 
stimuli are as well as the nature and magnitude of the response will depend upon the type 
of "smart structures" desired for a given application. This brief talk will focus on material 
parameter considerations with respect to generating smart structures and will illustrate a 
number of specific cases bv example. While a somewhat higher emphasis will be given to 
those concerned with pohmeric based materials, discussion will also consider those based 
on either metallic or ceramic systems as well Finally, some new directions toward the 
development of new h\brid or complex material structures will also be outlined. 
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Experimental Analysis of Smart Structure Materials by Full-Field Optical 

Techniques 

Dr Robert Czarnek 
Assistant Professor 
Engineering Science and Mechanics 

The concept of integration of sensing and controlling devices with load carrying elements 
has been implemented in the form of so-called "smart structures”. In a smart structure, 
sensors and actuators are embedded in the material that the elements of the structure are 
made of. The materials used to make these structures are mostly fiber-reinforced 
composites. Probably the most effective type of sensors used in smart structures are made 
of optical fibers embedded between the layers of the structural fibers. The optical fiber 
deforms together with the composite element and modulates light passing through the 
fiber, allowing the detection of dangerous strain levels in the structure as well as failure 
of the material 

The presence of the optical fiber embedded in the material makes sense only if it does not 

diminish the performance of the material. Composite materials on the micromecbamcs 

level are extremely complex structures and their behavior is still not fully understood. 

The interaction between the fibers in a layer and the layers in the material causes strain 

lex els much higher than those predicted analytically. Imperfections of the manufacturing 

process complicate the situation e\en further, introducing strong nonuniformity into the 

strain distribution in the material An optical fiber introduced as an element of such a 

complicated structure brings even more unknowns into the materials analytical modeL 

* 

o 

Some preliminary experiments have been performed on graphite epoxy coupons with 
optical fibers embedded in the direction perpendicular to the direction of the load. Two 
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Iuxups. [^Ah/fiber/fh/^Ch] and [(^/^(h/fiber/DCb/Cb], were tested using high 
sensitivity moire interferometry with a sensitivity resolution on displacement 
measurements of 100 nanometers The measurements were performed on the edges of 
the specimens Strains as high as 2 % were detected in the vicinity of the fibers at a load 
equal to 50% of the calculated failure load. Figure 1 illustrates fringe patterns 
representing: 

(a) the horizontal displacement field and 

(b) the vertical displacement field at this load level in the specimen with the optical 
fiber embedded between the fibers running m the 0 degree direction. The 
concentration of high strains near the optical fiber is clearly visible. 
Corresponding patterns were recorded for a specimen with the other layup. 

The results of the performed experiments show strong strain concentration around the 
embedded sensors and a very strong dependence of the strain concentration factor on the 
layup of a composite member. The magnitude of the effect of the embedded components 
on the strain distribution in the structural elements suggests a strong need for an extensive 
experimental and analytical program on the determination of the influence of embedding 
optical sensors on the behavior of composites and on providing a means for reducing 
undesired effects. 
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Large Deployable Antennas in Space as Smart Structures 

Dr Warren L Stutzman 
Satellite Coinmumcations Group 
Professor, Electrical Engineering 

Many space-based projects are planned into the next century that require large 
antennas These antennas are most often reflectors and occasionally phased arrays. In 
j either case, the antenna is up to many meters in extent and must be deployed in space 

Performance depends critically on ne accuracy of the antenna shape. Antennas do not 
deploy accurately and also are subject to distortions while in space, primarily m the 
form of thermal gradients Therefore, sensing and control are necessary. With large 
scale antennal structures the sensing, control, and computing functions are best 
j integrated into the design, creating a smart structure. This paper will focus on the 

j large space antenna as an ideal candidate for smart structure realization. 

I 

j 

! 

I 


f 
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Active Trusses for Structural Control 

Dr Harry H Robertshaw 
Associate Professor 
Mechanical Engineering 

and 

Dr Charles F Reinholtz 
Assistant Professor 
Mechanical Engineering 

Structures that can change their configuration, termed variable geometry trusses 
(VGTs) or active trusses, are being studied for their many applications’ active control 
of structural vibration, pointing of structures or substructures, or performing robotic- 
like tasks. In fact, given a structure that can move actively, the distinction between a 
robot and a structure becomes blurred 

We are presently using two forms of active structures to control the vibration of 
continua attached to them. We have performed planar and spatial vibration reduction 
experiments with mo configurations for the acme truss. 

The planar experiments have been conducted with one plane of a tetrahedral- 
tetrahedral (T-T) truss controlling a beam that is constrained to almost planar motion. 
The spatial experiments were performed with an octahedral-octahedral (O-O) truss 
(nicknamed the geodesic truss) controlling a rod. In each case the trusses had three 
degrees of freedom embodied in three links of variable length actuated by machine 
screws driven by dc motors The motion of the links were sensed with resistive 
transducers and the motions of the beam and the rod were sensed using strain gages. 
The equations of motion for the systems were linearized about an operating point and 
an LQR optimal control law was applied with individual estimation of most of the 
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needed rate terms. Significant vibration reduction was observed with reasonable 
correlation with the analyses. 

The spatial truss was also used to perform a rudimentary robotic task: writing letters 
witn a pen placed at the end of a rigid rod attached to the structure. This robotic task 
involved solving the inverse kinematic equations for the needed link lengths in order to 
follow the required "path" to write the letters. For the geodesic truss, the solution of 
both the forward and the inverse kinematics is iterative. A reasonable amount of effort 
has gone into optimizing theses iterative procedures. 

Additionally, we have studied kinematics of most of the different possible 
configurations of the basic jnit of variable geometry trusses. The advantages and 
disadvantages of these different forms are determined based on: ease of solution of 
the kinematic equations, possible workspace of the trusses, loads involved during 
quasi-static and dynamic maneuvers, and possible lands of motions available from the 
structure. 
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Development of a Novel 'Adaptive/Smart Composite Material' Utilizing 

Nitinol Fibers 

Craig A. Rogers 
Associate Professor 
Mechanical Engineering 

A novel class of adaptive materials has been developed and demonstrated at Virginia 
Polytechnic Institute and State University. This class of adaptive materials utilize a 
shape-memory-alloy (Nitinol) in a laminated, fiber-reinforced composite. Adaptive 
materials are a relatively new class of materials that have the capability of changing 
their physical geometry', or of altering their physical properties. The basic concept 
behind the adaptive material developed at VPI&SU is that the shape memory alloys 
are integrated m a bulk material (1 e„ laminated fiber- reinforced composite, or an 
elastomer) as an actuator for force, motion, and/or variable stiffness. Possible 
applications for this class of adaptive materials are: in structures that are part of Iong- 
duration, unattended space missions (for which the material must be able to 
compensate for damage by redistributing the load around failed portions of the 
structure); in biomedical applications (such as artificial heart valves and pumps with no 
moving parts), in active vibration control of large flexible structures, in active acoustic 
control for aircraft to reduce intenor sound levels, and m robotic manipulators. 

The adaptive response of the material has been demonstrated experimentally by 
controlling the motion of a fiberglass cantilever beam. The cantilever beam is a 
unique composite material containing several shape memory alloy fibers (or films) in 
such a way that the material can be stiffened or controlled by the addition of heat (i.e., 
apply a current across the fibers) A demonstation of the capabilities of this unique 
class of materials will be "resented. Shape memory alloys and the mechanism by which 


Smart Structures Workshop Abstracts ■ Virginia Tech ■ April 29, 1988 


10 


they exhibit the characteristic shape memory effect (SME) will also be explained in 
some detail. 

The basic physical design of the demonstration adaptive/smart composite cantilever 
beam is conceptually simple. The shape memory alloy fibers are embedded in a 
material off of the neutral axis of the beam which when activated introduce a large 
force and subseque bending of the structure. Before imbedding the Fibers, the shape 
memory alloy fibers are plastically elongated and constrained from contracting to their 
normal or 'memory' length upon curing the composite material with high- temperature. 
The fibers are therefore an integral part of the composite material and the structure. 
When the Fibers are heated, generally by passing a current through the shape memory 
alloy, the fibers 'try' to contract to their 'normal' length and therefore generate a 
uniformly distributed shear load along the length of the Fibers. The shear load offset 
from the neutral axis of the structure will then cause the structure to bend in a known 
and predictable manner. There are several other configurations that are currently 
being implemented and imestigated 

This presentation will describe the actual prototypes of the novel 'Adaptive/Smart 
Composite Material' and the capabilities of the devices. The basic concepts and 
mechanisms by which the adaptive response of the material is demonstrated will also 
be presented, several of the technical problems and solutions will be discussed. Lastly, 
the future direction and impact of this technology will be postulated. 
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Fiber Optic Methods for Structure Monitoring 

Kimberly D. Bennett 

Research .Associate. Fiber tk Electro-Optics Research Center 

Electrical Engineering 

and 

Dr. Richard O. Claus 

Director, Fiber & Electro-Optics Research Center 
Professor. Electrical Engineering 

Structures engineered for integrated active control impose demanding requirements on 
system components such as actuators, sensors, and intelligent control schemes. In 
particular, sensors are required to deterrune structural status, both in terms of normal 
operating parameters such as positio. . strain, or temperature, as well as actual 
material integrity. Perhaps the most me ure approach for embedded sensing, and one 
readily applicable to attached sensors, invokes the use of optical fibers. 

Much of the pioneering work in smart structures, specifically embedded fiber optic 
sensing in structural materials, was per formed in the early 19S0’s at NASA Langley 
and Virginia Tech. Presently a number of researchers are investigating the use of fiber 
optics for monitoring ground, air, and space based structures. In addition, a number of 
different mechanisms are being exploited to perform various sensing tasks. Among 
these are optical time and frequency domain techniques, modal domain sensing, 
absolute and differential fiber interferometry, polarimetry, and various intensity 
modulation methods. Targeted observables include both point-to-point and distributed 
strain measurements, v ibration and structural mode analysis, state of cure in composite 
materials, temperature monitoring, and damage assessment. 
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Ia this paper, several general sensing methods will be reviewed. Emphasis will be 

given to those which have had particular application in smart structures research at 
Virginia Tech. 
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Modal Interference Techniques for Strain Detection in Few-Mode Optical 

Fibers 

Bradley D Duncan 

Graduate Research Assistant. Fiber & Electro-Optics Research Center 

Electrical Engineering 

Interference between the modes of an optical fiber results in specific intensity’ patterns 
which can be modulated as a function of disturbances m the optical fiber system. 
These modulation effects are a direct result of the differential phase modulation which 
arises due to the difference in propagation constants of the constituent modes. Based 
on these concepts it has been shown how the modulated intensity patterns caused by 
\anous mode groups in few-mode optical fibers (V < 5.0) can be used to detect strain. 
A detailed discussion of the modal phenomena responsible for these strain induced 
pattern modulations is presented and it is shown that strain detection sensitivities on 
the order of 10-10 can be expected for simple systems based on these modal 
interference concepts. Actual strain detection systems based on the developed theory' 
are also proposed and evaluated. 
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Applications of OTDR Methods to Structural Analysis 

Russell G. May, Jr 

Technical Director, Fiber & Electro-Optics Research Center 
Electrical Engineering 


Embedded optical fibers have been proposed as a technique for measuring strain and 
displacement in materials, particularly composite materials. By measuring parameters 
of the light exiting the embedded fiber, the integrated strain along the fiber can be 
derived. A promising technique to resolve the strain at discrete locations in a material 
is optical time domain reflectometry (OTDR). 

OTDR is a measurement technique in which a short pulse of light is launched into an 
optical fiber and is reflected back to the input end by constituent atoms along the 
length of the fiber The intensity of the reflected light is plotted as a function of the 
time delay that the light pulse suffers in returning to the input end of the fiber. Since 
the time delay is linearly related by the speed of light in the fiber to the distance from 
the input to the reflection, OTDR may be used to accurately determine the location of 
a discrete reflector. 

Experiments at Virginia Tech in 1984 demonstrated the feasibility of using OTDR to 
sense deformation of composite coupons with embedded optical fibers. In those 
experiments, OTDR was used to measure the loss induced by microbends in an 
embedded fiber at a point where pressure was applied. The spatial resolution of the 
measurement system is limited by the optical pulse width. Although the results 
demonstrated the validity of the concept, derivation of actual fiber or composite strain 
from microbend loss is difficult. 
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Further research at the Fiber and Electro-Optics Research Center in 1987 extended 
the work and enabled the direct measurement of fiber (and composite) strain by 
placing reflective markers within the embedded fiber. The markers were fabricated by 
separating the fiber at a point, and splicing the two ends together with an air gap 
between them. The Fresnel reflection from the refractive index mismatch of the splice 
yields a reflection center that produces a readily observable mark on the OTDR trace. 
The longitudinal strain of the fiber section between two air gap splices is found by 
measuring the resulting increase in time delay between the two reflective markers 
corresponding to the two air gaps 

Limitations of this technique include a limit to the number of air gap splices that may 
be placed sequentially on a fiber due to the splice loss and fiber loss budget. A limit to 
the achievable spatial resolution exists, traceable to the low energy content of very 
narrow width laser pulses. 


•* 
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Acoustic Fiber Methods for Structural Monitoring 

Dr. Ahmad Safaai-Jazi 
Associate Professor 
Electrical Engineering 

Transmission properties of cladded fiber acoustic waveguides are briefly reviewed and 
compared with those of optical fibers. Applications of acoustic fibers in sensing and 
signal processing are addressed. Attention is focused on embedded acoustic fibers for 
structural monitoring Two schemes for the implementation of acoustic fiber sensing 
concept are discussed In one scheme the acoustic fiber sensor functions in much the 
same way as an optical fiber sensor, namely that an external disturbance perturbs a 
physical parameter of the fiber yielding a change m amplitude, phase, or polarization 
of the acoustic wa\e propagating in the fiber. The amount of the change is a measure 
of the external disturbance. In a second scheme, acoustic fiber is used in a listening 
mode picking up acoustic emissions and thus monitoring the health of the structure. 
Advantages of acoustic fiber sensors in certain special applications are pointed out. 
Research and development of acoustic fiber phase sensors at the Fiber and Electro- 
Optics Research Center of Virginia Tech are addressed. 
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Acoustic Emission Detection in Smart Materials 
Kimberly D. Bennett 

Research Associate, Fiber & Electro-Optics Research Center 
Electrical Engineering 

and 

Kevin D. Zehner 

Research Assistant. Fiber & Electro-Optics Research Center 
Electrical Engineering 

The monitoring of acoustic emission (AE) is an important technique for the 
nondestructive characterization of strained materials. Time and frequency domain 
analyses of AE events yield information about the type, geometry, and location of 
defects, as well as about material properties and how failure may occur. The 
quantita:i\e interpretation of AE event signatures is critically dependent upon the 
faithfulness of the acoustic transduction and signal processing system in reproducing 
localized stress wave amplitude as a function of time. Although presently employed 
piezoelectric sensors are sensitive and reasonably broadband, they are not suited to 
embedded measurements, and are highly susceptible to electrical noise. 

This paper reports the detection of acoustic emission signals generated in symmetric 
cross-ply composite laminates using internal optical fiber sensors. Few-mode optical 
fibers were embedded between different prepreg layers during lay-up, and the resulting 
specimens were cured under pressure between heated platens. During the original 
tests, samples containing optical fibers laid m the axial direction were loaded in tension 
in specially designed grips, allowing for continuous monitoring of light. Acoustic 
signals were obtained by detecting changes in an optimized modal domain feature 
inherent in the intermodal interference pattern Independent measurements were also 
supplied by a piezoelectric transducer made by Acoustic Emission Technology 


Smart Structures Workshop Abstracts ■ Virginia Tech ■ April 29, 1988 



Corporation Agreement between transducing methods was excellent, indicating that 
acoustic emissions may be measured using fiber optic sensing techniques. 

On-going experiments imoKe samples with optical fiber laid across the direction of 
applied load, allowing for the use of conventional grips, although requiring some care 
during cure. Collection of AE signals and computational analysis is simplified with the 
use of an inexpensive digital oscilloscope and associated computer. Also, the use of 
multiple fibers for the location of AE events is investigated. Results of these tests, as 
well as limitations of the method and implications for the composites industry are 
discussed. 


* 
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Fiber Optic Assessment of Impact Damage in Smart Materials 

Michael R. Dumais 

Research Assistant, Fiber & Electro-Optics Research Center 
Electrical Engineering 

Kimberly D. Bennett, 

Research Associate, Fiber & Electro-Optics Research Center 
Electrical Engineering 

This study investigates the possibility of incorporating optical fibers into a structure 
and using them to characterize impact damage. It is believed that the ability to 
correlate signal disturbances from optical fibers with physical disturbances would 
provide invaluable real-time information about the reliability of active engineering 
systems. 

Previous work in this area, conducted at Virginia Tech with sup port from NASA 
Langley, General Motors, and the Virginia Center for Innovative Technology, involved 
the use of an embedded optical fiber mesh. Continuous wave light signals propagated 
through the fibers in the x and y directions, and were collected by an array of discrete 
photodetectors. Disturbances in the cw output due either to fiber breakage or 
intermodal interference were related to various impacts. Disadvantages arise in the 
need for injection of light into multiple fibers, and for multiple detectors. 

In order to simplify the system and quantify impact signals, an exposed optical fiber 
having a V number of approximately 4.5 was attached to a thick aluminum plate with a 
rigid epoxy. The fiber was excited with monochromatic laser light and its output 
monitored with a pbotodetector. With. the plate inverted and lying flat, several trials 
were run in which a steel mass was dropped on the plate from different heights and 
positions with respect to the fiber sensor. Using a digital waveform analyzer, fiber 
output signal changes were captured o\er the course of an impact. With the aid of a 
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computer. Fourier analysis of the signals was performed, and attempts were made to 
correlate spectral signatures with different types of impacts. Also, output signal power 
densities were calculated and analyzed with respect to impact amplitude and location. 
After attaching a second fiber to the plate, similar drop tests were run and blind impact 
location information was extracted by examining the time lag between the disturbance 
signals from the two fibers 

In this paper the limitations of fiber optic impact and damage detection are discussed, 
and direction for future efforts proposed. 
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Control and Signal Processing Design Issues in Smart Structures 

John C. McKeeman 
Assistant Professor 
Electrical Engineering 

and 

Doug K. Lindner 
Assistant Professor 
Electrical Engineering 

The next generation of space structures will be significantly different from the current 
generation of spacecraft. The new space structures will be made of new composite 
materials which will allow them to be of significantly larger dimensions The increased 
size, in particular, along with tighter performance specifications will require that the 
structure be des.gned and built to function as a cohesive unit in which subassemblies 
strongly interact to achieve sxstem performance specifications. The system wide 
integration wall be achieved by control loops which (effectively) coordinate the 
subassemblies. Such a structure we call a smart structure. 

In this paper we examine the system configuration of smart structures from a systems 
perspective. To focus the discussion we consider as an example structure a large space 
antenna, the Hoop Column Antenna We describe the control loops required to 
implement and maintain the antenna in space. Special attention is gixen to the 
components required to implement the control system. Sex eral control design issues 
are highlighted. 

To successfully produce this control system, several new computing technologies must 
be tailored to the applications required by each subassembly. At the highest level in 
this hierarchical control system, .Artificial Intelligence coupled with Neural Networks 
will control the housekeeping actixities and gross structure control. At the lowest level. 
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application specific integrated circuits will be customized to handle the high 
throughput, high bandwidth requirements for the sensors and actuators. This paper 
focuses on new computers architectures required to implement advanced control 
structure within a intelligent structure. 




Smart Structures Workshop Abstracts ■ Virginia Tech ■ April 29, 1988 



APPENDIX 


- Abstracts from Papers To Be Presented at The SPIE 
Conference (Boston, MA) - September 1988 



TABLE OF CONTENTS 


Interferometric Measurements of Strain Concentration Induced 

by Fiber Embedded in Fiber Reinforced Composites 1 


Demonstration of A Smart Structure with Embedded Actuators 

and Sensors For Active Control 3 

Full Field Analysis of Modal Domain Sensor Signals 

for Structural Control 5 

Intensity Pattern Modulation in Optical Fiber Modal Domain 

Sensor Systems: Experimental Results 7 

Smart Structures Program at Virginia Tech 9 


Fiber Optic Sensor Systems for Smart Aerospace Structures . . 


.11 



r;.c*3 s 

CF r'GQK CUHJDJ 

INTERFEROMETRIC MEASUREMENTS OF STRAIN CONCENTRATIONS INDUCED BY 
AN OPTICAL FIBER EMBEDDED IN FIBER REINFORCED COMPOSITES 

R. Czarnek*, Y.F. Guo*, K.D. Bennett 4 and R.O. Claus T 

* Department of Engineering Science and Mechanics 
+ Department of Electrical Engineering 
Virginia Polytechnic and state University 
Blacksburg, Virginia 24061 

The concept of using of optical fibers for sensing the 
mechanical response of composite structures to applied loads has 
been implemented in the form of so-called "smart skins." An 
optical fiber imbedded between the layers of the structural fib- 
ers deforms together with the composite structure and modulates 
the light passing through the fiber. This allows for the detec- 
tion of dangerous strain levels m -he structure, as well as 
failure of the material. 

Tne presence of the cubical finer embedded m the material 
is viable only if it does not diminish the performance of the 
material. Composite material cn tr.e r icrcnechanics level is an 
extremely complex structure ar.d its behavior is still not fully 
understood. The interaction between the composite fibers in a 
layer and the layers m the material causes strain levels much 
higher than these predicted analytically. Imperfections of the 
manufacturing process cc-plicate the situation even further, 
introducing strong ncnar.ifcr~.ty into the strain distribution of 
the material. 

An optical fiber i -traduced as an element of such a compli- 
cated structure introduces more unknowns into the analytical 
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model of the material. Several important questions can and 
should be answered before the concept of embedding optical fibers 
can be commercialized. These include understanding how fibers 
influence the strain field in a composite element, and if the 
strain concentration factors introduced are significant when 
compared to the natural strain variations which arise from the 
nanufacturing process. Strain concentration dependence on the 
layup and position of the optical fiber is also unknown. 
Finally, understanding of whether the inclusion of optical fiber 
can cause premature failure of a composite element, and how to 
reduce any such effects is necessary. 

These is ues require serious experimental investigation. 
Moire interferoretry is a high sensitivity method yielding full 
field measurement of m-plane displacements . Its sensitivity is 
limited to 1/2 cf a wavelength of the light used per fringe 
order, and m practical applications a sensitivity resolution cf 
the order of IOC nm is standard m the VPI1SU Photomechanics Lab- 
oratory. The spatial resolution is limited only by the quality 
of the optical system used to capture the image and details as 
small as 1 micrometer can he successfully resolved. 

In this paper, ve report moire interferometry experiments 
performed cm. tr.e face cf a cc-posite, as well as cm. its edge. 
Measurements on the face indicate tr.e glooal influence of tr.e 
fiber cn the deformation of the specimen, and those on its edge 
give information about the very localized u-train disturbance i" 
the direct vicinity of the f.uer Represer native layup sequences 
have been investigated, e~d tr-: effect cf a— redding optical fiber 
has been analyzed. 
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DEMONSTRATION OF A SMART STRUCTURE UITH EMBEDDED 
ACTUATORS AND SENSORS FOR ACTIVE CONTROL 

C. A. Rogers*, K.D. Bennett + , H.H. Robertshav*, and J.C. McKeenan 4- 

* Department of Mechanical Engineering 
+ Department of Electrical Engineering 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 25061 

A novel class of "smart" materials having imbedded sensors 
and actuators and appropriate "intelligence" has been developed 
and demonstrated at Virginia Tech. This class of adaptive mate- 
rials utilize a shape rercry alloy (Nitinol) in a laminated, fib- 
er-reinforced composite as the embedded distributed actuator. 
The basic concept behind tne smart materials developed at Virgi- 
nia Tech is that the shape memory alley fibers are integrated 
into a bulk material such as a composite structural member, and 
used as actuators for force, notion, and, or variable stiffness. 
Sensing of material dynamics and strain is performed by embedded 
optical fiber sensors. Possible applications for this class of 
adaptive materials are im. structures that are part of lom.g- 
duration, unattended missions (for which the material rust be 
able to compensate for damage by redistributing the load around 
failed portions of the structure) , in biomedical applications 
(suer, as artificial heart valves and pumps vith no moving parts) , 
m active vibration control of large flexible structures, in 
acti -e acoustic control for aircraft to reduce interior sound 
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levels, and in robotic manipulators . 

The adaptive response of structural "members has been demon- 
strated in hardware by controlling the motion of a fiberglass 
cantilever beam. The bean" is a unique composite material con- 
taining several shape memory alloy fibers m such a way that the 
material can be stiffened or controlled by the addition of heat, 
created by applying a current across the fibers. Embedded opti- 
cal fibers monitor the notion of the cantilever bean and applica- 
tion-specific intelligent hardware has been designed and built to 
control the notion and/or vibration of the structure. 

Results from experimental studies showing active motion 
control with this novel smart structure will be presented. Shape 
memory alloys and the mechanism, by which they exhibit the charac- 
teristic shape memory effect (Sill) fcr use as an actuator will be 
explained, as well as the physical design and implementation of 
the embedded strain sensing fiber optics. The control strategy 
(i.e. "intelligence) will also be presented. Lastly, the future 
direction and impact of this technology will be postulated. 
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FULL FIELD ANALYSIS OF MODAL DOMAIN SENSOR 
SIGNALS FOR STRUCTURAL CONTROL 


K.D. Bennett, J.C. KcKeeman, and R.G. May 

Fiber and Electro-Optics Research Center 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24061 

when coherent light is launched into multimode optical 
fiber, a complex interference of the modes occurs’oalong the fiber 
length, and a speckle pattern results at the output. The exact 
arrangement of the speckles depends on the modal arrangement 
integrated along the fiber, and is modulated by local changes in 
geometry and refractive index profile. Interrogating the modal 
domain signal of an optical fiber embedded or attached to a 
structure is a well established means of determining pseudo- 
static strain and structural vibration frequency and mode shape, 
as veil as monitoring acoustic e*. arts occurring in the material. 

Prior research has depended upon spatially filtering the 
optical output to determine the speckle motion. In addition to 
requiring determination of an optical Q point, this method also 
recessanly results in a less cf optical pover and useful infer- 
rable-. Advancement m never collection has been, made through 
the use of lo” rode optical ficers, which share pev-or among, fe*- 
output speckles. Sc-o novel ir.-li'u fiber techniques have also 
recent'./ boor repoitia for si-plif/irg the process and gathrrin. 


ORIGfNAE PAGE 13 
OF POOR QUALITY 


5 



In this paper, we report the use of a detector array to 
process the modal domain output patterns. ’ The array is linked ro 
an IBM-PC, which reduces the data in real time. Full field det- 
ermination of low mode speckle rotation due to axial strain is 
noted, as well as gross speckle rotation ir. multimode fibers due 
to bulk fiber notion. Statistics of speckle notion in both low 
an multimode fibers has also been investigated. Finally, the 
importance of this method for use in smart structure dynamics and 
control will be discussed. 
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INTENSITY PATTERN MODULATION IN OPTICAL FIBER 
MODAL DOMAIN SENSOR SYSTEMS: EXPERIMENTAL RESULTS 

Bradley D. Duncan 

Fiber & Electro-Optics Research Center 
Virginia Polytechnic Institute and State University 
Blacksburg, Va 24061 

Brian W. Brennan 

Sinnonds Precision I.S.D. 

Panton Road, Vergennes, VT 05491 

(ABSTRACT) 

The past several years have given rise to nany optical fiber 
sensing techniques which ray be erployed in systems to detect 
'Ch observables as temperature, pressure, displacement and 
strain. As sensing devices, optical fibers are attractive due to 
their snail size, flexibility, ruggedness and intrinsic 
dielectric nature. In addition, their proclivity to be highly 
sensitive to external perturbations is quickly making optical 
floors the sensing devices of choice for rar.y passive 
rcr-intrusive sensor applications. 

Sirply put, optical fiber sensors exploit the effects of 
perturbations and external factors on the light in the fiber. 
Unparalleled in sensitivity are the phase modulated sensors, a 
classic exarple of which is the dual-fiber Mach-Zehnder 
_nt erf ercr oter . By bei"a based on an srsolute phase modulation 
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.. 4 e me, however,, such a sensor tends to be quite complex, with 
the reference phase and. quadrature point often needing to be 
actively controled. 

By comparison, systems utilizing single-fiber modal 
interference schemes, such as these to be discussed in this 
paper, essentially cause the two airs of the interferometer to be 
within the same environment. This provides for a high common 
node rejection capacity while affording the system the further 
advantages of simplicity and increased ruggedness, with no need 
for couplers. Though remaining highly sensitive, as will be 
shown, some sensitivity will be sacrificed in a single-fiber 
interferometer, or modal domain sensor, due to the fact that the 
sensor mechanism will now be based on differential phase 
modulation. 

In a companion paper a detailed theoretical analysis of the 
nodal phenomena responsible for ihe far leld intensity pattern 
-'cdulaticn of various nodal domain ser.so: systems is presented. 
This paper v/ill present experimental results which verify the 
expected system behavior and whicn seme as proof of the 
significance of the modal domain sensing technique. The results 
of both static and dynamic strain measurements are presented, as 
are data unich demonstrate the strain sensitivities achievable 
ith modal domain sensor systems. System design, stability, and 
ml it rat ion techniques are also presorted. 
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SMART STRUCTURES PROGRAM AT VIRGINIA TECH 

R. O. Claus and K. D. Bennett, Fiber and Electro-Optics Research 
Center, Virginia Tech, Blacksburg, VA 24061. 

Abstract 

This paper reviews the smart structures and avionics research 
and teaching program at Virginia Tech. Started in 1979 with 
support from the NASA Langley Research Center, this program has 
grown to include interdisciplinary participation between more 
than twenty faculty from seven departments in the Colleges of 
Engineering and Arts and Sciences, graduate and undergraduate 
engineering classes which specifically address smart structures 
issues, and directed research programs conducted in cooperation 
• ith more than a dozen aerospace companies and government 
organizations . 

Current smart structures research includes major efforts in 
1) the development of embedded and attached optical fiber and 
acoustic fiber sensors for a wide range of applications including 
cure ronitoring, in-service lifetime structural monitoring, 
nondestructive evaluation, and impact and damage detection and 
analysis, 2) sensor signal multiplexing, processing and data 
handling to achieve near real time districuted structural 
analysis, and 3) the combined use of embedded or attached sensors 
and actuators to achieve controlled structural - response. Special 
campus facilities which have been used for this work include an 
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structure fabrication and curing, and laboratories for optical 
fiber sensor development, materials response and nondestructive 
evaluation, structural control testing, and computer engineering. 
The work in all of these areas during the past few years by 
different Virginia Tech faculty, staff, graduate student and 
visiting scientist groups will be reviewed. 

[This work has been supported in particular by sustained 


funding from the NASA Langley Research Center, Simmonds 
Precision/Hercules, and the Virginia Center for Innovative 
Technology. ] 







Fiber Optic Sensor Systeas for 
Smart Aerospace Structures 

W. J. Rowe, R. o. Claus*, D. A- Martin+, and D. Tritesl 

Lockheed Aeronautical Systems Company 
Dept. 72-51, Zone 303 
Marietta, Georgia 30063 


ABSTRACT 

Current and projected advances in aerospace structures 
technology for such ongoirg programs as the Advanced Tactical 
Fighter, National Aerospace Plane, and NASA/DoD space platforms 
are used as the starting point for a top-leval structural 
monitoring systems analysis. Performance requirements are 
considered for such structures, and used to identify needs for 
integral sensors and systems. Results of recent work are 
discussed, which have demonstrated the potential for fiber optic 
sensor systems which can assess and report cr. structural 
condition in real time. The essential system elements are 
defined, and constraints imposed by venicle configurations and 
mission performance are considered. The resulting development 
needs, for fiber optic sensors and systems to achieve Smart 
Structures in these applications, are stated in conclusion. 


♦Fiber & Electro-Cptics Research Center, 648 Whittemore Hall, 
Virginia Tech, Blacksburg, VA 24061. 

+Kelly Johnson Research & Development Center, Rye Canyon, 
Lockneed Aeronautical Systems Cor.cany, P. 0. Box 551 
(7471/211/2), Buroank, CA 91520.' 

#Bldg. 104, Org. 62-13, Lockheed Missiles & Space Company, 
Space Systems Division, 1111 Lockheed Way, Sunnyvale, CA 
94036. 
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"Optical Fiber Fundamentals," Proc. Optical Fiber System 
Design (Madison, WI) February 1988 (R.O. Claus) . 

"Smart Skins and Avionics," LTV Aerospace (Dallas, TX) 
Februaury 1988 (R.O. Claus) . 

"Optical Fiber Sensors" (Short Course), MFOC *88 (Crystal 
City, VA) March 1988 (R.O. Claus) . 
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